求助，grads提取wrf数据

3120101008

请问我用grads提取wrfout的数据为什么不对？下面是我提取数据的gs文件，wrfout01.ctl是后处理出来的原始ctl文件。
'reinit'
'open e:\grads\wrfout1\wrfout01.ctl'
*'set x 70 140'
*'set y 15 55'
'set z 1'
'set t 31'
'set fwrite e:\grads\wrfout2.dat'
'set gxout fwrite'
'd RAINNC+RAINC'
'disable fwrite'
;
这是wrfout2.dat的ctl文件，和原始wrfout01.ctl的内容相比，除了变量名不一样，z和t都改成1了，其他都是一样的
dset e:\grads\wrfout2.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.4.1 MODEL
pdef  117  78 lcc  37.000  103.000   59.000   39.500  47.00000  25.00000  103.00000  60000.000  60000.000
xdef  389 linear   50.30641   0.27027029
ydef  190 linear    8.84747   0.27027029
zdef   1 linear 1 1  
tdef   1 linear 00Z01JUL2011    1440MN      
VARS  1
rain  1 0 ACCUMULATED TOTAL PRECIPITATION (mm)
endvars

river

看看是不是缺测值的问题，改成-9.99e+08试一下

yangcs

butaihui o

3120101008

river 发表于 2014-3-18 07:56

                               
登录/注册后可看大图

看看是不是缺测值的问题，改成-9.99e+08试一下

还是不对呢

river

3120101008 发表于 2014-3-19 16:01

                               
登录/注册后可看大图

还是不对呢

那就比较麻烦了，是不是模式转的有问题，所以资料有问题

黄小仙儿

我用wrf转出来画的图也和你的差不多，现在也在纠结是什么问题。http://bbs.06climate.com/forum.php?mod=viewthread&tid=21264

3120101008

river 发表于 2014-3-19 17:12

                               
登录/注册后可看大图

那就比较麻烦了，是不是模式转的有问题，所以资料有问题

我原始的dat和ctl是可以正常画出来图的，所以应该不是模式转的问题，是fwrite提取wrf数据生成dat文件的的时候，这个东西就出问题了

3120101008

黄小仙儿 发表于 2014-3-19 18:26

                               
登录/注册后可看大图

我用wrf转出来画的图也和你的差不多，现在也在纠结是什么问题。http://bbs.06climate.com/forum.php?mod=vi ...

额。。。你原始的生成的dat和ctl就画不出来？

黄小仙儿

3120101008 发表于 2014-3-19 19:15

                               
登录/注册后可看大图

额。。。你原始的生成的dat和ctl就画不出来？

画得出来，但是很怪异。。

3120101008

黄小仙儿 发表于 2014-3-19 19:41

                               
登录/注册后可看大图

画得出来，但是很怪异。。

我转出来的原始数据是对的，提取数据之后，出现了问题