- 积分
- 11182
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-10-29
- 最后登录
- 1970-1-1
|
发表于 2017-5-26 19:02:36
|
显示全部楼层
'define p=lev'
'define e=6.112*exp(17.67*(air-273.15)/(air-29.65))*(rhum.2/100)'
'define q=622*e/(p-0.378*e)'
*-- 水汽通量
'define q1=uwnd.3*q/9.8'
'define q2=vwnd.4*q/9.8'
*'define aq1=vint(p,q1,300.0)'
*'define aq2=vint(p,q2,300.0)'
*'d aq1;aq2'
*-- 水汽通量散度
'define qdiv=hdivg(q1,q2)'
*先画阴影
'set gxout shaded'
'set cmax 0'
'd qdiv'
*后画等值线
'set gxout contour'
'd qdiv'
*画0线加粗
'set clevs 0'
'set cthick 0'
'd qdiv'
'd q1;q2' |
|