【求助】怎么用Grads显示不同类型资料？求大神帮助！

Crazy_Merlin

各位大神好，我想利用Grads把晴空大气可降水量数据TPW和FY-2E IR1通道分区图同时显示出来，叠加到一张图~但是始终出不来~希望大家帮忙看看！
这是问我的gs文件：
'reinit'
'open F:\data\awx\ir1andtpw\tpw.ctl'
'set mpdset cnworld'
'set map 1 1 5'
'set lon 70 140'
'set lat 15 50'
'set grads off'
'set grid on'
'set xlint 5'
'set ylint 5'
'set gxout shaded'
*'set cint 10'
*'set cmin 0.0'
*'set cmax 75'
'd tpw'
*
pull dummy
*'c'
*'reinit'
'open F:\data\awx\ir1andtpw\ir1.ctl'
'set mpdset cnworld'
'set map 1 1 5'
'set lon 70 140'
'set lat 15 50'
'set grads off'
'set grid on'
'set xlint 5'
'set ylint 5'
'set gxout shaded'
'set display grey'
'd bt'
pull dummy
return
;


按照上面的gs文件操作的话，总会提示错误（如图）


附上数据和ctl文件：见附件

随缘

应该是d bt.2
因为你是要画第二个文件的变量

Crazy_Merlin

随缘 发表于 2012-12-7 10:16

                               
登录/注册后可看大图

应该是d bt.2
因为你是要画第二个文件的变量

你好，能具体点吗？我都是自学的，之前这个没用到~谢谢！

随缘

Crazy_Merlin 发表于 2012-12-7 10:19

                               
登录/注册后可看大图

你好，能具体点吗？我都是自学的，之前这个没用到~谢谢！

'd bt'
这边改成'd bt.2'

Crazy_Merlin

随缘 发表于 2012-12-7 10:28

                               
登录/注册后可看大图

'd bt'
这边改成'd bt.2'

额，还是没画出来~依然提示错误~我用的数据数组大小不同，TPW是1001*1001，IR1是1900*1300的~

随缘

Crazy_Merlin 发表于 2012-12-7 10:53

                               
登录/注册后可看大图

额，还是没画出来~依然提示错误~我用的数据数组大小不同，TPW是1001*1001，IR1是1900*1300的~

这个没影响的，你最好把数据的ctl贴出来

Crazy_Merlin

随缘 发表于 2012-12-7 10:57

                               
登录/注册后可看大图

这个没影响的，你最好把数据的ctl贴出来

数据和ctl都在附件里~

dset F:\data\awx\ir1andtpw\FY2E_ANI_IR1_R04_20120802_0030.AWX
title FY-2E_TPW data
options yrev
undef 99
*99为无效值
xdef 1900 linear  50.02 0.05
ydef 1300 linear -4.95 0.05
zdef 1 levels 850
tdef 11 linear 00:30Z02Aug2012 3hr
vars 1
bt  0 -1,40,1 FY-2E_TPW (mm)
endvars


dset F:\data\awx\ir1andtpw\FY2E_TPW_MLT_OTG_20120802_0000.AWX
title FY-2E_TPW data
options yrev
undef 99
*99为无效值
xdef 1001 linear  55 0.1
ydef 1001 linear -50 0.1
zdef 1 levels 850
tdef 11 linear 00Z02Aug2012 3hr
vars 1
tpw  0 -1,40,1 FY-2E_TPW (mm)
endvars

随缘

Crazy_Merlin 发表于 2012-12-7 11:54

                               
登录/注册后可看大图

数据和ctl都在附件里~

dset F:\data\awx\ir1andtpw\FY2E_ANI_IR1_R04_20120802_0030.AWX

你这个真正的问题是初始时间不一致造成的，将ctl里面时间设置改一致就可以了

淡水河

随缘 发表于 2012-12-7 12:24

                               
登录/注册后可看大图

你这个真正的问题是初始时间不一致造成的，将ctl里面时间设置改一致就可以了

两个都是shaded也会覆盖的吧。

river

淡水河 发表于 2012-12-7 12:52

                               
登录/注册后可看大图

两个都是shaded也会覆盖的吧。

两个都是shaded肯定会覆盖的，后画的会把先画的覆盖掉。而且你这个范围都不一样也不好比较啊，可以计算一个差值什么的再画图是不是更好一点