【一页多图】-利用GrADS自带脚本函数实现

风云翩跹

'reinit'
rc = gsfallow("on")
nrows=2
ncols=2
args=nrows%' '%ncols
*'sdfopen d:\pres.sfc.2013.nc'
*取消经纬比例自动调整
'set mproj scaled'
panels(args)
*p 表示当前所在的虚页是第几个虚页
p = 1
*ptot表示总共的虚页数量
ptot = nrows * ncols
*按panel循环并画图
while (p <= ptot)
*设置第p个虚页
   'open D:\ECMF\13110520.0'p*24'.ctl'
    _vpg.p
   'set parea 1 9.8 1 7.5'
   'set t 'p
   'set grads off'
*  'set gxout shaded'
   'set font 2'
   'set  xlopts 1 4 0.2';'set  ylopts 1 4 0.2'  
   'd m2g'
   'set font 0'
*  'cbar_interp 1 1 1'
*  'zb 30 30 0.1'
    'set font 5';'set strsiz 0.3 0.3'
   'draw string 3.5 8 pres(t='p*24',panel='p*24')'
   p = p + 1
   'close D:\ECMF\13110520.0'p*24'.ctl'
endwhile
  'set vpage off'
'printim D:\ECMF\pres2.png png white x1000 y800'
;
清风大大，我现在想把四个预报场画在同一张图上，它们的CTL文件是分开的，我现在用这个循环只能画出一张，您说的步骤一和步骤三调换并及时关闭我还是没太懂，您看看我这个文件哪里有问题

mofangbao

风云翩跹 发表于 2013-11-7 10:07

                               
登录/注册后可看大图

'reinit'
rc = gsfallow("on")
nrows=2

先检查你ctl文件名拼凑的是否有问题，你可以先用一个变量保存文件名，如果要动态生成，请先输出这个文件名变量，确保拼凑的文件名是正确的。
还有 如果你的文件时分开的，下面的set t要确保取到的是正确的时次

忘却

恩 很好 正好用得到！

jeffggg

清风版主威武！

倚天照海

楼主辛苦，感谢分享

wudi1201

不错，不错

yiluo521wzh

学习了

藍鶄

楼主，本人智力有限，能给个简单一点的例子么？谢了顺便问下函数是不是这个

鬼王

谢谢楼主分享

GWJjiayou

好棒~真在学习