;+
; NAME:
;       temp_along_satadiabat
;
; PURPOSE:
;       To calculate the temperature at a point along a saturated adiabat
; CATEGORY:
;       Function
; CALLING SEQUENCE:
;
; INPUTS:
;
; OPTIONAL INPUTS:
;	
; KEYWORD PARAMETERS:
;
; OUTPUTS:
;
; OPTIONAL OUTPUTS:
;
; COMMON BLOCKS:
;
; SIDE EFFECTS:
;
; RESTRICTIONS:
;
; PROCEDURE:
;
; EXAMPLE:
;
; MODIFICATION HISTORY:
;       Written by:     Wayne F. Feltz, CIMSS/SSEC 
;
; $Date:$       May 10, 1998
; $Id:$
; $Log:$
;
;
;-

function temp_along_satadiabat, thetae, p, p_surf

pcon=(p_surf/p)^.286
delta=20.0
t_k=273.13
eps=0.62197058

while (abs(delta) gt .1) do begin
   mixr=satmix(p,t_k,eps)
   th=t_k*pcon*exp(2.5*mixr/t_k)
   if (((th-thetae(0))*delta) gt 0.0) then delta=-0.5*delta
   t_k=t_k+delta

endwhile

return,t_k
end