<div dir="ltr"><div class="gmail_default" style="font-size:small">It&#39;s hard for us to look at a code with lots of calculations to try to figure out why a plot is not correct, if we don&#39;t have access to the data.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Have you tried debugging your code, to see if any of the values you&#39;ve calculated are incorrect?</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I suggest using printVarSummary and printMinMax to make sure you have the values and metadata that you expect:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">printVarSummary(corr_dec)</div><div class="gmail_default">printMinMax(corr_dec,0)</div><div class="gmail_default">printVarSummary(prob_dec)</div><div class="gmail_default" style="font-size:small">printMinMax(prob_dec,0)</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">etcetera.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">As a side, why do you have this line:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px">e=transpose(transpose(JJAS_</span><span style="font-size:12.8px">1967_2012))</span><br style="font-size:12.8px"></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px"><br></span></div><div class="gmail_default"><span style="font-size:12.8px">If you transpose an array and then transpose it again, you are going to get the same array back.  This is an unnecessary operation.</span></div><div class="gmail_default"><span style="font-size:12.8px"><br></span></div><div class="gmail_default"><span style="font-size:12.8px">--Mary</span></div><div class="gmail_default"><span style="font-size:12.8px"><br></span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Feb 2, 2016 at 11:35 AM, Ipshita Majhi <span dir="ltr">&lt;<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>&gt;</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 calculated correlation and am trying to plot the significance plot. The significance plot , is not plotting correctly, it is giving no significance where there is high correlation and vice versa. It will be great if you could point out where I am going wrong.<br></div>Best Regards<br></div>Ipshita<br clear="all"><div><div><div><div><br>;=============================================<br>; Reading data<br>;=============================================<br><br>a=addfile(&quot;~/Documents/SST/Monthly/<a href="http://sst.mnmean.v4.nc" target="_blank">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>;Extracting JJAS for 1967 to 2012<br>;============================================<br><br>JJAS_1967_2012=JJAS(96:141,0)<br><br>e=transpose(transpose(JJAS_1967_2012))<br><br>;=============================================<br>;The dataset goes from 1854-2015<br><br>sst=a-&gt;sst<br><br> copy_VarAtts(a,sst)   <br> copy_VarCoords_1(a,sst)<br>;============================================<br><br>;Extracting sea surface temperature<br><br>sst=a-&gt;sst<br><br>;Extracting 1967-2012 for SST<br>sst_1967_2012=sst(1356:1907,:,:)<br><br>copy_VarAtts(sst,sst_1967_2012)   <br>copy_VarCoords_1(a,sst_1967_2012)<br>   <br>;============================================<br>;Extracting 1871 to 2012 SST<br><br>Dec_sst = new((/46,89,180/),float)<br>Dec_sst@_FillValue = -9.9621e+36<br><br><br>;==========================================<br>;Extracting SST for each month<br>;========================================<br><br>do nyr=0,551,12<br><br>Dec_sst(nyr/12,:,:) =sst_1967_2012(nyr+11,:,:)<br>end do<br><br>copy_VarAtts(sst_1967_2012,Dec_sst) <br>copy_VarCoords_1(sst_1967_2012,Dec_sst)<br><br>;=======================================<br>;Reordering<br>;=======================================<br><br>Dec_reorder=Dec_sst(lat|:,lon|:,time|:)<br><br>Dec_reorder@_FillValue = -9.9621e+36<br><br>copy_VarAtts(Dec_sst,Dec_reorder) <br>copy_VarCoords_1(Dec_sst,Dec_reorder)<br><br>;=====================================<br><br>;======================================<br>;Calculating Correlation<br>;======================================<br><br>corr_dec=escorc(Dec_reorder,e)<br>copy_VarCoords_1(Dec_reorder,corr_dec)<br><br>;===========================================<br><br>;Calculating significance<br>prob_dec   = rtest(corr_dec(:,:,0),46,0)<br>copy_VarCoords_1(Dec_reorder,prob_dec)<br>prob_dec@_FillValue = -9.9621e+36<br><br>;==============================================<br><br>; Now plotting the Significance<br>  wks  = gsn_open_wks(&quot;pdf&quot;,&quot;Significance_test&quot;)  <br>  gsn_define_colormap(wks,&quot;testcmap&quot;)              ; choose colormap<br> <br>  res                      = True                ; make plot mods<br>  res@gsnDraw = True<br>  res@gsnFrame = True<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 = &quot;ManualLevels&quot; ; manually set cnlev<br>  res@cnMinLevelValF       = -1            ; min level<br>  res@cnMaxLevelValF       =  1                 ; max level<br>  res@cnLevelSpacingF      = 0.1                 ; contour level spacing<br><br> <br><br>  res@tiMainString         = &quot;Dec SST with JJAS rainfall Regular 1967-2012&quot;<br>  plot= gsn_csm_contour_map_ce(wks,prob_dec,res)<br>   <br> ;===========================================<br><br></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></div>