grads作图问题

凌寻

请问一下  原维数环境是lon 60 150 lat 0 40 lev 850 850 time 1982:1:1:0 1982:1:1:0
然后我用编了一个gs文件  用set 重新设置维数 lat 0 40 lon 0 20 z 1 t 1
然后生成图形后显示为 Entire Grid Undefined
我想知道为什么？
是不是维数设置有什么要求的？

wet510

不知道你的源数据是什么格式的，不过你原位数环境是lev 850，在gs中你把set z 1，比较下两者是否一致。你可以把set z1 试试改成set lev 850

凌寻

wet510 发表于 2012-3-16 22:45

                               
登录/注册后可看大图

不知道你的源数据是什么格式的，不过你原位数环境是lev 850，在gs中你把set z 1，比较下两者是否一致。你可 ...

为什么是3呢  我觉得我需要深刻挖掘课本

wet510

这个也不一定是3，不过nc数据的一般高度层是1000 925 850,700....，具体怎样还要看你的数据，如果你是用ctl打开，建议你核对下你的ctl文件描述的环境和你的数据是否对应

mofangbao

贴出你的转换前后的ctl和转换的gs

冰枫

可能是你在转的时候没有把数据写进去，或者没对应上

易小凯

原来的lon范围是60~150，你set lon 0 20？！

凌寻

mofangbao 发表于 2012-3-16 23:32

                               
登录/注册后可看大图

贴出你的转换前后的ctl和转换的gs

data.ctl文件
dset d:\grads\data.grd
undef -9.99E+33
title NCEP/NCAR REANALYSIS PROJECT
xdef   37 linear   60.000  2.500
ydef   17 linear   0.000  2.500
zdef    2 levels   850  200
tdef 48 linear   JAN1982      1mo
vars   4
U  2 99 u wind (m/s)
V  2 99 v wind (m/s)
H  1 99 H500
TSFC 1 99 TSFC DATA
endvars
下边是我打算改的维数编的gs
'renint'
'open d:\grads\DATA.ctl'
'enable print d:\grads\b.gmf'
'set lat 0 40'
'set lon 0 20'
'set z 3'
'set t 1'
'd u'
'print'
'disable print'
;

然后图像显示的就是Entire Grid Undefined

mofangbao

你的Z根本就没有第三层啊 哪来的 set z 3呢 还有你楼上那位说的问题，楼主先看看基础知识理解下grads的ctl文件吧

凌寻

mofangbao 发表于 2012-3-17 17:55

                               
登录/注册后可看大图

你的Z根本就没有第三层啊 哪来的 set z 3呢 还有你楼上那位说的问题，楼主先看看基础知识理解下grads的ctl文 ...

恩 好吧 那还不是我设置的在所给区域范围就可以呢~