Home > atmlab > h2o > thermodynamics > e_eq_ice.m

e_eq_ice

PURPOSE ^

E_EQ_ICE Equilibrium water vapor pressure over ice

SYNOPSIS ^

function ei = e_eq_ice(T)

DESCRIPTION ^

 E_EQ_ICE  Equilibrium water vapor pressure over ice

 Calculate the equilibrium water vapor pressure over a plane surface
 of ice according to Sonntag, 1994.
 
 The coefficients for units of K for temperature and Pa for the partial
 pressure are quoted in Helten et al., 1999.  I have verified that
 these coefficients are consistent with the ones given in Leiterer et
 al., 1997, with the only difference that the ones in the Leiterer et
 al. paper are for partial pressure in units of hPa.

 The temperature of 0C corresponds to 273.15K. (Not 273.16K, as stated
 in the Leiterer paper.)

 The formula (T in K, Ei in Pa):
 
 Ei(T) = exp[ a/T + b + cT + dT^2 + e ln(T) ]

 References:
 
 Sonntag, D., Advancements in the field of hygrometry, Meteorologische
 Zeitschrift, 3, 51-66, 1994.
 
 Helten, M. et al, In-flight comparison of MOZAIC and POLINAT water
 vapor measurements, JGR, 104, 26.087-26.096, 1999.
 
 Leiterer U. et al, Improvements in Radiosonde Humidity Profiles Using
 RS80/RS90 Radiosondes of Vaisala, Beitr. Phys. Atmosph., 70(4),
 319-336, 1997. 

 FORMAT   ei = e_eq_ice(T)

 OUT       ei = Equilibrium water vapor pressure in [Pa].
 IN       T  = Temperature in [K].

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

e_eq_ice.m

SOURCE CODE ^

0001 % E_EQ_ICE  Equilibrium water vapor pressure over ice
0002 %
0003 % Calculate the equilibrium water vapor pressure over a plane surface
0004 % of ice according to Sonntag, 1994.
0005 %
0006 % The coefficients for units of K for temperature and Pa for the partial
0007 % pressure are quoted in Helten et al., 1999.  I have verified that
0008 % these coefficients are consistent with the ones given in Leiterer et
0009 % al., 1997, with the only difference that the ones in the Leiterer et
0010 % al. paper are for partial pressure in units of hPa.
0011 %
0012 % The temperature of 0C corresponds to 273.15K. (Not 273.16K, as stated
0013 % in the Leiterer paper.)
0014 %
0015 % The formula (T in K, Ei in Pa):
0016 %
0017 % Ei(T) = exp[ a/T + b + cT + dT^2 + e ln(T) ]
0018 %
0019 % References:
0020 %
0021 % Sonntag, D., Advancements in the field of hygrometry, Meteorologische
0022 % Zeitschrift, 3, 51-66, 1994.
0023 %
0024 % Helten, M. et al, In-flight comparison of MOZAIC and POLINAT water
0025 % vapor measurements, JGR, 104, 26.087-26.096, 1999.
0026 %
0027 % Leiterer U. et al, Improvements in Radiosonde Humidity Profiles Using
0028 % RS80/RS90 Radiosondes of Vaisala, Beitr. Phys. Atmosph., 70(4),
0029 % 319-336, 1997.
0030 %
0031 % FORMAT   ei = e_eq_ice(T)
0032 %
0033 % OUT       ei = Equilibrium water vapor pressure in [Pa].
0034 % IN       T  = Temperature in [K].
0035 
0036 % 2003-03-28 Created by Stefan Buehler
0037 
0038 function ei = e_eq_ice(T)
0039 
0040 % Coefficients for Ei:
0041 a = -6024.5282;
0042 b = 29.32707;
0043 c = 1.0613868e-2;
0044 d = -1.3198825e-5;
0045 e = -0.49382577;
0046 
0047 ei = exp( a./T + b + c*T + d*T.^2 + e*log(T) );

Generated on Mon 15-Sep-2014 13:31:28 by m2html © 2005