FeMetalMELTS.h

/* 
 * FeMetalMELTS class header file
 * Data from Tomiska
 * Mark Ghiorso
 */
#ifndef FEMETALMELTS_H
#define FEMETALMELTS_H

#include "Phase.h"

class FeMetalMELTS: public Phase {
 public:
  FeMetalMELTS();
  virtual ~FeMetalMELTS();
  virtual void setPa(double lpa);
  virtual double getGibbs();
  virtual double getEnthalpy();
  virtual double getEntropy();
  virtual double getCp();
  virtual double getdCpdT();
  virtual double getVolume();
  virtual double getdVdT();
  virtual double getdVdP();
  virtual double getd2VdT2();
  virtual double getd2VdTdP();
  virtual double getd2VdP2();
 private:
  virtual void init();
  const double pbr;
  double pb;
  double k0, k1, k2, k3, ct, v1, v2, v3, v4, h0, s0, v0;
};

#endif