<div dir="ltr"><div><div><div>Dear NCL,<br><br></div>I don&#39;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 &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;<br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl&quot;<br>;*******************************************************<br><br>;=============================================<br>; Reading data<br>;=============================================<br><br>a=addfile(&quot;~/Documents/SST/Monthly/<a href="http://sst.mnmean.v4.nc">sst.mnmean.v4.nc</a>&quot;,&quot;r&quot;)<br>JJAS =asciiread(&quot;~/Documents/SST/Monthly/jjas_1871_2012.txt&quot;,(/142,1/),&quot;float&quot;)<br><br><br>;=============================================<br>;The dataset goes from 1854-2015<br>sst=a-&gt;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(&quot;pdf&quot;,&quot;June_sst_jjas_cor_prob&quot;)              ; open a ps plot<br>  gsn_define_colormap(wks,&quot;BlWhRe&quot;)              ; 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 = &quot;ManualLevels&quot;      ; 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         = &quot;June SST with JJAS rainfall&quot;<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 = &quot;ManualLevels&quot; ; 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>;***************************************************************<br clear="all"><br>-- <br><div class="gmail_signature"><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>(907)978-4220 <a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a><br></div></div>
</div></div></div></div></div>