<div dir="ltr"><div><div><div>Hi,<br><br></div>I wrote a program to convert weekly data to daily. It works, but the output has only -9 throughout the matrix.<br>Can anyone tell me where I am going wrong<br><br>******************************************************************<br>;This is to create daily values from weekly values<br>;Created on Nov. 24th 2014<br><br><br>;Reading in the snow covered extent original file for time<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>;*******************************************<br><br>;Reading in the original data<br><br>a=addfile(&quot;~/Documents/Snow_cover_IMS/<a href="http://nhsce_v01r01_19661004_20140602.nc">nhsce_v01r01_19661004_20140602.nc</a>&quot;,&quot;r&quot;)<br><br>;Now extracting time from the file<br>time_sce=a-&gt;time<br><br>;Now we are going to extracting snow cover extent in area<br>b =addfile(&quot;~/Documents/Snow_cover_IMS/NCL_codes/sce_area_B.nc&quot;,&quot;r&quot;)<br>sce_area=b-&gt;B<br>sce_area_avg=wgt_areaave_Wrap(sce_area,1.0,1.0,0)<br>x=sce_area_avg(11:dimsizes(time_sce)-21)<br><br>;Now we want to convert the snow covered extent to daily values<br><br>;Since first year has only 12 data points we would like to remove that to have full years and last year 2014 also have 22 points<br><br>time=time_sce(11:dimsizes(time_sce)-21)<br><br>;Since we are reading data from 1967-2013, we have 46 years<br><br>tNew=new(46*365,&quot;float&quot;)<br><br>xo =linint1_n (time, x, False, tNew, 0, 0)<br>print(xo)<br><br>********************************************************************<br></div>Thank you<br>Best Regards<br></div>Ipshita<br></div>