[ncl-talk] Gap in a plot
Anahita Amiri Farahani
aamir003 at ucr.edu
Tue May 22 16:19:13 MDT 2018
Hi all,
I plotted change in the low cloud as a difference between two simulations
and added the 95% confidence level but the problem is there is a band that
significant values are not shown. it is around 180 degrees longitude.
Attached is the figure, and here is the code that I used:
wks = gsn_open_wks("eps","low_cloud_summer")
gsn_define_colormap(wks,"BlWhRe") ; choose color
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
res = True ; plot mods desired
res at cnFillOn = True ; turn on color
res at gsnSpreadColors = True ; use full color map
res at cnLinesOn = False ; no contour lines
res at cnLineLabelsOn = False ; no line labels
;res at gsnAddCyclic = False
res at gsnDraw = False ; don't draw
res at gsnFrame = False ; don't advance frame
res at cnInfoLabelOn = False ; turn off cn info label
res at cnFillOn = True ; turn on color
res at gsnSpreadColors = True ; spread out color table
res at gsnSpreadColorStart = 2 ; color to startt
res at lbLabelBarOn = False ; turn off individual cb's
res at cnLevelSelectionMode = "ManualLevels" ; AutomaticLevels
res at cnMinLevelValF = -20. ; min level
res at cnMaxLevelValF = 20. ; max level
res at cnLevelSpacingF = 2 ; interval
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
res at tiMainString = " Summer changes in low clouds" ;Cloud Fraction
Derivative fall
plot= gsn_csm_contour_map(wks,low_summer,res)
;delete(wks)
; ========================= PLOT 2 ==============================
res2 = True ; res2 probability plots
res2 at gsnDraw = False ; Do not draw plot
res2 at gsnFrame = False ; Do not advance frame
res2 at cnInfoLabelOn = False ; turn off info label
res2 at cnLinesOn = False ; do not draw contour lines
res2 at cnLineLabelsOn = False ; do not draw contour labels
res2 at cnFillScaleF = 0.6 ; add extra density
;; delete(alpha at long_name)
;; delete(alpha at units)
plot2 = gsn_csm_contour(wks,alphat_low_summer, res2)
plot2 = ShadeGtContour(plot2, 0.05, 17) ; shade all areas less than the
; 0.05 contour level
overlay (plot, plot2)
draw (plot)
frame(wks)
delete(wks)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180522/c08e29b1/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: temp_summer700.eps
Type: application/postscript
Size: 830952 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180522/c08e29b1/attachment-0001.eps>
More information about the ncl-talk
mailing list