grads处理grib2数据

hhhhh

我用g2ctl处理grib2数据。我在grads里输入的命令是：               ! g2ctl (-v) e:/cfsr/pgbh06.gads.icec.sfc.grb2 >e:/cfsr/icec.ctl
   把(-v)加上括号是表示我加上-v还是去掉结果是一样的，我出来的ctl为：
dset E:/cfsr/pgbh06.gdas.ICEC.SFC.grb2
index E:/cfsr/pgbh06.gdas.ICEC.SFC.grb2.idx
undef 9.999E+20
title E:/cfsr/pgbh06.gdas.ICEC.SFC.grb2
*  produced by g2ctl v0.0.3d
* griddef=1:0:(720 x 361):grid_template=0: lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 0.500000 lon 0.000000 to 359.500000 by 0.500000 #points=259920

dtype grib2
ydef 361 linear -90.000000 0.5
xdef 720 linear 0.000000 0.500000
tdef 1 linear 00Z00dec0000 1mo
zdef 1 linear 1 1
vars 1
ICECsfc  0,1,0   10,2,0,193 ** surface none Ice Cover [Proportion]
ENDVARS


在这个ctl中tdef只有1个时次，而我的数据是从1979年到2011年的月平均资料，请高手帮忙，应该怎样设置才会有正确的时间维。谢谢各位大神。

hhhhh

补充一下，数据都在一个文件中。就是从cfsr下载的pgbh06.gdas.ICEC.SFC.grb2

传说中的谁

使用论坛的搜索功能，查找关键字“template”

river

把时间改成你的所有资料的长度应该就可以了

hhhhh

传说中的谁 发表于 2013-4-3 21:50

                               
登录/注册后可看大图

使用论坛的搜索功能，查找关键字“template”

我查了一下，template应该是用一个ctl描述多个文件的控制参数，我的所有数据都在一个文件里，g2ctl生成的ctl中应该有时间维的变化，我的ctl没有，时间维为1，我也不知道为什么。这个ctl应该是错了，我打开这个ctl画图时出现
all undefined values
数据肯定不都是缺省值，
谢谢 传说中的谁

hhhhh

river 发表于 2013-4-3 22:33

                               
登录/注册后可看大图

把时间改成你的所有资料的长度应该就可以了

改了，不可以。画图时出现错误信息
  GRIB2 I/O error： message not found
data request error： error for variable 'icecsfc'
     eorror ocurred at column 1
DISPLAY error: Invalid expression
   Expression = icecsfc
只要set t 大于1时，就会出现这个错误信息
谢谢 river

river

hhhhh 发表于 2013-4-4 08:43

                               
登录/注册后可看大图

改了，不可以。画图时出现错误信息
  GRIB2 I/O error： message not found
data request error： erro ...

既然能够生成ctl还能画一个时次的图，如果图示正确的那就说明ctl应该不会有错。你在仔细看看你的原始资料的说明，到底里面有多少个时次的

hhhhh

river 发表于 2013-4-4 08:56

                               
登录/注册后可看大图

既然能够生成ctl还能画一个时次的图，如果图示正确的那就说明ctl应该不会有错。你在仔细看看你的原始资料 ...

你误解我的意思了，生成的ctl不能画图，只是在
set t 大于1时会出现上述的错误，set t 1时 会有all undefined values

mofangbao

你用meteoinfo打开试试

river

hhhhh 发表于 2013-4-6 10:01

                               
登录/注册后可看大图

你误解我的意思了，生成的ctl不能画图，只是在
set t 大于1时会出现上述的错误，set t 1时 会有all und ...

没太用过这个资料，实在不行这个资料要是不大的话，你传上来，大家都去试试说不定就有办法了