想把四个时次的图放在一张上 但是这个循环不对。。。。

caojunwei

                               
登录/注册后可看大图

'reinit'
rc = gsfallow("on")
nrows=2
ncols=2
args=nrows%' '%ncols
panels(args)
p=1
ptot=nrows*ncols
while (p<=ptot)
   
'sdfopen E:\ncep-data\data\2016\uwnd.2016.nc'
'sdfopen E:\ncep-data\data\2016\vwnd.2016.nc'
'sdfopen E:\ncep-data\data\2016\hgt.2016.nc'
'sdfopen E:\ncep-data\data\2016\air.2016.nc'
'E:/ncep-data/red_blue3.gs'
*''disk'/color/colornum.gs'
'enable print E:\pic\300hpa.gmf'
'set mproj scaled '
*------------
_vpg.p
'set parea 1 9.8 1 7.5'
year='2016'
mon='Apr'
lev=300
t=1
while (t<=2)
hr=0
while (hr<=23)
*'c'
'set time 'hr'z't''mon''year''
'set z 3'
'set lon 80 150'
'set lat 30 70'
'define div=hdivg(uwnd.1*100000,vwnd.2*100000)'
*'set mproj nps'
*-------draw-----
'set lon 80 150'
'set lat 30 70'
*'set vpage 0 11 0 8.5'
*'set parea 1 10.5 1 7.5'
'set grads off'
'set grid off'
'set mpdset cn lowres'
*'set mpdset cn cnriver lowres'
'set map 1 1 3'
'set xlint 10'
'set ylint 10'
'set xlopts 1 4 0.16'
'set ylopts 1 4 0.15'
'set gxout shaded'
'set clevs -2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2'
'd div'
cbarn.gs
'set gxout contour'
'set cint 4'
'set ccolor 2'
'set cstyle 2'
'set cthick 5'
'set clskip 2'
*'d air.4-273'
'set gxout contour'
'set cint 4'
'set ccolor 1'
'set cstyle 1'
'set cthick 5'
'set clskip 2'
'd hgt.3/10'
'draw title t ,h , 'lev'hPa,  'hr'z't''mon''year''
*****************************
hr=hr+12
'print'
endwhile
t=t+1
endwhile
p=p+1
endwhile
'set vpage off '
'disable print'
*-----------------------------

codd

{:eb513:}{:eb513:}