- 积分
- 787
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-8-9
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
以下是本人的gs文件,请高手帮忙!谢谢
'reinit'
'open e:\11\2014012700.ctl'
'set mpdset cnworld guangdong_q'
'set lon -30 130'
'set lat -10 80'
'set lev 700'
tt=1
while(tt<=30)
'set t 'tt''
'set grads off'
'set grid off'
'define p=lev'
'define e=exp(17.67*(TMPprs-273.15)/(TMPprs-29.65))*(RHprs/100)*6.112'
'define q=622*e/(p-0.378*e)'
'define q1=UGRDprs*q/9.8'
'define q2=VGRDprs*q/9.8'
'set lon -20 120'
'set lat 0 70'
'set lev 850'
'set gxout shaded'
'set cmin 2'
'd mag(q1,q2)'
*'d q1;q2'
'set gxout vector'
'set clab on'
'cbarb 1 0'
*'d UGRDprs;VGRDprs'
'd UGRDprs*11;VGRDprs*11'
'run d:\opengrads\contents\resources\scripts\cbarn.gs'
'q dims'
rec=sublin(result,5)
recs=subwrd(rec,6)
' draw title 700hPa vapourflux 'recs''
'enable print e:\tu\2014012700-700vf.gmf '
'print'
'c'
tt=tt+1
endwhile
'disable print'
|
|