#include "frost_point.h"

#include <cmath>

float h2o_vapr_ice_rapp05( const float& t) 
{
/*
   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) )

	return 0.01*std::exp( 9.550426 - 5723.265/t + 3.53068*std::log(t) - 0.00728332*t ) ;
/*
   endfor

   return,p_sat

   end
*/

}