;************************************************ load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl" load "$HAO_NCL_LIB/halfzp2full.ncl" load "$HAO_NCL_LIB/axis.ncl" ;************************************************ begin ;************************************************ ; before the data can be plotted, there is a ; sequence of things that have to be done to it: ;************************************************ in = addfile("/chewy/d/murphys/files/pyrj001.nc","r") x = in->UN x = halfzp("UN",x) ; convert to full zp t2plt = 2 ; time index to plot lat2plt = -43 ; lat value to plot ;************************************************ ; create plot ;************************************************ wks = gsn_open_wks("ps" ,"zplon") ; open a ps file gsn_define_colormap(wks,"BlWhRe") res = True ; plot mods desired res@cnLineLabelPlacementMode = "constant" ; best line label mode res@cnFillOn = True ; add color res@gsnSpreadColors = True ; use full range of color map res@cnLineLabelsOn = True ; contour line labels res@lbLabelAutoStride = True res@tiMainString = "UT = "+in->ut(t2plt)+" Lat = "+lat2plt ; adjust the aspect ratio of plot. with vpXF, we slide the plot over a little ; to make room for the expanded x dimension. res@vpWidthF = 0.8 res@vpHeightF = 0.4 res@vpXF = 0.1 ; manually select contour levels res@cnLevelSelectionMode = "ManualLevels" ; set manual cn levels res@cnMinLevelValF = -100. ; set min contour level res@cnMaxLevelValF = 100. ; set max contour level res@cnLevelSpacingF = 20. ; set contour spacing plot = gsn_csm_contour(wks,x(t2plt,:,{lat2plt},:), res) end