- 积分
- 1184
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-3-28
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
'reinit'
'open D:/grads/gs/202007.ctl'
i=20
while(i<=40)
'set t 'i''
'q time'
aa=subwrd(result,3)
hour=substr(aa,1,2)
day=substr(aa,4,2)
month=substr(aa,6,3)
year=substr(aa,9,12)
'set mpdset cnworld'
'set map 1 1 4'
'set grads off'
'set lev 1000 200'
'set parea 1.5 10.5 1 8'
'set zlog off'
'define div=hdivg(ugrdprs,vgrdprs)*100000'
'define vor=hcurl(ugrdprs,vgrdprs)'
'set lat 42.6'
'set lon 123.4 124'
'set gxout shaded'
'set annot 1 3'
'set csmooth on'
'set clab on'
*'d vvelprs'
'd vor'
'cbarn'
'set gxout contour'
'set cthick 5'
*'d vvelprs'
'd vor'
'draw title vor 'year''month''day''hour''
'printim D:/grads/gs/202007/vor'year''month''day''hour'.png white x1000 y1000'
'c'
i=i+1
endwhile
;
|
|