function h2o_vapr_ice_rapp05,t

;------------------------------------------------------------------
; Purpose:
;   Compute the saturation vapor pressure of h2o over ice according
;   to Rapp and Thomas, JASTP, 2005.
;
;   The R&T relationship is valid for T > 110 K  
;
; Input:
;     t.......temp in K
;
; Output:
;     p_sat...saturation vapor pressure (mbar)
;
; Source:  Mark Hervig
;------------------------------------------------------------------


   nt    = n_elements(t)
   p_sat = fltarr(nt)

   for i = 0, nt-1 do begin

     p_sat(i) = .01*exp( 9.550426 - 5723.265/t(i) + 3.53068*alog(t(i)) - 0.00728332*t(i) )

   endfor

   return,p_sat

   end