function goff_gratch_ice,t ;------------------------------------------------------- ; ; purpose: ; The Goff-Gratch formulation is used to compute ; saturation water vapor pressure over ice ; (Smithsonian Met. Tables, 5th ed., pp. 350, 1984) ; ; input: ; t.......temp in K ; ; output: ; result..saturation vapor pressure (mb) ; ; source: Mark Hervig ; ;------------------------------------------------------------- term = -9.09718 * (273.16/T - 1.0) $ -3.56654 * alog10(273.16/ T) + 0.876793 $ *(1.0 - T/ 273.16) + alog10(6.1071) result = 10.0^(term) return,result end