- 积分
- 75
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-8-11
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
最近看到一个求水汽通量的gs文件如下:
* This is a script for displaying moisture convergence
* Written by Michael Maxwell
*
* 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
* u = U-wind in m/s
* v = V-wind in m/s
* mconv = moisture convergence/divergence. convergence is positive and divergence is negative.
'reinit'
'open d:\12090.ctl'
'set lev 700'
'set lat 20 35'
'set lon 100 120'
'set t 5'
'set mpdset cnworld'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set grads off'
'set timelab off'
'set grid off'
#'tc=(tmpprs-273.16)'
'td=tc-((14.55+0.114*tc)*(1-0.01*rh) + pow((2.5+0.007*tc)*(1-0.01*rh),3) + (15.9+0.37*tc)*pow((1-0.01*rh),14))'
'vapr=6.112*exp((17.67*td)/(td+243.5))'
'e=vapr*1.001+(lev-100)/900*0.0034'
'define mixr=0.62137*(e/(lev-e))*1000/9.8'
#通过此gs文件算出的水汽通量的单位是g/(cm.s),一般为 十几到几十
'define qx=u*mixr'
'define qy=v*mixr'
'define mconv=hdivg(qx,qy)*1e6'
'enable print d:\2lb-shuiqi900.gmf'
'set arrscl 1 500'
#'set lon 124'
#'set cmin 30'
'set gxout shaded'
'd mag(qx,qy)'
#'set cmin 30'
'set gxout contour'
'd mag(qx,qy)'
#'d theta'
#'d skip(u,8);skip(v,8)'
'print'
'disable print'
谁能帮忙解释一下里面对露点温度和水汽压的公式计算由何处导出?是什么意思?
|
|