function temsat,press
; **** version of 22.09.81
;     to get kelvin temperature at which 'press' is the
;     saturation vapor pressure(mb) over water
;     zero is returned if the given vapor pressure is out of range of
;     the approximation (.0636 to 123.3972)

      vp = alog10(press)
      temsat = -.225896152438e+2 + vp*(.261012286592e+2   $
       + vp*(.30206720594e+1 + vp*(.370219024579e+0       $
       + .72838702401e-1 * vp))) + 273.16
      if((press lt .0636) or (press gt 123.3972)) then temsat=0.0        
          
      return,temsat
      end