- 积分
- 6146
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-9-27
- 最后登录
- 1970-1-1
|
发表于 2013-8-19 22:41:54
|
显示全部楼层
今天在学习楼主的脚本,有点不明白,想请教一下。下面是在grads里面看见的他人的脚本:
- #饱和水汽压的计算,单位:hPa
- 'define es=6.112*exp(17.67*(t-273.15)/(t-29.65))'
- #饱和比湿的计算,单位:克/千克
- 'define qs=1000*0.622*es/(prs-0.378*es)'
- #求实际比湿
- 'define q=rh*qs/100'
- #计算水汽通量 |V|*q/g
- 'define q1=q*u/g'
- 'define q2=q*v/g'
下面是楼主的,想请问楼主的公式是从哪里来的?求指导
tt 代表温度
rhprs 代表相对湿度
tc = tt-273.16
lev = tt@lv_ISBL3
td=tc-((14.55+0.114*tc)*(1-0.01*rhprs)+(((2.5+0.007*tc)*(1-0.01*rhprs))^3)+(15.9+0.37*tc)*((1-0.01*rhprs)^14))请问这一行是什么意思?
vapr=6.112*exp((17.67*td)/(td+243.5)) ;饱和水汽压(为什么是td+243.5,请问楼主的公式是从哪来的?)
ee=vapr*1.001+(lev-100)/900*0.0034 ;这个是计算?
mixr=0.62137*(ee/(lev-ee))*1000/9.8
qu=u*mixr
qv=v*mixr
|
|