BermanPhase.h

#ifndef BERMANPHASE_H
#define BERMANPHASE_H

#include "Phase.h"
struct BermanData {
  char *label;  /* label for solid phase or solid phase components   */ 
  char *formula;/* character string formula                          */
  double mw;       /* molecular weight pure phase or component          */
  double hr;       /* ref thermodynamic data, constants at TR, PR       */
  double sr;
  double k[4];         /* k terms (0-3)*/
  double v[5];         /* v terms (0-3)*/
  double Tt;           /* transition temperature */
  double Ht;           /* transition enthalpy */
  double l[2];         /* l terms */
};

class BermanPhase: public Phase {
 public:
  BermanPhase();
  BermanPhase(BermanData* dat);
  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();

  friend class BermanLambdaPhase;
 protected:
  BermanData *dat;
  virtual void init(BermanData* bd);
 private:
  const double pbr;
  double pb;
  double G(double tk,double pb);
  double H(double tk,double pb);
  double S(double tk,double pb);
  double Cp(double tk,double pb);
  double dCpdT(double tk,double pb);
  double V(double tk,double pb);
  double dVdT(double tk,double pb);
  double dVdP(double tk,double pb);
  double d2VdT2(double tk,double pb);
  double d2VdTdP(double tk,double pb);
  double d2VdP2(double tk,double pb);
};

#endif