function tvpice,press 
;     ice is at saturation.  this is the inverse of :vpice: 
      a0=-.203188177e+2
      a1=.2394167436e+2 
      a2=.2252719878e+1
      a3=.1914055442e+0
      a4=.9636593860e-2  

      vp = press  
      v =alog10(vp)  
      tvpce= a0+v*(a1+v*(a2+v*(a3+v*a4)))
      tvpice=tvpce+273.16 
  
      if ((vp lt 1.403e-5) or (vp gt 6.108)) then tvpice=0.
      return,tvpice
      end