Grads处理nc资料时的问题

稚气未脱

首先我下载了ECMWF全球的再分析产品，然后用Grads编写文件如下：
'reinit'
'sdfopen F:\batch-netcdf-web239-20130619020130-HEewjg.cdf'
'set gxout fwrite'
'set fwrite F:\results\EC_2011080100.dat'
'set lon 0 360'
'set lat -90 90'
'set t 1'
'd t'
'disable fwrite'

问题一：发现与直接在Grads敲打命令得到的.dat资料大小不符，会不会是Grads软件的问题呢？
问题二：编写ctl文件如下：
dset F:\results\EC_2011080100.dat
title ECMWF xa field
undef -9.99e+33
xdef 360 linear 0 1.0
ydef 181 linear -90 1.0
zdef 1 linear 1000 1
tdef 1 linear 00z01AUG2011 360mn
vars 1
t=>t 17 t,z,y,x Temperature
endvars
发现画出来的图形有漂移，与直接用nc画出来的图对不上，找了很久都没发现问题所在，nc文件给出的头文件如下：

想请教各位大神是什么原因，多谢了!

river

图怎么漂移的啊，贴个图出来看看对比一下

稚气未脱

分别用二进制格式和nc格式读取的1000hPa的位势高度图  用二进制格式读的图似乎是在nc格式上向右平移了一些

做个霸气的木头

把'set lon 0 360'换成‘set x 1 360’试试

Aires

不要怀疑是grads的问题，先检查你自己的代码再说。
gs文件把层也设置上就可以了，grads必须都要set保证运算不出错。

jeffggg

问下楼主这个资料在哪里下载的？“我试试

稚气未脱

jeffggg 发表于 2013-6-23 09:37

                               
登录/注册后可看大图

问下楼主这个资料在哪里下载的？“我试试

http://apps.ecmwf.int/datasets/data/interim_full_daily/ 先注册才能下载

稚气未脱

做个霸气的木头 发表于 2013-6-22 23:10

                               
登录/注册后可看大图

把'set lon 0 360'换成‘set x 1 360’试试

嗯 两个问题都解决了 第一个 gs 文件结尾分号没有加  第二个 'set lon' 和'set lat'换成'set x'和'set y'就行了

dzxconan

应该set lon 0 259吧~？不然360和0是重复的 显示两次

qxtlyf

dzxconan 发表于 2013-6-28 17:07

                               
登录/注册后可看大图

应该set lon 0 259吧~？不然360和0是重复的 显示两次

应该set lon 0 259吧~？中259，应该是359