使用shgetnp函数报错，看不懂报错

Nostalgia

请问，我运行了shgetnp这个函数，函数具体介绍如下链接：http://www.ncl.ucar.edu/Document/Functions/Built-in/shgetnp.shtml
然后这是我的部分代码

lat_1 = ispan(lat_min, lat_max+1, 1) ;26)
lon_1 = ispan(lon_min, lon_max+1, 1) ;28)
printVarSummary(lat_1)
printVarSummary(lon_1)

lat_11 = new(728, float)
lon_11 = new(728, float)
  z_11 = new(728, float)
  a    = new(nsta,integer)
lat_11@_FillValue = -32767
lon_11@_FillValue = -32767
  z_11@_FillValue = -32767
     a@_FillValue = -32767

do i = 0,lat_cha-1
    do j = 0,lon_cha-1
        lat_11(i*lon_cha+j) = lat_1(i)
        lon_11(i*lon_cha+j) = lon_1(j)
    end do
end do
z_11 = 1.
print(lat_11(720:727))
print(lon_11(720:727))
; print(z_11)

print(lat(0))
print(lon(0)) ;lat(0), lon(0), lon(0)
a(0) = shgetnp(39.8673, 116.366, 1.0,  tofloat(lat_11), tofloat(lon_11), tofloat(z_11), 0)
print(a(0))


他报错内容：

IERR=    7
SHSTORE3 - cell grid dimensions must be positive
IERR=    9
SHQSHEP - collinear input, no unique solution
fatal:shgetnp: Error number 9.
fatal:["Execute.c":8640]:Execute: Error occurred at or near line 283 in file try.ncl


麻烦大佬们解答一下

doice

请问您解决了吗，我碰到了和您一样的问题。

Nostalgia

doice 发表于 2021-1-31 09:51
请问您解决了吗，我碰到了和您一样的问题。

没有解决。。。最后我用的取整函数做的