rcm2points函数使用

lideqin

setfileoption("nc", "FileStructure", "Advanced")
fn = "SAND5min.nc"

fi = addfile(fn, "r")

setfileoption("nc", "Format",  "NetCDF4Classic")
printVarSummary(fi)

depth = fi->depth
lat2d = fi->lat
lon2d = fi->lon

；sand3d(depth,lat,lon)
sand3d = fi->SAND

lat0 = (/34.9/)
lon0 = (/-97.9/)

sandpoint = rcm2points(lat2d,lon2d,sand3d(0,:,:),lat0,lon0,0)


如上想把格点数据插值到格点上去，但是始终运行不对。提示错误如下：
fatal:Number of dimensions in parameter (0) of (rcm2points) is (1), (2) dimensions were expected
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 29 in file readsand.ncl

请各位看看是哪里的错？？

longlivehj

rcm2points不合适
试试linint2_points

lideqin

sandpoint = linint2_points(sand3d&lon,sand3d&lat,sand3d,False,lon0,lat0,0)

使用linint2_points函数以后提示：
warning:linint2_points: xi and yi must be monotonically increasing

longlivehj

lideqin 发表于 2014-12-22 13:52
sandpoint = linint2_points(sand3d&lon,sand3d&lat,sand3d,False,lon0,lat0,0)

使用linint2_points函 ...

经纬度单调递增，你看看纬度是不是从高到低排列了，把顺序反一下，比如sand3d&lat(::-1)。

lideqin

版主你好，在你的指导下，我仔细查看了原来的经纬度的数据sand3d&lon和sand3d&lat
但是发现sand3d&lon本来就是递增的顺序，sand3d&lat为递减的顺序
调整以后变为：
  sandpoint = linint2_points(sand3d&lon,sand3d&lat(::-1),sand3d,False,lon0,lat0,0)

可以完成插值，但是插值出来的结果貌似有问题的，我想问下，只需要设置经纬度就可以了吗？sand3d数据本身不需要调整吗？

lideqin

longlivehj 发表于 2014-12-22 14:56
经纬度单调递增，你看看纬度是不是从高到低排列了，把顺序反一下，比如sand3d&lat(::-1)。

版主，虽然我直接写程序(fortran)读出来了，但是还想知道NCL这个插值的问题。

Liujiaxin

你好，请问解决了吗？