- 积分
- 364
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-10-10
- 最后登录
- 1970-1-1
|
发表于 2014-10-10 21:34:34
|
显示全部楼层
这是我最近修改的,不知道对不
* rh = Relative Humidity in %
* t = Temp at *set level in degrees Kelvin
* tc = Temp in degrees C
* td = Dewpoint at *set level in degrees C
* e = Vapor pressure
* mixr = Mixing ratio
* ugreprs = U-wind in m/s
* vgrdprs = V-wind in m/s
* mconv = moisture convergence/divergence. convergence is positive and divergence is negative.
'tc = (tmpprs-273.16)'
'td = tc-((14.55+0.114*tc)*(1-0.01*rhprs) + pow((2.5+0.007*tc)*(1-0.01*rhprs),3) + (15.9+0.117*tc)*pow((1-0.01*rhprs),14))'
'vapr = 6.112*exp((17.67*td)/(td+243.5))'
'e = vapr*1.001+(lev-100)/900*0.0034'
'mixr = 0.62197*(e/(lev-e))*1000'
'mconv = hdivg(ugrdprs*mixr/9.8,vgrdprs*mixr/9.8)*1e-2'
'd mconv'
'd mconv'
|
|