function goff_gratch_liq,t

;-------------------------------------------------------
;
; Purpose:
;   The Goff-Gratch formulation is used to compute
;   saturation vapor pressure over liquid water
;   (Smithsonian Met. Tables, 5th ed., pp. 350, 1984)
;
; Parameters:
;
;  input:
;     t.......temp in K
;
;  output:
;     result...saturation vapor pressure (mb)
;
;  source:  Mark Hervig
;
;-------------------------------------------------------------

   term = -7.90298 * (373.16/T - 1.0) $ 
          + 5.02808 * alog10(373.16/T) $
          - 1.3816e-7 * (10.0^(11.344 * (1. - T/373.16)) -1.0) $
          + 8.1328e-3 * (10.0^(-3.49149 * 373.16/T) -1.0) $
          + alog10(1013.246)

   result = 10.0^(term) 

   return,result
   end