grads 编写ctl处理nc文件

言墨

grads 编写ctl处理nc文件
grads小白，求问各位大佬，我针对nc文件编写了ctl文件，画图时却出现这个报错（已经搜索了好久相关问题，气象家园也有人问，但是好像都没说怎么解决）：
gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds dimension bound
Data Request Error:  Error for variable 'sfh'
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = sfh


这个是原nc文件：
      lev = 12
      lat = 180
   variables:
      float sfh ( lev, lat )
         average_op_ncl :       dim_avg_n over dimension(s): model
         _FillValue :   1e+20

      double lev ( lev )
         standard_name :        air_pressure
         long_name :    pressure
         positive :     down
         axis : Z
         units :        hPa

      float lat ( lat )
         standard_name :        latitude
         long_name :    latitude
         units :        degrees_north
         axis : Y


这是我编写的ctl：
dset sfh.nc
title correlation
dtype netcdf
undef 9999
xdef 1 linear 0 1
ydef 180 linear -89.5 1.0
zdef 12 levels 100 150 200 250 300 400 500 600 700 850 925 1000
tdef 1 linear JUN1960 1mon
vars 1
sfh 12 99  z,y ssffhh
endvars

言墨

有大佬知道吗  求帮助

伽蓝鸟

nc文件不是自描述的咩，为啥要给它编写ctl啊。dat文件才需要ctl啊

言墨

伽蓝鸟 发表于 2019-12-13 09:36
nc文件不是自描述的咩，为啥要给它编写ctl啊。dat文件才需要ctl啊

sdfopen打不开 所以我根据网上有些大佬的教程写了ctl   然后open 这个ctl  文件能打开了  但是一画图就报这个错误

香草拿铁

xdef 1 linear 0 1     x就一个？
ydef 180 linear -89.5 1.0   这个是从-89.5度到1.0度？
zdef 12 levels 100 150 200 250 300 400 500 600 700 850 925 1000
tdef 1 linear JUN1960 1mon
vars 1
sfh 12 99  z,y ssffhh  感觉这个也有点奇怪

感觉你是画一个剖面或者时间序列的东西。
建议楼主再看看写sdfopen的nc文件的ctl该怎么写

或者试试matlab

言墨

香草拿铁 发表于 2019-12-13 09:49
xdef 1 linear 0 1     x就一个？
ydef 180 linear -89.5 1.0   这个是从-89.5度到1.0度？
zdef 12 level ...

是的 我是用别的软件处理好的数据变量sfh是两维的（lev，lat）   

言墨

解决了 我别的软件处理后输出成grd就可以了  不过还是不知道为什么nc文件会报这个错

苏拉苏拉

ctl一般zdef是从低层到高层吧。你的ctl从100到1000，不知道可不可以

大萌钰

我想问怎么改ctl时间啊

天宇2020

GRADS2.0之前的版本打不开nc文件，2.0，2.2和2.2可以打开