【已解决】Grads同时读取三个ctl，做时间序列与场相关及其检验

边缘人糊涂乐

感谢river的细心指导，终于把“时间序列与场相关及其检验”的填色图和等值线图画出来了~

gs和ctl如下：

'reinit'
'open e:\tskt.CTL'
'open e:\uwnd.CTL'
'open e:\qv.CTL'
'set parea 2.50 8.50 3.00 8.00'
'enable print e:/pingjun.gmf'
'set grads off'
'set grid off'
'set z 1'
'set t 1'
'set lon 60 150'
'set lat 0 50'
'set gxout shaded'
'set rgb 16 175 175 175'
'set clevs   -2.101 2.101'
'set ccols  16 0 16'
'd ttest'
*'cbarn 0.8 1 8.65 5.7'
'close 1'
*填色图是检验图

'set dfile 2'
'set x 1'
'set y 1'
'set z 1'
'set t 1 10'
'define u=uwnd'

'set dfile 3'
'set lon 60 150'
'set lat 0 50'
'set z 1'
'set t 1'
'set gxout contour'
'd tcorr(u,qv.3,t=1,t=10)'
*等值线图是时间序列与场相关图
'print'
'disable print'
'reinit'
;


DSET E:\tskt.grd
UNDEF -9.96921E+36
TITLE daily means ttst from the NCEP Reanalysis
XDEF 144 LINEAR 0 2.5
YDEF 73  LINEAR -90 2.5
ZDEF 1   LINEAR 0 1
TDEF 1   LINEAR 00Z01JAN1983 1dy
VARS 1
ttest 1 99 t-test dataset
ENDVARS

DSET E:\qv.grd
undef -9.96921E+36
TITLE daly means q1 from the NCEP Reanalysis
XDEF 144 LINEAR 0 2.5
YDEF 73  LINEAR -90 2.5
ZDEF 1  LINEAR 0 1
TDEF 10 LINEAR 00Z01JAN1983 1dy
VARS 1
qv 1 99 qv dataset
ENDVARS


DSET E:\a.grd
undef -9.96921E+36
TITLE daly means uwnd from the NCEP Reanalysis
XDEF 1 LINEAR 0 1
YDEF 1  LINEAR 0 1
ZDEF 1  LINEAR 0 1
TDEF 10 LINEAR 00Z01JAN1983 1dy
VARS 1
uwnd 1 99 u dataset
ENDVARS

画出的图形如下：


下面的错的图和数据的附件删不了，就放着吧~

river

边缘人糊涂乐 发表于 2014-3-21 22:24

                               
登录/注册后可看大图

好的，river大神记得有时间看看哈，非常感谢！

做出来没？如果还没做出来我就要批评你几句了······你画图的时候运行 gs 只提示了数据超出范围？没有别的错误提示？
无论做什么都要仔细认真，一步一步把能解决的问题都解决掉，才能出来结果······没有耐心的话那就先锻炼一下自己的耐心好吧。
我把正确答案给你，你自己看问题在哪。找到的话，我希望你能回复一个信息给我。
'reinit'
'open D:\data\tskt.CTL'
'open D:\data\Uwnd.CTL'
'open D:\data\Q1.CTL'
'set parea 2.50 8.50 3.00 8.00'
'enable print D:\data\pingjun.gmf'
'set grads off'
'set grid off'
'set z 1'
'set t 1'
'set lon 60 150'
'set lat 0 50'
'set gxout shaded'
'set rgb 16 175 175 175'
'set clevs   -2.101 2.101'
'set ccols  16 0 16'
'd ttest'
'cbarn 0.8 1 8.65 5.7'
*单独画这个填色图可以出图

'set dfile 2'
'set x 1'
'set y 1'
'set z 1'
'set t 1 10'
'define u=uwnd'

'set dfile 3'
'set lon 60 150'
'set lat 0 50'
'set z 1'
'set t 1'
'set gxout contour'
'd tcorr(u,qv,t=1,t=10)'
*单独画这个等值线图也可以出图
'printim D:\data\pingjun.gif white'
'print'
'disable print'
'reinit'
;

nunu18

楼主你想要什么样的图，说清楚呀，什么叫同时出图？你的意思是等值线图叠加在阴影图上吗？

lqouc

把你的ctl的时间都对应上再试。

边缘人糊涂乐

nunu18 发表于 2014-3-21 17:49

                               
登录/注册后可看大图

楼主你想要什么样的图，说清楚呀，什么叫同时出图？你的意思是等值线图叠加在阴影图上吗？

是啊 同时画在一张图上

边缘人糊涂乐

lqouc 发表于 2014-3-21 17:56

                               
登录/注册后可看大图

把你的ctl的时间都对应上再试。

时间，有两个是10，有一个是1，怎么对应呢？

nunu18

撇开时间怎么对应，首先要这两张图出现在一张图上，你要设置的经纬度阴影图和等值线图必须一样，再阴影图，最后画等值线图，看你gs文件也没问题，应该可以画出来的

lqouc

边缘人糊涂乐 发表于 2014-3-21 18:45

                               
登录/注册后可看大图

时间，有两个是10，有一个是1，怎么对应呢？

这没关系，步长和起始一样就行？还是画不出来？

边缘人糊涂乐

nunu18 发表于 2014-3-21 18:58

                               
登录/注册后可看大图

撇开时间怎么对应，首先要这两张图出现在一张图上，你要设置的经纬度阴影图和等值线图必须一样，再阴影图， ...

先画阴影图，再画等值线图，经纬度一样，然后就是画不出来，不知道是怎么回事儿。。。

边缘人糊涂乐

lqouc 发表于 2014-3-21 19:00

                               
登录/注册后可看大图

这没关系，步长和起始一样就行？还是画不出来？

嗯，嗯，画不粗来...单独画都可以画出来，一起画就是画不粗来..

nunu18

边缘人糊涂乐 发表于 2014-3-21 19:02

                               
登录/注册后可看大图

先画阴影图，再画等值线图，经纬度一样，然后就是画不出来，不知道是怎么回事儿。。。

要不楼主尝试，先把阴影图的变量按你的要求存成二进制文件，这样时间部长只有一维，再打开此文件，画阴影图，再画等值线图