<div dir="ltr"><div><div>Hi Ipshita,<br></div>I think you need to add some additional resources to your res2 resource list to get stippling to be applied correctly. I took the graphics portion of your script and read in some data I had here to test this. Try the following:<br><br> wks = gsn_open_wks("png","June_sst_jjas_cor_prob") <br> gsn_define_colormap(wks,"BlWhRe") ; choose colormap<br> <br> res = True ; make plot mods<br> res@gsnDraw = False<br> res@gsnFrame = False<br> res@cnFillOn = True ; turn on color<br> res@gsnSpreadColors = True ; use full colormap<br> res@lbLabelAutoStride = True ; automatic lb label stride<br><br> res@cnLinesOn = False ; turn off contour lines<br> res@cnLevelSelectionMode = "ManualLevels" ; manually set cnlev<br> res@cnMinLevelValF = -1. ; min level<br> res@cnMaxLevelValF = 1. ; max level<br> res@cnLevelSpacingF = .1 ; contour level spacing<br><br> lag = 0<br> res@tiMainString = "June SST with JJAS rainfall"<br> plotA = gsn_csm_contour_map_ce(wks,corr_june(:,:,0),res)<br><br> ; Significance plot on top of correlation plot<br> ;========================= PLOT 2 =========<br> res2 = True ; res2 probability plots<br><br> res2@gsnDraw = False ; Do not draw plot<br> res2@gsnFrame = False ; Do not advance frome<br> res2@cnFillOn = True ; turn on color/pattern fill<br> res2@cnMonoFillPattern = False ; allow different patterns<br> res2@cnMonoFillColor = True ; only use one color (black)<br><br> res2@cnLevelSelectionMode = "ExplicitLevels" ; set explicit cnlev<br> res2@cnLevels = (/.95/) ; only have 1 contour level<br> res2@cnFillPatterns = (/-1,17/) ; don't fill <0.95, stipple >=0.95<br> <br> res2@gsnAddCyclic = True ; may or may not be needed<br> <br> res2@cnInfoLabelOn = False ; turn off info label<br><br> res2@cnLinesOn = False ; do not draw contour lines<br> res2@cnLineLabelsOn = False ; do not draw contour labels<br><br> res2@cnFillScaleF = 0.5 ; add extra density<br><br> printMinMax(prob_june,0) ; check to make sure values are<br></div> ; between 0 and 1.<br><div><br> plot2 = gsn_csm_contour(wks,prob_june, res2)<br> overlay (plotA, plot2)<br><br> draw (plotA) <br> frame(wks)<br>;-----------------------------------------<br><br></div><div>A few notes about what I added:<br></div><div>1 - I added gsnDraw/gsnFrame = False to both resource lists as you are overlaying the two plots, and added a frame(wks) call at the end.<br><br></div><div>2 - cnFillOn must = True when color or pattern filling.<br></div><div>3 - If you are using more than one fill pattern you must set cnMonoFillPattern = False.<br><br></div><div>For more pattern fill examples see examples 10 and 11 here:<br><a href="http://www.ncl.ucar.edu/Applications/coneff.shtml#ex10">http://www.ncl.ucar.edu/Applications/coneff.shtml#ex10</a><br></div><div><br></div><div>Hopefully that all makes sense.. If you have any further questions please respond to the ncl-talk email list.<br></div><div>Adam<br><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jul 21, 2015 at 10:17 PM, Ipshita Majhi <span dir="ltr"><<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Dear NCL,<br><br></div>I don't get any errors when I run this code, but the stippling does not show up. Could you guide me on this.<br></div>Best Regards<br></div>Ipshita<br><div><div><div><div><br><br>;This is for correlating JJAS monthly precp with JJAS temperature<br><br>;******************************************************<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br>;*******************************************************<br><br>;=============================================<br>; Reading data<br>;=============================================<br><br>a=addfile("~/Documents/SST/Monthly/<a href="http://sst.mnmean.v4.nc" target="_blank">sst.mnmean.v4.nc</a>","r")<br>JJAS =asciiread("~/Documents/SST/Monthly/jjas_1871_2012.txt",(/142,1/),"float")<br><br><br>;=============================================<br>;The dataset goes from 1854-2015<br>sst=a->sst<br>; copy meta data and coordinate variables using contributed functions<br> copy_VarAtts(a,sst) <br> copy_VarCoords_1(a,sst)<br>;============================================<br><br>;Extracting 1871-2012 for SST<br> sst_1871_2012=sst(216:1908,:,:)<br><br>copy_VarAtts(sst,sst_1891_2012) <br>copy_VarCoords_1(a,sst_1891_2012)<br> <br>;============================================<br>;Extracting 1871 to 2012 in JJAS rainfall<br><br><br>June = new((/142,89,180/),float)<br>June@_FillValue = -9.9621e+36<br><br>Feb = new((/142,89,180/),float)<br>Feb@_FillValue = -9.9621e+36<br>;==========================================<br><br>;Extracting SST for month of June<br><br>do nyr=0,1692,12<br><br>June(nyr/12,:,:)=sst_1871_2012(nyr+5,:,:) <br>Feb(nyr/12,:,:)=sst_1871_2012(nyr+1,:,:) <br>end do<br><br>copy_VarAtts(sst_1871_2012,June) <br>copy_VarCoords_1(sst_1871_2012,June)<br><br>copy_VarAtts(sst_1871_2012,Feb) <br>copy_VarCoords_1(sst_1871_2012,Feb)<br><br>;==========================================<br>;Calculating correlation<br><br><br>e=transpose(JJAS)<br><br>print(dimsizes(JJAS))<br>June_reorder=June(lat|:,lon|:,time|:)<br><br>copy_VarAtts(June,June_reorder) <br>copy_VarCoords_1(June,June_reorder) <br><br>corr_june=escorc(June_reorder,e)<br><br>copy_VarAtts(June_reorder,corr_june) <br>copy_VarCoords_1(June_reorder,corr_june)<br><br><br>prob_june = rtest(corr_june(:,:,0), 142,0)<br>copy_VarAtts(June_reorder,prob_june) <br>copy_VarCoords_1(June_reorder,prob_june)<br>;===========================================<br><br>; Now plotting the correlation<br><br><br> wks = gsn_open_wks("pdf","June_sst_jjas_cor_prob") ; open a ps plot<br> gsn_define_colormap(wks,"BlWhRe") ; choose colormap<br> <br> res = True ; make plot mods<br><br> res@cnFillOn = True ; turn on color<br><br>;---These two resources not needed in NCL V6.1.0<br> res@gsnSpreadColors = True ; use full colormap<br><br>;---This resource defaults to True in NCL V6.1.0<br> res@lbLabelAutoStride = True ; automatic lb label stride<br><br> res@cnLinesOn = False ; turn off contour lines<br> res@cnLevelSelectionMode = "ManualLevels" ; manually set cn levels<br> res@cnMinLevelValF = -1. ; min level<br> res@cnMaxLevelValF = 1. ; max level<br> res@cnLevelSpacingF = .1 ; contour level spacing<br><br> lag = 0<br> res@tiMainString = "June SST with JJAS rainfall"<br> plotA = gsn_csm_contour_map_ce(wks,corr_june(:,:,lag),res) <br> <br> ; Significance plot on top of correlation plot<br> ;<br> <br> ; Significance plot on top of correlation plot<br> ;========================= PLOT 2 ==============================<br> res2 = True ; res2 probability plots<br><br> res2@gsnDraw = False ; Do not draw plot<br> res2@gsnFrame = False ; Do not advance frome<br><br> res2@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels<br> res2@cnMinLevelValF = 0.00 ; set min contour level<br> res2@cnMaxLevelValF = 1.05 ; set max contour level<br> res2@cnLevelSpacingF = 0.05 ; set contour spacing<br><br> res2@cnInfoLabelOn = False ; turn off info label<br><br> res2@cnLinesOn = False ; do not draw contour lines<br> res2@cnLineLabelsOn = False ; do not draw contour labels<br><br> res2@cnFillScaleF = 0.6 ; add extra density<br><br><br> plot2 = gsn_csm_contour(wks,prob_june, res2) <br> ;plot2 = ShadeLtContour(plot2, 0.07, 17) ; shade all areas less than the<br> ; 0.05 contour level<br> overlay (plotA, plot2)<br><br> draw (plotA)<br>;***************************************************************<span class="HOEnZb"><font color="#888888"><br clear="all"><br>-- <br><div><div dir="ltr"><div><div><div><div>Ipshita Majhi<br></div>PhD Candidate<br></div>University of Alaska , Fairbanks<br></div>Atmospheric Science Department<br></div><a href="tel:%28907%29978-4220" value="+19079784220" target="_blank">(907)978-4220</a> <a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a><br></div></div>
</font></span></div></div></div></div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a> </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</div>