- 积分
- 202
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-7-17
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
想用gs文件一次性画出1958年5-9月,逐候500hPa环流场,这是我的gs文件,请大家看看有什么问题,可以输出但是一个文件有N个图叠加,不知道原因,请大家指点!
'reinit'
iyear=1958
while(iyear<=1958)
'sdfopen D:\lw2011\huanliu\hgt.1958.nc'
'set lat 0 90'
'set lon 40 180'
*'set t 1'
'set lev 500'
'set xlopts 1 1 0.17'
'set ylopts 1 1 0.17'
'set grid off'
'set grads off'
'set map 1 1 10'
'set mpdset hires'
'set mpdset cnworld'
mon.5=may
mon.6=jun
mon.7=jul
mon.8=aug
mon.9=sep
*-----------------------------------------------
a=5
while(a<=9)
say mon.a' '
b=1
while(b<=5)
t1=(b-1)*5+1
t2=(b-1)*5+5
'define zz=ave(hgt,time='t1''mon.a''iyear',time='t2''mon.a''iyear')'
'd zz/10.0'
pull dummy
'enable print D:\lw2011\huanliu\'t1''mon.a'-'t2''mon.a''iyear'.500.gmf'
'print'
'disable print'
b=b+1
endwhile
t1=26
t2=31
if(a=6);t2=30;endif
if(a=9);t2=30;endif
'define zz=ave(hgt,time='t1''mon.a''iyear',time='t2''mon.a''iyear')'
'set cint 2'
'set gxout contour'
'set cthick 6'
'set cterp on'
'set csmooth on'
'd zz/10.0'
pull dummy
'enable print D:\lw2011\huanliu\'t1''mon.a'-'t2''mon.a''iyear'.500.gmf'
'print'
'disable print'
a=a+1
endwhile
'close 1'
iyear=iyear+1
endwhile
'reinit'
|
|