GizmoGas.H

Go to the documentation of this file.

 
#include <iostream>
#include <vector>
#include <string>
#include "../Definitions.H"
#include "../Core/FieldQty.H"
#include "../MPI/MPIUtil.H"
#include "../Utils/Constants.H"
#include "../Utils/Units.H"
#include "../Core/CRTree.H"
#include "../Losses/Losses.H"
#include "../Prob/Prob.H"
#include "../Propagation/Propagation.H"
#include "../Utils/Types.H"
#include "../Gas/Gas.H"
#include "GasTree.H"
#include "hdf5.h"
 
#ifndef _GIZMOGAS_H_
#define _GIZMOGAS_H_
 
namespace criptic {
  class GizmoReader {
 
  public:
  GizmoReader(std::string filename);
  
  template <typename T>
  static void giveUnits(std::vector<T>& unscaled, Real unit) {
    IdxType len = unscaled.size();
    for (IdxType i = 0; i<len; i++) {
        unscaled[i] *= unit;
    }
  }
  
  // Internal Data
 
  std::string file; 
  std::vector<double> smooth; 
  std::vector<double> density; 
  std::vector<double> ionDensity; 
  // Element abundances except metals are all fractions compared to all hydrogen
  std::vector<double> hIfrac; 
  std::vector<double> hIIfrac; 
  std::vector<double> heIfrac; 
  std::vector<double> heIIfrac; 
  std::vector<double> heIIIfrac; 
  std::vector<double> elAbundance; 
  std::vector<double> metals; 
  std::vector<RealVec> velocity; 
  std::vector<RealVec> position; 
  std::vector<RealVec> magField; 
  };
    
  namespace gas{
    class GizmoGas : public Gas {
      public:
      // CONSTRUCTORS & DESTRUCTORS
 
      GizmoGas(std::string filename, Real hLimit = 0, int leafSize_ = 32);  
      
      GizmoGas(const ParmParser& pp,
          const Geometry& geom_,
          const IdxType nBB,
          const int nGhost);
        
      ~GizmoGas() {
          delete data;
          delete particleTree;
      }
      // Important functions for loading data
 
      virtual GasData gasData(const RealVec &x, const Real t) const override;
      
      static GizmoGas* loadData(const std::string filename, 
         Real hLimit=0, int leafSize = 32) {
        return new GizmoGas(filename, hLimit, leafSize);
      };
      static GizmoGas* loadData(const std::string filename,
         const ParmParser& pp, const Geometry& geom_,
         const IdxType nBB_, const int nGhost_ = 1) {
        Real hMax = 0;
        int lSize = 32;
        int verbosity = 1;
        pp.get("verbosity", verbosity);
        pp.get("gas.tree.maxH", hMax);
        pp.get("gas.tree.leafSize", lSize);
        if (verbosity > 1 && MPIUtil::IOProc){
          std::cout << "GizmoGas: reading " << filename.c_str() << std::endl;
        }
        return GizmoGas::loadData(filename, hMax, lSize);
      };
      
      static Real loadTime(std::string filename);
      
      // Internal data
      private:
      std::string filename; 
      Real hLim; 
      int leafSize; 
      int verbosity; 
      GizmoReader* data; 
      GasTree* particleTree; 
    };
  }
}
 
//_GIZMOGAS_H_
#endif