<div dir="ltr"><div>Not sure of the issue. I have attached month_to_annual_weighted.ncl.<br></div><div>Make sure you are using this one.<br></div><div><br><br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>load "./month_to_annual_weighted.ncl"<br></div>begin<br><div> <br> diri = "/Users/nmb1441/research/papers/china/monsoons/cru_data/pet/"<br> srcFileName = "<a href="http://cru_ts3.21.1951.1960.pet.dat.nc">cru_ts3.21.1951.1960.pet.dat.nc</a>" <br> <br> refile= addfile(diri + srcFileName,"r")<br><br> pet = refile->pet<br> printVarSummary(pet)<br> <br> yearfrac = cd_calendar(pet&time, 4) ; same as time = f->timeNot sure of the issue<br> yyyymm = cd_calendar(pet&time, -1)<br><br> print(yyyymm+" "+yearfrac)<br> <br> pet_annual_total1= month_to_annual_weighted(yyyymm,pet,0)<br> printVarSummary(pet_annual_total1)<br> <br>end<br><br></div><div>Seems to work fine.<br><br></div><div><br>Variable: pet_annual_total1<br>Type: double<br>Total Size: 20736000 bytes<br> 2592000 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [year | 10] x [lat | 360] x [lon | 720]<br>Coordinates: <br> year: [1951..1960]<br> lat: [-89.75..89.75]<br> lon: [-179.75..179.75]<br>Number Of Attributes: 3<br> _FillValue : 9.969209968386869e+36<br> long_name : potential evapotranspiration<br> units : mm/day<br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 10, 2015 at 2:08 PM, Lei Meng <span dir="ltr"><<a href="mailto:dream916@gmail.com" target="_blank">dream916@gmail.com</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">Hi,<div><br></div><div>I am still unable to get the function month_to_annual_weighted.ncl to work. I used the function code that you sent. Here is my script:</div><div><br></div><div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div>load "./month_to_annual_weighted.ncl"</div><span class=""><div>begin</div><div> </div><div> diri = "/Users/nmb1441/research/papers/china/monsoons/cru_data/pet/"</div><div> srcFileName = "<a href="http://cru_ts3.21.1951.1960.pet.dat.nc" target="_blank">cru_ts3.21.1951.1960.pet.dat.nc</a>" </div><div> refile= addfile(diri + srcFileName,"r")</div><div><br></div><div> pet = refile->pet</div><div> lat = refile->lat</div><div> lon = refile->lon</div><div> time = refile->time</div><div> </div><div> time@calendar = "gregorian"</div><div> yearfrac = cd_calendar(time, 4)</div><div> </div></span><div> yyyymm = cd_calendar(time, -1)</div><div> </div><div> dym = cd_calendar(time, -1) </div><div> </div><div> print(dym)</div><div><br></div><div> printVarSummary(pet)</div><div> printVarSummary(dym)</div><div> pet_annual_total1= month_to_annual_weighted(yyyymm,pet,0)</div><div> </div><div>end</div></div><div><br></div><div>The error message I got is : </div><div>
<p><span>fatal:Variable (yrFrac) is undefined, can not assign attribute (calendar)</span></p>
<p><span>fatal:["Execute.c":8578]:Execute: Error occurred at or near line 52 in file ./month_to_annual_weighted.ncl</span></p>
<p><span></span><br></p>
<p><span>fatal:["Execute.c":8578]:Execute: Error occurred at or near line 28 in file test-temp1.ncl</span></p><p><span><br></span></p><p><span>I am using ncl 6.2.1.</span></p><p><span>Thanks!</span></p><span class="HOEnZb"><font color="#888888"><p><span>Lei</span></p><p><span><br></span></p></font></span></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jan 28, 2015 at 1:33 PM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.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>[1] Please indicate which version of NCL you are using.<br><br></div>[2] I speculate the following will work<br><div><br> yyyymm = cd_calendar(time, -1) ; <===<br><div> pet_annual_total1= month_to_annual_weighted(yyyymm,pet,0)<br><br></div><div>[3] Since, I am not sure what version of NCL you have, please load the attached<br><br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br></div><div>load "./month_to_annual_weighted.ncl"<br><br></div><div>... your script ...<br><br></div><div>[4] Note: from 6.2.0 onward, users do not need to explicitly load the gsn_code/gsn_csm/contributed libraries<br><br></div><div>Goos Luck<br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Wed, Jan 28, 2015 at 9:54 AM, Lei Meng <span dir="ltr"><<a href="mailto:dream916@gmail.com" target="_blank">dream916@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr">Hi,<div><br></div><div>I tried to convert monthly values to annual valued weighted by the number of days. Here is my code:</div><div><br></div><div>----</div><div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"</div><div>begin</div><div> </div><div> diri = "/Users/nmb1441/research/papers/china/monsoons/cru_data/pet/"</div><div> srcFileName = "<a href="http://cru_ts3.21.1951.1960.pet.dat.nc" target="_blank">cru_ts3.21.1951.1960.pet.dat.nc</a>" </div><div> refile= addfile(diri + srcFileName,"r")</div><div><br></div><div> pet = refile->pet</div><div> lat = refile->lat</div><div> lon = refile->lon</div><div> time = refile->time</div><div> </div><div> time@calendar = "gregorian"</div><div> yearfrac = cd_calendar(time, 4)</div><div> </div><div> yyyymm1 = cd_calendar(time, 1)</div><div> </div><div> dym = cd_calendar(time, 1) </div><div> </div><div> print(dym)</div><div><br></div><div> printVarSummary(pet)</div><div> printVarSummary(dym)</div><div> pet_annual_total1= month_to_annual_weighted(dym,pet,0)</div><div> </div><div>end</div></div><div>--</div><div>The summaries of variables pet and dym are as follows:</div><div><br></div><div>
<p><span>Variable: pet</span></p>
<p><span>Type: double</span></p>
<p><span>Total Size: 248832000 bytes</span></p>
<p><span> 31104000 values</span></p>
<p><span>Number of Dimensions: 3</span></p>
<p><span>Dimensions and sizes:<span> </span>[time | 120] x [lat | 360] x [lon | 720]</span></p>
<p><span>Coordinates: </span></p>
<p><span> time: [18642..22264]</span></p>
<p><span> lat: [-89.75..89.75]</span></p>
<p><span> lon: [-179.75..179.75]</span></p>
<p><span>Number Of Attributes: 5</span></p>
<p><span> long_name :<span> </span>potential evapotranspiration</span></p>
<p><span> units :<span> </span>mm/day</span></p>
<p><span> correlation_decay_distance :<span> </span>-999</span></p>
<p><span> _FillValue :<span> </span>9.969209968386869e+36</span></p>
<p><span> missing_value :<span> </span>9.969209968386869e+36</span></p><p>Variable: dym</p><p><span></span></p>
<p><span>Type: double</span></p>
<p><span>Total Size: 960 bytes</span></p>
<p><span> 120 values</span></p>
<p><span>Number of Dimensions: 1</span></p>
<p><span>Dimensions and sizes:<span> </span>[120]</span></p>
<p><span>Coordinates: </span></p>
<p><span>Number Of Attributes: 1</span></p>
<p><span> calendar :<span> </span>gregorian</span></p>
<p><span>(0)<span> </span>month_to_annual_weighted: mm must be 1-to-12 inclusive</span></p><p>The error message is " mm must be 1-to-12 inclusive". </p><p>dym values are from 195101 to 196012.</p><p>This is a very simple code. I couldn't figure out what was wrong with my code.</p><p>The data is too big and I will send it if necessary. Thanks for any suggestions. </p><span><font color="#888888"><p>Lei</p><p><br></p><p><span><br></span></p></font></span></div><div><br></div></div>
<br></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>