grads中nc资料转换grd出现的画图问题

可乐加盐

我在用1982年的一年nc资料，提取其中一个时次的500hPa数据，画出一个图
然后我用下面这个gs文件 'reinit'

'sdfopen D:\dailyncep\hgt.1982.nc'
'set fwrite D:\dailyncep\500hgt19820901.grd'
'set gxout fwrite'
'set t 244'
'set lev 500'
'set lon 0 360'
'set lat -90 90'
'd hgt'
'disable fwrite'
'reinit'
然后grd的ctl如下
dset D:\dailyncep\500hgt19820901.grd
title hgt
undef -9.99e+33
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 linear 500 1
tdef 1 linear 00Z001SEP1982 1440mn
vars 1
hgt 0 99 hgt
endvars
然后用ctl画出500hPa的图，但是不对。

可乐加盐

图反正就是有区别 不知道怎么挂上来。

可乐加盐

已解决了
是ctl中的格点读取有问题

Jungle

可乐加盐 发表于 2012-9-3 20:34

                               
登录/注册后可看大图

图反正就是有区别 不知道怎么挂上来。

不应该是 set lon 0 360  应该是set lon 0 357.5 或者set x 1 144
这里有个链接，你看看，希望对你有帮助http://bbs.lasg.ac.cn/bbs/thread-58575-1-1.html

long10001

学习过程中

njzqxt

楼主，能看看你正确的CTL文件吗？能不能说说你写CTL的思路呢，undef是根据什么得出来的呢？谢谢。

可乐加盐

Jungle 发表于 2012-9-3 20:41

                               
登录/注册后可看大图

不应该是 set lon 0 360  应该是set lon 0 357.5 或者set x 1 144
这里有个链接，你看看，希望对你有帮助 ...

对，是这个格点问题，在ctl中把格点改了也行。

可乐加盐

njzqxt 发表于 2012-9-4 08:45

                               
登录/注册后可看大图

楼主，能看看你正确的CTL文件吗？能不能说说你写CTL的思路呢，undef是根据什么得出来的呢？谢谢。

正确的ctl是把x的144改成145. 在读取的时候少读取了一个格点。

夜深沉

想问一下，为什么大虾们做物理量图时的gs文件，没有将某一天的nc变化为grd，而是直接打开呢？菜鸟，咨询~~

夜深沉

我按着做的，怎么不出图