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

tulalang

看看楼主的成果，学习学习~

泡沫学天气

看了好多清风大神的帖子 太牛啦

熊猫支队007

学习啦，好棒！grads虽然古老技能还是满分的！

福禄寿囍财

zym246810

写论文需要这样的图，过来学习

zym246810

想请教一下，在写到_vpg.p这个语句的时候，grads显示“Unkown command：”这是什么原因呢？
都是按照楼主的例子写的，只有调用分页函数中的全局变量：_vpg.p这里有问题，求大神帮助，这是我的gs文件

'reinit'
'open F:\lw\Module\pingfen\data\biatu\map.ctl'
'open F:\lw\Module\pingfen\ABCD\pingfen87shixiaohecheng.ctl'
'set grads off'
'set mproj scaled'
'set lon 114 123'
'set lat 29 35.5'
'set poli off'
nrows=3
ncols=4
f=1
args=nrows%' '%ncols
while(f<=1)
if(f=1);o=TS;endif;
'H:\IEdownload\panels.gsf 'args
k=1
while(k<=3)
j=1
while(j<=4)
p=j+(k-1)*4
_vpg.p
'set grads off'
'set t 'j+(k-1)*4
'define aa'1'=oacres(g.1,'o'.2)'
'define bb'1'=smth9(aa'1')'
'set csmooth on'
'set gxout shaded'
'set lon 114 123'
'set lat 29 35.5'
'set grid off'
'd bb'1
'set grid off'
'cbar'
'set gxout contour'
'set ccolor 1'
'set clab off'
'H:\tools\OpenGrADS\lib\province-basemap.gs jiangsu bb'1
'set line 4 1 6'
'draw shp H:\IEdownload\he\hyd1_4l.shp'
'set line 4 1 3'
'draw shp H:\IEdownload\he\hyd2_4l.shp'
'set line 0 3 4'
'draw shp H:\IEdownload\dijishi_2004.shp'
'set line 1 1 10'
'draw shp H:\IEdownload\bou2_4m\bou2_4l.shp'
if(k=1)
n=6h
else
if(k=2)
n=12h
else
n=24h
endif
endif
if(j=1)
m=Light
else
if(j=2)
m=moderate
else
if(j=3)
m=heavy
else
m=storm
endif
endif
endif
'draw title Precipitation accumulated in 'n%' hours, '%m%' rain, 018~072'
j=j+1
endwhile
k=k+1
endwhile
'set vpage off'
'printim F:\lw\Module\pingfen\ABCD\ABCDtu_PNGshixiaohecheng\'o'\'o'.png white'
'c'
f=f+1
endwhile

福尔摩斯特酷团

谢谢楼主分享！正好用到

langit

多谢

大气小白

看起来比较复杂啊，研究一下

zmy1548487100

好东西，感谢楼主分享