singularStreamFriction.H

Go to the documentation of this file.

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

// Virtual base class encapsulating coefficient of basal friction 
// for ice-sheet model

#ifndef _SINGULARSTREAMFRICTION_H_
#define _SINGULARSTREAMFRICTION_H_

#include "BasalFriction.H"
#include "LevelData.H"
#include "FArrayBox.H"
#include "LevelSigmaCS.H"

#include "NamespaceHeader.H"


class singularStreamFriction : public BasalFriction
{
public:
  singularStreamFriction(const Real& a_slippyC,
                         const Real& a_stickyC,
                         const Real& a_width,
                         const Real& a_twistNumber,
                         const Real& a_twistAmplitude,
                         const RealVect& a_domainSize);
    
  virtual ~singularStreamFriction();


  virtual BasalFriction* new_basalFriction() const;


  virtual void setBasalFriction(LevelData<FArrayBox>& a_C,
                                LevelSigmaCS& a_coordSys,
                                Real a_time,
                                Real a_dt);

  void setParameters(const Real& a_slippyC,
                     const Real& a_stickyC,
                     const Real& a_width,
                     const Real& a_twistNumber,
                     const Real& a_twistAmplitude,
                     const RealVect& a_domainSize);

  Real m_slippyC, m_stickyC, m_width, m_twistNumber, m_twistAmplitude;
  RealVect m_domainSize;

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

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


#include "NamespaceFooter.H"
#endif