AMRMelange.H

Go to the documentation of this file.

#ifdef CH_LANG_CC
/*
 *      _______              __
 *     / ___/ /  ___  __ _  / /  ___
 *    / /__/ _ \/ _ \/  V \/ _ \/ _ \
 *    \___/_//_/\___/_/_/_/_.__/\___/
 *    Please refer to Copyright.txt, in Chombo's root directory.
 */
#endif

#ifndef _AMRMELANGE_H_
#define _AMRMELANGE_H_

#include "AmrIce.H"
#include "LevelSigmaCS.H"
#include "SurfaceFlux.H"
#include "NamespaceHeader.H"

#define MELANGE_N_GHOST 1
#define MELANGE_N_COMP 1


class AMRMelange
{

  int m_time_step;
  Real m_time;
  Real m_diffusion_factor;
  
  //Vector<LevelData<FArrayBox>*> m_melange;
  int m_finestLevel;
  Vector<DisjointBoxLayout> m_grids;
  Vector<int> m_ratio;
  Vector<RealVect> m_dx;
  SurfaceFlux* m_external_source;

public:
  ~AMRMelange();
  AMRMelange();

  const LevelData<FArrayBox>* melange(int a_level) const;
  int finestLevel() const
  {
    return m_finestLevel;
  }
  void define(const Vector<DisjointBoxLayout>& a_grids, 
              const Vector<int>& a_ratio,
              int a_finestLevel,
              const RealVect& a_crseDx);

  void timestep(Real a_dt, AmrIce& a_amrIce);

  void increment(AmrIce& a_amrIce, Real a_scale);

  
#ifdef CH_USE_HDF5 
  //plot file methods

  void addPlotVars(Vector<std::string>& a_vars);

  void writePlotData(LevelData<FArrayBox>& a_data, int a_level);

  //check point methods
  
  void addCheckVars(Vector<std::string>& a_vars);
  
  void writeCheckData(HDF5Handle& a_handle, int a_level);
  
  void readCheckData(HDF5Handle& a_handle,  HDF5HeaderData& a_header, int a_level, const DisjointBoxLayout& a_grids);
#endif

  const LevelData<FArrayBox>* melangeThickness(int lev)
  {
    return m_melange[lev];
  }
  
protected:
  //cell-centered melange
  mutable Vector<LevelData<FArrayBox>* > m_melange;

private:
  
  void computeSource(Vector<LevelData<FArrayBox>* >& a_source,
                     AmrIce& a_amrIce,   Real a_dt);

  void updateMelange(Vector<LevelData<FArrayBox>* >& a_melange,
                     const Vector<LevelData<FArrayBox>* >& a_source,
                     const Vector<RefCountedPtr<LevelSigmaCS> >& geometry,
                     Real a_dt);

};
  
  

class MelangeIceObserver : public AmrIce::Observer 
{
  AMRMelange* m_melangePtr;
  bool  m_next_increment_positive; //used to make sure that m_melangePtr->increment calls occur in pairs
  
public:

  MelangeIceObserver();
  ~MelangeIceObserver();

  AMRMelange& melange() const;
  
  void notify(AmrIce::Observer::Notification a_n, AmrIce& a_amrIce);

  void addPlotVars(Vector<std::string>& a_vars);

  void writePlotData(LevelData<FArrayBox>& a_data, int a_level);

  //check point methods
  
  void addCheckVars(Vector<std::string>& a_vars);
  
  void writeCheckData(HDF5Handle& a_handle, int a_level);
  
  void readCheckData(HDF5Handle& a_handle,  HDF5HeaderData& a_header, int a_level, const DisjointBoxLayout& a_grids);

};


class MelangePhysIBC : public PhysIBC
{

public:

  MelangePhysIBC()
  {
  }

  void define(const ProblemDomain& a_domain,
              const Real&          a_dx);
  

  PhysIBC* new_physIBC();

  void initialize(LevelData<FArrayBox>& a_U);

  void primBC(FArrayBox&            a_WGdnv,
              const FArrayBox&      a_Wextrap,
              const FArrayBox&      a_W,
              const int&            a_dir,
              const Side::LoHiSide& a_side,
              const Real&           a_time);


  void setBdrySlopes(FArrayBox&       a_dW,
                     const FArrayBox& a_W,
                     const int&       a_dir,
                     const Real&      a_time);
  

  void artViscBC(FArrayBox&       a_F,
                 const FArrayBox& a_U,
                 const FArrayBox& a_divVel,
                 const int&       a_dir,
                 const Real&      a_time);   
};



#include "NamespaceFooter.H"

#endif