Gill (1982)の巻末付録にのっている式で計算しました。
tcは温度で単位は摂氏
esは飽和蒸気圧で単位はhPa
として、Fortranで書くとこうなります。
! Saturation vapor pressure computed with the equation described in Gill (1982)
POW=(0.7859 + 0.03477*tc)/(1.0 + 0.00412*tc) !tc in degree Celsius
es=10.0**POW !es in hPa
より正確には、下記の補正項fwをかける必要があります。pは気圧で単位はhPaです。
POW=(0.7859 + 0.03477*tc)/(1.0 + 0.00412*tc) !tc in degree Celsius
es=10.0**POW !es in hPa
fw=1.0+1.E-6*p*(4.5+0.0006*tc**2) !p in hPa
es=es*fw
図はこうなります。Petty (2008)のFig. 7.3 (p. 172)と比較してチェックしました。
satr_vapor_press.eps - 2011/07/20 17:52、A M (バージョン 1) 削除
17KB 表示 ダウンロード
! References
Atmosphere-ocean dynamics. A. E. Gill. Academic Press (London). 1982, 663 pp.
Petty (2008): A First Course in Atmospheric Thermodynamics. この本は大変良いがAmazonからは入手不可（直接著者のサイトから注文する）．