MultiLevelDataIBC.H

Go to the documentation of this file.

#ifdef CH_LANG_CC
/*
*      _______              __
*     / ___/ /  ___  __ _  / /  ___
*    / /__/ _ \/ _ \/  V \/ _ \/ _ \
*    \___/_//_/\___/_/_/_/_.__/\___/
*    Please refer to Copyright.txt, in Chombo's root directory.
*/
#endif
//
// MultiLevelDataIBC.H
// ============
//
// PhysIBC-derived class which stores initial topography and thickness data
// and imposes either periodic or reflection boundary conditions
#ifndef _MULTILEVELDATAIBC_H_
#define _MULTILEVELDATAIBC_H_

#include "IceThicknessIBC.H"
#include "NamespaceHeader.H"


class MultiLevelDataIBC : public IceThicknessIBC
{
public:
  MultiLevelDataIBC(const Vector<RefCountedPtr<LevelData<FArrayBox> > >& a_thck,
                    const Vector<RefCountedPtr<LevelData<FArrayBox> > >& a_topg, 
                    const RealVect& a_dxCrse, const Vector<int> & a_refRatio);
                    
  virtual ~MultiLevelDataIBC();

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

  virtual IceThicknessIBC* new_thicknessIBC();

  virtual void initialize(LevelData<FArrayBox>& a_U);


  virtual void initializeIceGeometry(LevelSigmaCS& a_coords,
                                     const RealVect& a_dx,
                                     const RealVect& a_domainSize,
                                     const Real& a_time, 
                                     const LevelSigmaCS* a_crseCoords,
                                     const int a_refRatio);
   
  bool regridIceGeometry(LevelSigmaCS& a_coords,
                         const RealVect& a_dx,
                         const RealVect& a_domainSize,
                         const Real& a_time,
                         const LevelSigmaCS* a_crseCoords,
                         const int a_refRatio);

  virtual 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);


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


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


  virtual RefCountedPtr<CompGridVTOBC> velocitySolveBC();

  virtual void modifyVelocityRHS(LevelData<FArrayBox>& a_rhs,
                                 LevelSigmaCS& a_coords,
                                 const ProblemDomain& a_domain,
                                 Real a_time, Real a_dt)
  {;}


  virtual void setSurfaceHeightBCs(LevelData<FArrayBox>& a_zSurface,
                                   LevelSigmaCS& a_coords,
                                   const ProblemDomain& a_domain,
                                   const RealVect& a_dx, 
                                   Real a_time, Real a_dt);

  virtual void setGeometryBCs(LevelSigmaCS& a_coords,
                              const ProblemDomain& a_domain,
                              const RealVect& a_dx,
                              Real a_time, Real a_dt);
private:
  
  // Disallowed for all the usual reasons
  void operator=(const MultiLevelDataIBC& a_input)
  {
    MayDay::Error("invalid operator");
  }

  // Disallowed for all the usual reasons
  MultiLevelDataIBC(const MultiLevelDataIBC& a_input)
  {
    MayDay::Error("invalid operator");
  }

private:
  RefCountedPtr<CompGridVTOBC> m_velBCs;
  bool m_isBCsetUp, m_verbose;
  //levels of thickness data
  Vector<RefCountedPtr<LevelData<FArrayBox> > > m_thck;
  //levels of topography data
  Vector<RefCountedPtr<LevelData<FArrayBox> > > m_topg;
  // count of levels, should be == m_thck.size() == m_topg.size()
  int m_nLevel;
  //spacing of the coarsest mesh
  const RealVect m_dxCrse;
  //refinement ratios of the usual kind
  const Vector<int>  m_refRatio;

};
#include "NamespaceFooter.H"
#endif