怎样用GrADS用两个CTL文件各自要画的图画在一张上

tbag

在画第二张图时候前加上set dfile 2

tbag

open 1.ctl
open 2.ctl

set dfile 1
set parea 图1坐标
d variable.1

set dfiel 2
set parea 图2坐标
d variable.2

1126492257

tbag 发表于 2018-2-27 00:16
open 1.ctl
open 2.ctl

按照你的方法加了仍然不行，还是只能读取第一个CTL文件的数据，画出的两张图都是第一个CTL文件的

qxtjls

什么都要贡献，是提问呀？还是挣。。。。

river

1126492257 发表于 2018-2-26 20:14
加上之后提示这样的错误，也只能画出一张图

'reinit'
'open f:\bht\3\fnl_20180105_12_00.ctl'
'open f:\bht\3\fnl_20180131_12_00.ctl'
'enable print f:\bht\3\div.gmf'
'set parea 1.0 5.2 0.8 7.8'
'set grads off'
'set csmooth on'
'set lev 900 100'
'set t 1 18'
'define qdiv=hdivg(ugrdprs.1,vgrdprs.1)*100000'
'set ccolor 1'
'set cint 1'
'set clskip 2'
'set cthick 5'
'set lat 27.2'
'set lon 102.9'
'd qdiv'

'set dfile 2''reset'
'set parea 5.5 10 0.8 7.8'
'set lev 900 100'
'set t 1 15'
'define qdiv=hdivg(ugrdprs,vgrdprs)*100000'
'set ccolor 1'
'set cint 1'
'set clskip 2'
'set cthick 5'
'set lat 27.2'
'set lon 102.9'
'd qdiv'

'print'
'disable print'
;
试试这个

柚子皮

'reinit'
'open f:\bht\3\fnl_20180105_12_00.ctl'
'open f:\bht\3\fnl_20180131_12_00.ctl'
'enable print f:\bht\3\div.gmf'
'set parea 1.0 5.2 0.8 7.8'
'set grads off'
'set csmooth on'
'set lev 900 100'
'set t 1 18'
'define qdiv=hdivg(ugrdprs.1,vgrdprs.1)*100000'
'set ccolor 1'
'set cint 1'
'set clskip 2'
'set cthick 5'
'set lat 27.2'
'set lon 102.9'
'd qdiv'

'set parea 5.5 10 0.8 7.8'
'set lev 900 100'
'set t 1 15'
'define qdiv.2=hdivg(ugrdprs.2,vgrdprs.2)*100000'   读取的变量要区别开来，同时定义的变量也区别一下比较好
'set ccolor 1'
'set cint 1'
'set clskip 2'
'set cthick 5'
'set lat 27.2'
'set lon 102.9'
'd qdiv.2'

'print'
'disable print'

任性小王子

怎么把观测数据插值成格点，然后用GRADS画图，请教大神！