DMWvapor.h

/* DMWvapor header file
 * $Id: DMWvapor.h,v 1.7 2007/07/23 19:55:50 ghiorso Exp $
 * by Victor Kress
 */

#ifndef DMWVAPOR_H
#define DMWVAPOR_H


#include "Phase.h"
#include "SpecSolNID.h"

class DMWvapor: public SpecSolNID {

public:
  // public variables ////////////////////////////////
  const int NCOMP;
  const int NSPEC;
  const int H2DEX;
  const int N2DEX;
  const int O2DEX;
  const int S2DEX;
  const int CL2DEX;
  const int CODEX;
  const int H2SDEX;
  const int CH4DEX;
  const int CO2DEX;
  const int H2ODEX;
  const int SO2DEX;

  // creation/destruction ////////////////////////////////
  DMWvapor();
  virtual ~DMWvapor();

  // setting methods ////////////////////////////////////
  virtual void setTk(double ltk);
  virtual void setPa(double lpq);
  virtual double getGibbs();
  virtual double getEnthalpy();
  virtual double getEntropy();
  virtual void getVolume();

 protected:
  double pb
  double Z;
  double B,C,D,E;
  double dZdVm,dZdTm;
  double Tm;
  double Pm;
  double Vm;
  double sigma;
  double epsilon;
  virtual void init();
  double sigmavec[11];
  double epsilonvec[11];
  static double a[15];
  static double **k[2];
  static double allstoich[11][6];
  virtual double mufunc(int p,double *x);    //implements speciateNID form
  virtual double dmufunc(int m,int p,double *x);//implements speciateNID form
  virtual double muxs(int ispec,double *x);
  virtual double dmuxs(int m,int p,double *x);
  virtual void update() {
    if (!updated) {
      updateMvars();
      SpecSolution::update(); //updated variable reset here, so keep this last
    }
    return;
  }
  virtual void updateMvars()  throw(PhaseError *);
};

#endif