<div dir="ltr"><div><div><div><div>TYour comments do not make any sense to me. Actually, they are rather frustrating. <br><br></div>[1] The CDO  and NCL are not returning the correct results. They are two independent software tools. Something does not seem right here.<br><br></div>[2] I am sure that neither were designed to work with &#39;weekly data&#39;.  &#39;Weekly data&#39; can span across neighboring months.<br>NCL&#39;s function requires that &#39;time&#39; be recognized by &#39;cd_calendar&#39;. Are the units of your variable  &#39;weeks since ...&#39;? Certainly NCL and, likely, the CDO are not going to partition a week that spans two months. <br><br></div>[3] You make statements like &quot;<font size="2"><span style="font-weight:normal">detrend not working right and calculate_monthly_va</span></font><font size="2"><span style="font-weight:normal">lues&quot;. You give no particulars.<br></span></font></div><font size="2"><span style="font-weight:normal">Why do you say these are not working correctly?<br><br><br></span></font></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Feb 16, 2016 at 1:14 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 have been working on making a plot with correlation between detrended data (snow cover extent) and I had first</div><div>used CDO to convert weekly data to monthly data. In the process the lat and lon cordinates became y and x.</div><div>It works fine if till I reorder it. But it is not calculating dtrend correctly, all the values are fill values. </div><div><br></div><div>I worked on using calculate_monthly_values to see if I could use it convert weekly to monthly in NCL and I could not use it. I would be grateful for any input</div><div><br></div><div>Here is the code I am using</div><div><br></div><div><div>;============================================</div><div>; set up parameters</div><div>;============================================</div><div><br></div><div>;This is extracting latitude of Western Europe<br></div><div>;Reading in the data</div><div>a=addfile(&quot;~/Documents/NCL_files/Snow_cover_IMS/NCL_codes/<a href="http://sce_mon_mean.nc" target="_blank">sce_mon_mean.nc</a>&quot;,&quot;r&quot;) ; open output netCDF file</div><div>b=addfile(&quot;~/Documents/NCL_files/Snow_cover_IMS/<a href="http://nhsce_v01r01_19661004_20140602.nc" target="_blank">nhsce_v01r01_19661004_20140602.nc</a>&quot;,&quot;r&quot;)</div><div><br></div><div>;Extracting  snow cover extent (sce) from respective data files</div><div><br></div><div>sce=byte2flt(a-&gt;sce)</div><div><br></div><div>sce@_FillValue = 9.9621e+36</div><div><br></div><div>lat=b-&gt;latitude</div><div>lon=b-&gt;longitude</div><div><br></div><div>sce@lat2d=lat</div><div>sce@lon2d=lon</div><div><br></div><div>;==============================================</div><div>;Reading in JJAS rainfall </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>;Extracting JJAS for 1967 to 2012</div><div><br></div><div>JJAS_1967_2012=JJAS_1871_2012(96:141,0)</div><div><br></div><div>;============================================</div><div>;Extracting monthly data for sce</div><div><br></div><div>Jan_sce = new((/46,88,88/),float)</div><div><br></div><div>Jan_sce@_FillValue = 9.9621e+36</div><div><br></div><div>Jan_sce@lat2d=lat</div><div>Jan_sce@lon2d=lon</div><div><br></div><div>;==========================================</div><div><br></div><div>do nyr=0,551,12</div><div><br></div><div>Jan_sce(nyr/12,:,:) =sce(nyr,:,:)</div><div><br></div><div>end do</div><div><br></div><div>Jan_reorder=Jan_sce(y|:,x|:,time|:)</div><div><br></div><div>Jan_dt=dtrend(Jan_reorder,False)</div><div><br></div><div><br></div><div>Jan_reorder@_FillValue = 9.9621e+36</div><div><br></div><div>corr_jan=escorc(Jan_dt,JJAS_1967_2012)</div></div><div><div><br></div><div><div dir="ltr"><br></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>