Re: fatal:natgrids: z cannot contain any missing values

Hi Leo,

z is the name that function "natgrids" gives to your values to be gridded.
y_ano(:,iy) is your z variable here. It cannot contain missing values. Take a
look at

http://www.ncl.ucar.edu/Document/Functions/Built-in/natgrids.shtml

Regards,

-Ozan

Quoting LEO_ARIES <Leo_Aries_at_tom.com>:

> Hi, all
>
> I met a error when I read some scattered time series and want to
> derive a grid data before plotting. But I met a
> error , with a message as noted below:
>
> fatal:natgrids: z cannot contain any missing values
> fatal:Execute: Error occurred at or near line 71 in file
> Anomaly.contour_Gridding.ncl
>
> I haven't set any variable with a name of "z", and I don;t know
> where and why my code cause error.
>
> I also attached my script as below :
>
> ...........
> proxyindex =
> (/"21.2","30","28","27","24","22","52","53","01","02","77","13","70","14","11","51","36","40","45","44","46","530.25","530.28","530.18","530.19"/)
> lon =
> (/87.,88.,81.5,82.5,87.2,93.,81.,85.1,90.,98.,92.5,96.2,100.6,98.6,100.1,96.,103.8,106.1,110.,110.5,112.5,103.9,105.8,106.3,109.5/)
> lat =
> (/46.,48.,44.6,43.5,43.3,43.,35.,28.2,30.,30.,34.6,33.8,36.6,37.,38.9,38.,37.5,39.1,34.3,39.5,40.5,36.1,34.5,38.5,38.3/)
>
> ni = dimsizes(proxyindex)
>
> all_files = systemfunc ("ls " + "*10a.txt.nc")
> nf = dimsizes(all_files)
>
> y = new((/ni,50/),"float")
> y!0 = "proxyindex"
> y!1 = "time"
> y&time = fspan(1500,1990,50)
> y = -999.
> y@_FillValue = -999.
>
> do i = 0,ni-1 ; read all proxy data
> fi = addfile(all_files(i),"r")
> ;printVarSummary(fi)
> x = fi->px
> t1 = min(x&time)
> t2 = max(x&time)
> if(t1.lt.1500)then
> t1 = 1500
> end if
> y(i,{t1:t2}) = x({t1:t2})
> delete(x)
> end do
>
>
> mean = dim_avg(y(:,{1680:1970})) ; calculate the std
> anomalies for 1680-1970
> std = dim_stddev(y(:,{1680:1970}))
> print(mean)
> y_ano = new((/ni,50/),"float")
> do ii = 0,ni-1
> y_ano(ii,:) = (y(ii,:)-mean(ii))/std(ii)
> end do
> printVarSummary(y_ano)
>
> pgrid = new((/50,81,51/),"float") ; initialize the
> proxy gridded field
> do iy = 0,49
> ; =================================
> ; define new grid
> ; =================================
> nlon = fspan(75,115,81)
> nlon_at_long_name = "longitude"
> nlon_at_units = "degree_east"
> nlat = fspan(25,50,51)
> nlat_at_long_name = "latitude"
> nlat_at_units = "degree_north"
> time = fspan(1500,1990,50)
> time_at_long_name = "time"
> time_at_units = "year"
> ; =================================
> ; convert Proxy. --> lon-lat grid
> ; =================================
> pgrid(iy,:,:) = natgrids(lon,lat,y_ano(:,iy),nlon,nlat)
> end do
>
> ; =================================
> ; output to NetCDF file
> ; =================================
> pgrid!0 = "time"
> pgrid!1 = "lon"
> pgrid!2 = "lat"
> pgrid&time = time
> pgrid&lon = nlon
> pgrid&lat = nlat
> pgrid_at_long_name = "Decadal Standard Anomaly of PRCP"
> pgrid_at_units = " "
> system("rm -f p160.nc")
> f = addfile("Decadal_Std.Anomaly_proxy.PRCP_WestChina.nc","c")
> f->Pd = pgrid(time|:,lat|:,lon|:)
>
>
> Thanks !
>
>
>
> Lin
>
>
>
>
> LEO_ARIES
> 2007-04-11
>

----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Apr 10 2007 - 14:06:49 MDT