VieliPayneIBC.H

Go to the documentation of this file.

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

//
//  PayneVieliIBC.H
// ============
//
// PhysIBC-derived class for grounding-line testing
//

#ifndef _VIELIPAYNEIBC_H_
#define _VIELIPAYNEIBC_H_

#include "IceThicknessIBC.H"

#include "NamespaceHeader.H"



class VieliPayneBCFunction : public IceVelBCFunction
{
public:
  VieliPayneBCFunction() {;}
  
  VieliPayneBCFunction(Vector<RefCountedPtr<LevelSigmaCS > >& a_vectCS,
                       Vector<ProblemDomain>& a_vectDomain);
  
  
  ~VieliPayneBCFunction(); 

  virtual void define(Vector<RefCountedPtr<LevelSigmaCS > >& a_vectCS,
                      Vector<ProblemDomain>& a_vectDomain);
  
  

  virtual void operator()(FArrayBox&           a_state,
                          const Box&           a_valid,
                          const ProblemDomain& a_domain,
                          Real                 a_dx,
                          bool                 a_homogeneous);

  virtual void operator()(FArrayBox&           a_state,
                          const Box&           a_valid,
                          const ProblemDomain& a_domain,
                          Real                 a_dx,
                          const DataIndex&     a_index,
                          bool                 a_homogeneous);


protected:
  
  Vector<RefCountedPtr<LevelSigmaCS > > m_vectCS;
  Vector<ProblemDomain> m_vectDomain;
};



class VieliPayneIBC : public IceThicknessIBC
{
public:

  VieliPayneIBC();


  virtual ~VieliPayneIBC();


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

  void setParameters(const Real& a_thickness,
                     const RealVect& a_slope,
                     const Real& a_originElevation,
                     const Real& a_seaLevel);


  virtual IceThicknessIBC* new_thicknessIBC();

#if 0
  void setBoundaryThickness(Real a_boundaryThickness)
  { m_boundaryThickness = a_boundaryThickness; }
#endif


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



  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 velocityGhostBC(LevelData<FArrayBox>& a_velocity,
                               LevelSigmaCS& a_coords,
                               const ProblemDomain& a_domain,
                               Real a_time);


  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 setGridHierarchy(Vector<RefCountedPtr<LevelSigmaCS > >& a_vectCS,
                                Vector<ProblemDomain>& a_vectDomain);
    


  static RealVect s_edgeThickness;

protected:

  RealVect m_slope;

  Real m_originElevation;

  Real m_thickness;

  RealVect m_domainSize;

  Real m_seaLevel;

  // have bc's been set up?
  bool m_isBCsetUp;

  // have parameters been set
  bool m_paramsSet;

  void setupBCs();

  RefCountedPtr<CompGridVTOBC> m_velBCs;
  RefCountedPtr<VieliPayneBCFunction> m_BCfunction;


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

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




#include "NamespaceFooter.H"
#endif