fnl资料做高度位势图总是报错

丫丫yang

想ncl画850hpa的位势高度图，变量修改如下：
begin
  f=addfile("./fnl_20151224_12_00.grib2","r")
         u = f->UGRD_P0_L100_GLL0({850},:,:)
         v = f->VGRD_P0_L100_GLL0({850},:,:)
         t = f->TMP_P0_L100_GLL0({850},:,:)
        hgt= f->HGT_P0_L100_GLL0({850},:,:)
  
  newt = t - 273.15

报错信息为fatal:NclOneDValGetClosestIndex: finish coordinate index out of range, can't continue
fatal:Could not obtain coordinate indexes, unable to perform subscript
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 7 in file hgt_temp.ncl
求大神指点，应该怎么修改啊。。。

zongqing123

你读取的数据好像没有坐标信息，你给附一个坐标信息试一下

丫丫yang

zongqing123 发表于 2017-1-2 10:52
你读取的数据好像没有坐标信息，你给附一个坐标信息试一下

附了坐标还是不行，可能是我坐标的格式不对，不知道应该怎么给

jianglegejiang

{850}改为1到26中的数字

yeah...

丫丫yang 发表于 2017-1-2 12:28
附了坐标还是不行，可能是我坐标的格式不对，不知道应该怎么给

你先什么也不加 printVarSummary一下变量看看