UniformGas.H

Go to the documentation of this file.

 
#ifndef _UNIFORMGAS_H_
#define _UNIFORMGAS_H_
 
#include "Gas.H"
#include "../IO/ParmParser.H"
 
namespace criptic {
  namespace gas {
 
    class UniformGas : public Gas {
 
    public:
 
      UniformGas(const ParmParser &pp) {
 
    // Read data from ParmParser
    pp.get("gas.density", gd.den);
    pp.get("gas.ionDensity", gd.ionDen);
    pp.get("gas.magField", gd.B);
    pp.get("gas.velocity", gd.v);
    gd.dx = 0.0;
    pp.query("gas.dx", gd.dx); // No intrinsic length scale, but
                   // user can choose to set one --
                   // this is necessary in some test
                   // problems where the problem
                   // itself is completey scale-free
 
    // Set derivatives of all gas quantities to zero
    gd.denGrad = gd.ionDenGrad = zeroVec;
    gd.vGrad = gd.BGrad = zeroTensor;
 
    // See if user as specified a composition, and, if so, use
    // that to set abundances
    std::string comp;
    pp.query("gas.comp", comp);
    if (comp.size() == 0) { 
      // Not specified; default to atomic composition, but let
      // user override abundances manually
      gd.setAtomicComposition();
      pp.query("gas.xH0", gd.xH0);
      pp.query("gas.xHp", gd.xHp);
      pp.query("gas.xHe0", gd.xHe0);
      pp.query("gas.xHep", gd.xHep);
      pp.query("gas.xHep2", gd.xHep2);
      pp.query("gas.xe", gd.xe);
      pp.query("gas.Z", gd.Z);
    } else {
      if (comp == "atomic") gd.setAtomicComposition();
      else if (comp == "molecular") gd.setMolecularComposition();
      else if (comp == "ionized") gd.setIonizedComposition();
      else {
        std::string err = "unknown gas composition: " + comp;
        MPIUtil::haltRun(err, errBadInputFile);
      }
    }
 
    // Check for radiation fields
    pp.query("gas.TBB", gd.TBB);
    pp.query("gas.WBB", gd.WBB);
    if (gd.TBB.size() != gd.WBB.size())
      MPIUtil::haltRun("gas.TBB and gas.WBB must have same number "
               "of entries", errBadInputFile);
         
    // Set TNB basis data
    b = TNBVectors(gd.B, gd.BGrad);
      }
 
      virtual ~UniformGas() { };
 
      virtual void frame(const RealVec &x,
             const Real t,
             RealVec& v,
             TNBBasis& tnb) const override {
    v = gd.v;
    tnb = b;
      }
 
      virtual GasData gasData(const RealVec &x,
                  const Real t) const override {
    return gd;
      }
 
    private:
 
      TNBBasis b;      
      GasData gd;      
    };
 
  }
}
 
#endif
// _UNIFORMGAS_H_