SurfaceFlux which has a piecewise linear dependence on ice thickness. More...

#include <ComplexSurfaceFlux.H>

Inheritance diagram for PiecewiseLinearFlux:

Collaboration diagram for PiecewiseLinearFlux:

Public Member Functions
	PiecewiseLinearFlux (const Vector< Real > &a_abscissae, const Vector< Real > &a_ordinates, const Real a_minWaterDepth=-1.0)

	~PiecewiseLinearFlux ()

virtual SurfaceFlux *	new_surfaceFlux ()
	factory method: return a pointer to a new SurfaceFlux object More...

virtual void	surfaceThicknessFlux (LevelData< FArrayBox > &a_flux, const AmrIceBase &a_amrIce, int a_level, Real a_dt)
	define source term for thickness evolution and place it in flux More...

Public Member Functions inherited from SurfaceFlux
virtual	~SurfaceFlux ()
	virtual destructor More...

virtual void	evaluate (LevelData< FArrayBox > &a_data, const AmrIceBase &a_amrIce, int a_level, Real a_dt)
	Sub-classes of this class turned out to be useful for various data apart from thickness fluxes... More...

Protected Attributes
Vector< Real >	m_abscissae

Vector< Real >	m_ordinates

Real	m_minWaterDepth

Additional Inherited Members
Static Public Member Functions inherited from SurfaceFlux
static SurfaceFlux *	parse (const char *a_prefix)
	assemble a SurfaceFlux* object from ParmParse inputs More...

Detailed Description

SurfaceFlux which has a piecewise linear dependence on ice thickness.

Usually used in conjunction with MaskedFlux to impose a depth dependent ice shelf melt rate. Provided that the cavity water thickness $ d = s-h-b > d_0$ , the surface flux is computed from $ n $ pairs of abissiciae and ordinates $ ( h_i, f_i ) $

$ f = f_0 $ for $ h < h_0 $

$f = f_{i}(h-h_{i+1})/(h{i}-h_{i+1}) + f_{i+1}(h-h_{i})/(h_{i+1}-h_{i}) +$ for $h_i \leq h < h_{i+1}$

$f = f_{n-1}$ for $h > h_{n-1}$

Run-time configuration is carried out in SurfaceFlux::parse. If the string A.B.C.type = piecewiseLinearFlux is found, then the expected form is:

A.B.C.n = <int> (number of abscissae and ordinates)
A.B.C.abscissae = <float> <float> ... (list of thickness absissae)
A.B.C.ordinates = <float> <float> ... (list of flux ordinates)
`A.B.C.minWaterDepth = <float (optional, default -1.0)>

Constructor & Destructor Documentation

◆ PiecewiseLinearFlux()

PiecewiseLinearFlux::PiecewiseLinearFlux	(	const Vector< Real > &	a_abscissae,
		const Vector< Real > &	a_ordinates,
		const Real	a_minWaterDepth = `-1.0`
	)

Parameters

a_abscissae	list of thickness abscissae
a_abscissae	list of flux ordinates
a_minWaterDepth	set flux to zero in sub-shelf cavities thinner than a_minWaterDepth