<div dir="ltr"><div><div><div><div>Hi, <br></div>I checked the output numbers. Those values are  fine and ranging from 0-4. But the plot is empty. I would appreciate any help and suggestions.<br><br></div>Thanks<br></div>Dz<br></div>;;;;;;;;;;;Plot properties;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br><div>   wks  = gsn_open_wks("x11",Out)       ; specifies a ps plot<br>gsn_define_colormap(wks,""+col+"")<br> plot = new(1,graphic)<br>   res                       = True<br>res@gsnDraw             = False<br>   res@gsnMaximize           = True             ; make large<br>   res@gsnFrame           =  False            ; make large<br>   res@cnFillOn              = True             ; turn on color<br>   res@cnLinesOn             = False            ; turn off contour lines<br>   res@cnLineLabelsOn        = False            ; turn off contour line labels<br>   res@cnInfoLabelOn       = False<br>res@lbLabelBarOn        = False           ; turn off individual cb's<br><br>res@gsnTickMarksOn = False<br>res@mpGridAndLimbOn       = False<br><br>   res@cnLevelSelectionMode  = "ManualLevels"   ; set manual contour levels<br>   res@cnMinLevelValF        =  0            ; set min contour level<br>   res@cnMaxLevelValF        =   4           ; set max contour level<br>   res@cnLevelSpacingF       =   0.5          ; set contour interval<br>   res@gsnAddCyclic          =  False<br>   res@mpFillOn              = False            ; turn off default background gray<br> res@mpLimitMode         = "Corners"<br><br>  res@mpRelativeCenterLon = True<br><br>  res@mpCenterLonF        = -152.2044<br><br>  res@mpRelativeCenterLat = True<br><br>  res@mpCenterLatF        = 90.<br><br>  res@mpRightCornerLonF   = -103.99817<br><br>  res@mpRightCornerLatF   = 67.58854<br><br>;  res@mpLeftCornerLonF    = -181.7859<br>  res@mpLeftCornerLonF    = -159.0859<br><br>  res@mpLeftCornerLatF    = 43.23277<br>res@mpGeophysicalLineThicknessF = 3.0<br>res@mpGeophysicalLineColor = "Black"; (/22/)<br>;res@mpNationalLineThicknessF  = 3.0<br>res@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; add state boundaries<br>res@mpNationalLineColor  = res@mpGeophysicalLineColor<br>res@mpUSStateLineThicknessF = 3.0<br>res@mpUSStateLineColor  = res@mpGeophysicalLineColor<br>res@gsnStringFont   = "helvetica-bold"<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>; res@gsnCenterString       = year(0)+"-"+year(nyrs-1) <br>  res@gsnLeftString   = ""               ; add the gsn titles<br>  res@gsnRightString  = ""<br>;   res@tiMainString          = "Trend, Reanalysis-2m Temperature(~S~o~N~C)"    ; fili<br>;res@tiMainFontHeightF  = 0.02<br>res@gsnCenterStringFontHeightF     = 0.020<br>res@gsnCenterString     = "A-Clim-"+mon+""+yr+""<br><br>plot = gsn_csm_contour_map_ce(wks,zAvg,res)<br><br>;************************************************<br>; create panel<br>;************************************************<br> resP                  = True                   ; modify the panel plot<br> resP@gsnPanelLabelBar = True<br> resP@lbLabelFontHeightF = 0.020                ; set font height of Label Bar labels<br> resP@gsnPanelBottom   = 0.2                    ; shrink panel plot by setting bottom edge of plot<br> resP@gsnPanelTop      = 0.9                    ; shrink panel plot by setting top edge of plot<br> resP@pmLabelBarWidthF  = 0.55                 ; make thinner<br>;  resP@pmLabelBarHeightF = 0.05<br> resP@pmLabelBarHeightF = 0.1<br> resP@gsnPanelYWhiteSpacePercent = 5.           ; increase spacing along Y-axis between panel plots<br> resP@cnLinesOn       = False                   ; turn off contour lines<br>;  resP@lbTitleString    = "2m-temperature(~S~o~N~C)"<br>  resP@lbTitleString    = ""+unit+""<br>  resP@lbLabelFont    = "helvetica-bold"          ; add common colorbar<br>  resP@lbTitleFont    = "helvetica-bold"          ; add common colorbar<br>  resP@lbTitleFontHeightF= .025                 ; make title smaller<br>   resP@cnLineLabelsOn  = False<br><br>  resP@txFontHeightF   = 0.03<br>  resP@txFont   = "helvetica-bold"<br>  ;resP@txString   = ""+Ind+"-"+Initial+""<br>  gsn_panel(wks,plot,(/1,1/),resP)               ; now draw as one pl<br>end<br></div></div>