- 积分
- 8083
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-12-27
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2019-8-16 22:11:04
|
显示全部楼层
- ;假相当位温
- pre_levels = nlev(i)
- tmp_ave = tmp
- rhum_ave = 10^((7.5*6)/(237.3+6)) / 10^((7.5*10)/(237.3+10)) *100
- es = where( tmp_ave.lt.0, 6.112*exp(17.67* tmp_ave /( tmp_ave +243.5)),\
- 6.11*10^(7.5* tmp_ave /( tmp_ave +237.3)) )
- q = rhum_ave * (0.62197*es/( pre_levels - 0.378*es))/100.
- e = pre_levels * q/(0.62197+q)+10^(-10)
- tlcl = 55.0+2840.0/(3.5*log(tmp_ave +273.15)-log(e)-4.805)
- theta = ( tmp_ave +273.15)*(1000. / pre_levels )^(0.2854*(1.0-0.28*q))
- eqt = theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q)) ;假相当位温
复制代码 |
|