- 积分
 - 11183
 
	- 贡献
 -  
 
	- 精华
 
	- 在线时间
 -  小时
 
	- 注册时间
 - 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' |   
 
 
 
 |