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

大石头

感谢分享

wuzhx777

好东西哈，学习了！

大石头

请教楼主，我在运行的过程中循环不停止，下面的报错：
Unknown command: _vpg.78
SET Error: parea values beyond page limits
Data Request Warning:  Request is completely outside file limits
Cannot contour grid - all undefined values
Data Request Warning:  Request is completely outside file limits
Cannot contour grid - all undefined values
Drawing the sticks of the mproj latlon
incr lon = 30 incr lat = 30 length = 0.1
sticks will draw inside the frame
----------------------------------------------------------------
改过的脚本如下：
'reinit'

rc = gsfallow("on")
nrows=4
ncols=2
args=nrows%' '%ncols
'sdfopen E:\ncep\jan_w.nc'

*取消经纬比例自动调整
'set mproj scaled'
panels(args)
*p 表示当前所在的虚页是第几个虚页
p = 1
*ptot表示总共的虚页数量
ptot = nrows * ncols
*按panel循环并画图
while (p <= plot)
*'sdfopen E:\necp\jan_w.nc'

*设置第p个虚页
    _vpg.p
   'set parea 1 10 1 8'
   'set t 'p
   'set lon 102.5 111.5'
   'set lat 31 35'
   'set grads off'
   'set gxout shaded'
   'set font 2'
   'set  xlopts 1 4 0.2';'set  ylopts 1 4 0.2'
   'set csmooth on'
   'd w'
   'set gxout contour'
   'd w'
   '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 w(t='p',panel='p')'
   p = p + 1
endwhile
  'set vpage off'
'printim e:\ncep\test.png png white x800 y600'

本人grads刚入门，望楼主指点啊！！谢谢！！

kzpyezi

好东西要好好学习

拜伦大大

学习了

fashion

lnsqxt

好东西，虚页的麻烦死了

fzmmarco

最近刚好在学习用得到

niyan

好棒啊，正是我需要用的，大神啊

djl21

学习一下。