<DIV>Dear all,</DIV>
<DIV><BR>I want to&nbsp;interpolate my station timeseries data to grids,which is with missing values.I&nbsp;find&nbsp;althought "natgrids" "<A href="http://www.ncarg.ucar.edu//ngmath/natgrid/nnhome.html">Natgrid</A>" and "natgridd" can deal with timeseries data,but they don't&nbsp;allow&nbsp;missing values in the input data arrays."obj_anal_ic"&nbsp;allow&nbsp;missing values in the input data arrays,but it design for one dimension array,and must do loop.Is there any other function&nbsp;</DIV>
<DIV>I&nbsp;didn't&nbsp;find which&nbsp;is&nbsp;according with my goal?</DIV>
<DIV>&nbsp;</DIV>
<DIV>Before find other function,I&nbsp;write a sript with "obj_anal_ic"&nbsp;(the attached),but it pop up&nbsp;following error messeges when run in ncl 5.1.1,and I don't know how to modify it.Can somebody help me please?</DIV>
<DIV>&nbsp;</DIV>
<DIV>Thanks in advance.</DIV>
<DIV>sunflower</DIV>
<DIV>&nbsp;</DIV>
<DIV>error messege as following(line 36 is&nbsp;where the function "<FONT color="#f00000">obj_anal_ic</FONT>" lie,in red color&nbsp;):</DIV>
<DIV>-----------------------------------------------------------------------------------</DIV>
<DIV>fatal:wgt_runave: The length of wgt must be less than or equal to the last dimen<BR>sion of x<BR>fatal:Execute: Error occurred at or near line 10354 in file $NCARG_ROOT/lib/ncar<BR>g/nclscripts/csm/contributed.ncl</DIV>
<DIV>fatal:Execute: Error occurred at or near line 10673 in file $NCARG_ROOT/lib/ncar<BR>g/nclscripts/csm/contributed.ncl</DIV>
<DIV>fatal:Execute: Error occurred at or near line 36 in file obj_anal_ic.ncl </DIV>
<DIV>---------------------------------------------------------------------------------------</DIV>
<DIV>&nbsp;</DIV>
<DIV>---------------------------------------------------------------------------------------</DIV>
<DIV>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<BR>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<BR>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</DIV>
<DIV>&nbsp;</DIV>
<DIV>begin<BR>&nbsp; npts = 5 <BR>&nbsp; ntim = 480&nbsp;&nbsp; ;extract 1961-2000 timeseries <BR>&nbsp; xx = asciiread ( "Yangtze_river_1961_2000_tas_observed_cof.txt", (/npts, ntim/), "float")<BR>&nbsp; xx!0&nbsp; = "pts"<BR>&nbsp; xx!1&nbsp; = "time"<BR>&nbsp; xObs = xx(time|:,pts|:)<BR>&nbsp; <A href="mailto:xObs@_FillValue">xObs@_FillValue</A> = -99.99<BR>&nbsp; printVarSummary(xObs )&nbsp;<BR></DIV>
<DIV>&nbsp; stain= asciiread("Yangtze_river_temperature_station.txt",(/5,3/), "float")&nbsp;&nbsp; <BR>&nbsp; printVarSummary(stain ) <BR>&nbsp; lat1 = stain(:,1)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </DIV>
<DIV>&nbsp; lon1 = stain(:,2)&nbsp;&nbsp;&nbsp; </DIV>
<DIV>&nbsp; <BR>&nbsp; lon = fspan(89, 98, 10)&nbsp;&nbsp; ;I output is 1.0*1.0&nbsp; <BR>&nbsp; nlon = dimsizes( lon)<BR>&nbsp; print(lon)<BR>&nbsp; lat = fspan(32, 36, 5)<BR>&nbsp; nlat = dimsizes( lat)<BR>&nbsp; print(lat)<BR>;============================================================== ;</DIV>
<DIV>&nbsp; grid =new((/nlat,nlon,ntim/),typeof(xObs)) <BR>&nbsp; do nmo =0,ntim-1<BR>&nbsp;&nbsp;&nbsp; zVal = xObs(nmo,:)<BR>&nbsp;&nbsp;&nbsp; rscan = (/5,3,2,1,0.5/)&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp; opt&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp; <A href="mailto:opt@zonal">opt@zonal</A> = True&nbsp;&nbsp; ; use zonal averages as a 1st guess<BR>&nbsp;&nbsp;&nbsp;<FONT color="#f00000"> grid(:,:,nmo) = obj_anal_ic(lon1,lat1,zVal,lon,lat, rscan, opt)</FONT>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;;line 36</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;printVarSummary(grid)&nbsp;&nbsp; <BR>&nbsp; end do</DIV>
<DIV>end</DIV>
<DIV>-------------------------------------------------------------------------------------</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;</DIV><br><br><span title="neteasefooter"/><hr/>
<a href="http://www.yeah.net/?from=footer">没有广告的终身免费邮箱,www.yeah.net</a>
</span>