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

边缘人糊涂乐

nunu18 发表于 2014-3-21 19:15

                               
登录/注册后可看大图

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

不行啊，都做过了...

river

把close 1 去掉看看

边缘人糊涂乐

river 发表于 2014-3-21 19:59

                               
登录/注册后可看大图

把close 1 去掉看看

river 你终于出现啦，我都试过了，加不加‘close 1’结果都是画不错等值线图...求拯救。

river

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

                               
登录/注册后可看大图

river 你终于出现啦，我都试过了，加不加‘close 1’结果都是画不错等值线图...求拯救。

'set x 1'
'set y 1'
'set z 1'
'set t 1 10'
'define u=uwnd.2'
把这几句都去掉，后面的'd tcorr(u,qv.3,t=1,t=10)' 改成'd tcorr(uwnd.2,qv.3,t=1,t=10)'  试试

边缘人糊涂乐

river 发表于 2014-3-21 20:17

                               
登录/注册后可看大图

'set x 1'
'set y 1'
'set z 1'

试过了，这个改了等值线图就出错了，显示的错误是：DISPLAY error：Invalid expression  expression=tcorr(uwnd.2,qv.3,t=1,t=10) 之前我就试过这样改了，可是不对...

river

边缘人糊涂乐 发表于 2014-3-21 20:27

                               
登录/注册后可看大图

试过了，这个改了等值线图就出错了，显示的错误是：DISPLAY error：Invalid expression  expression=tcor ...

那就整体改造，使用set dfile来区分每个文件里的变量
······
去掉'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 contuor'
'd tcorr(u,qv.3,t=1,t=10)'
'print'
'disable print'
'reinit'
;

边缘人糊涂乐

river 发表于 2014-3-21 20:46

                               
登录/注册后可看大图

那就整体改造，使用set dfile来区分每个文件里的变量
······
去掉'close 1'

按照你的方法试过了，第二个等值线图变成了填色图，第一个填色图没有了~求拯救...

river

边缘人糊涂乐 发表于 2014-3-21 20:55

                               
登录/注册后可看大图

按照你的方法试过了，第二个等值线图变成了填色图，第一个填色图没有了~求拯救...

不好意思，明天再看看，今天要睡觉了

边缘人糊涂乐

river 发表于 2014-3-21 22:18

                               
登录/注册后可看大图

不好意思，明天再看看，今天要睡觉了

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

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'
;