- 积分
- 60
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-4-12
- 最后登录
- 1970-1-1
|
GrADS
| 系统平台: |
|
| 问题截图: |
 |
| 问题概况: |
'reinit'
'open D:/FNL/fnl_2013.ctl'
tt=1
while(tt<=16)
'set t'tt
'set lev 975'
'set lat 0 40'
'set lon 60 150'
'enable print D:/FNL/picture/MPV/mpv975.gmf'
'define it=TMPprs'
'define rh=RHprs'
'define prs=lev'
'define es=(6.112*exp(17.67*(it-273.15)/(it-29.65)))'
'define q=rh*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(it)-log(e)-4.805)'
'define theta=it*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
'define vo=hcurl(UGRDprs,VGRDprs)'
'define f=2*7.292*sin(lat*3.14159/180.0)'
'define g=9.8'
'define dp=lev(1000)-lev(950)'
'define deqt=eqt(1000)-eqt(950)'
'define du=UGRDprs(1000)-UGRDprs(950)'
'define dv=VGRDprs(1000)-VGRDprs(950)'
'define dx=2.0*6370949.0*cos(lat*3.14159/180.0)*3.14159/180.0'
'define dy=2.0*6370949.0*3.14159/180.0'
'define dtx=cdiff(eqt,x)'
'define dty=cdiff(eqt,y)'
'define mpv1=-g*(vo+f)*deqt/dp'
'define mpv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define mpv=pv1+pv2'
'set gxout shaded'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set csmooth on'
'd mpv1'
'q dims'
it=sublin(result,5)
it=subwrd(it,6)
iz=sublin(result,4)
iz=subwrd(iz,6)
'draw title 'it' mpv for 'iz'hPa'
'cbar'
'print'
'c'
tt=tt+1
endwhile
'disable print'
; |
| 我看过提问的智慧: |
看过 |
| 自己思考时长(天): |
1 |
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看，没有帐号？立即注册
x
|
|
窥视卡
雷达卡
发表于 2014-5-2 15:46:29
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
发表于 2014-5-2 17:56:48
楼主