<div dir="ltr"><span class="" style="border-collapse:collapse;font-size:13px">Dear NCL,<div><br></div><div>I am able to get a reasonable plot with correlation with out detrending the data</div><div>.When I detrend it then the whole plot is one color.</div><div><br></div><div>Here is the printVarSummary of Snow Water Equivalent (SWE_82_12)</div><div><br></div><div><div>Variable: swe_82_12</div><div>Type: float</div><div>Total Size: 773523408 bytes</div><div>            193380852 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:<span class="" style="white-space:pre">        </span>[time | 372] x [lat | 721] x [lon | 721]</div><div>Coordinates: </div><div>            time: [ 672..271656]</div><div>            lat: [  30.. 750]</div><div>            lon: [   0.. 720]</div></div><div><br></div><div><div>snow water equivalent: min=-2   max=354.662</div></div><div><br></div><div>when the data is detrended the min max is </div><div><span class="" style="white-space:pre">        </span>min=-9.4341e+08   max=1.03525e+09<br></div><div><br></div><div>;Listed below is the whole code</div><div><div>a=addfile(&quot;/Users/ipshita/Desktop/GlobsnowV2.0/version02/<a href="http://monthly_82_13.nc">monthly_82_13.nc</a>&quot;,&quot;r&quot;)</div><div>swe=a-&gt;swe</div><div><br></div><div>time=a-&gt;time</div><div>lat=a-&gt;lat</div><div>lon=a-&gt;lon</div><div><br></div><div>print(dimsizes(lat))</div><div>print(dimsizes(lon))</div><div><br></div><div>swe_FillValue=-9.99e+08</div><div><br></div><div>utc_date = cd_calendar(time, 0)</div><div>year=utc_date(:,0)</div><div>month=utc_date(:,1)</div><div><br></div><div>;********************************************</div><div>;This is to extract 1982-2012 of swe data</div><div>;********************************************</div><div>swe_82_12=swe(0:371,:,:)</div><div>printVarSummary(swe_82_12)</div><div><br></div><div>;********************************************</div><div><br></div><div>;********************************************</div><div>;This is to read in All India Rainfall</div><div>;For June,July,August and Sept</div><div>;********************************************</div><div><br></div><div>JJAS_1871_2012 =asciiread(&quot;~/Documents/NCL_files/SST/Monthly/jjas_1871_2012.txt&quot;,(/142,1/),&quot;float&quot;)</div><div><br></div><div>;*******************************************</div><div>;Now this is to extract 1982 to 2012</div><div>;*******************************************</div><div><br></div><div>JJAS_82_12=JJAS_1871_2012(111:141,0)</div><div><br></div><div>;*******************************************</div><div>;============================================</div><div>;Extracting monthly data for sce</div><div><br></div><div>Dec_swe = new((/31,721,721/),float)</div><div><br></div><div>Dec_swe@_FillValue = -9.99e+08</div><div><br></div><div>do nyr=0,371,12</div><div><br></div><div>Dec_swe(nyr/12,:,:) =swe_82_12(nyr+11,:,:)</div><div><br></div><div>end do</div><div><br></div><div>Dec_swe@_FillValue = -9.99e+08</div><div><br></div><div>print(dimsizes(Dec_swe))</div><div><br></div><div>;=========================================</div><div>Dec_reorder=Dec_swe(lat|:,lon|:,time|:)</div><div><br></div><div>printMinMax(Dec_reorder,False)</div><div>Dec_reorder@_FillValue =-9.99e+08</div><div><br></div><div>Dec_dt=dtrend(Dec_reorder,False)</div><div><br></div><div>printMinMax(Dec_dt,False)</div><div>corr_Dec=esccr(Dec_dt,JJAS_82_12,1)</div><div><br></div><div><br></div><div>;************************************************</div><div>; create plot</div><div>;************************************************</div><div>  </div><div>wks = gsn_open_wks(&quot;x11&quot;,&quot;SWE_Dec_JJAS&quot;)               ; </div><div>gsn_define_colormap(wks,&quot;BlueRed&quot;)</div><div> </div><div> res = True              </div><div>  res@tiMainString         = &quot;SWE&quot; ; plot title</div><div> res@gsnAddCyclic        = False          ; Don&#39;t add a cyclic point.    </div><div>  res@cnFillOn             = True               ; turns on the color</div><div>  res@cnLinesOn            = False              ; turn off contour lines</div><div> res@cnLevelSelectionMode = &quot;ManualLevels&quot;     ; set manual contour levels</div><div>   </div><div>  plot = gsn_csm_contour(wks,corr_Dec(:,:,1),res)</div><div>;=================================================</div><div><br></div><div>printMinMax(corr_Dec,False)</div><div><br></div></div></span></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 30, 2016 at 2:41 PM, 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">Dear NCL,<div><br></div><div>I am able to get a reasonable plot with correlation with out detrending the data</div><div>.When I detrend it then the whole plot is one color.</div><div><br></div><div>Here is the printVarSummary of Snow Water Equivalent (SWE_82_12)</div><div><br clear="all"><div><br></div><div><div dir="ltr"><br></div></div>
</div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <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>