function phaseo(t)
*
*  is it liquid metal, or is it xtal lattice?
*  ao(6) contains the melting coefficients for oxygen.
*
      implicit double precision(a-h,o-z)

      common/ao/ao(6)

      phaseo=ao(6)
      do 1 i=1,5
         phaseo = phaseo * t + ao(6-i)
1     continue

      return
      end 

************************************************************************