BISICLES AMR ice sheet model
0.9
|
Zwinger Rate Factor. More...
#include <ConstitutiveRelation.H>
Public Member Functions | |
ZwingerRateFactor (Real a_seconds_per_unit_time) | |
void | setDefaultParameters (Real a_seconds_per_unit_time) |
void | setParameters (Real a_E, Real a_A0, Real a_T0, Real a_R, Real a_Qm, Real a_Qp) |
void | computeA (FArrayBox &a_A, const FArrayBox &a_thetaPC, const FArrayBox &a_pressure, const Box &a_box) const |
RateFactor * | getNewRateFactor () const |
![]() | |
virtual | ~RateFactor () |
Zwinger Rate Factor.
Another Arrhenius law, this one from Zwinger (2007) Annals of Glaciology vol 45 p 29 and used in (older versions of?) Elmer/Ice.
A = E*A0(T)*exp(-Qc / (R * T) ] where Qc = Qm = if theta < theta0; Qc = Qp otherwise. and A0 = Am = if theta < theta0; A0 = Ap otherwise. Am is computed so that A is continuous with T
Th is the usual pressure-corrected temperature Th = T + 7e−8 * P T0* = T0 + 7e−8 * P (T0 = 263, P = pressure) Qc = Qm = if Th < T∗; Qc = Qp if Th > T∗.
Defaults:
A0 = 1.916e3 Pa−3 s−1 * secondsperyear Qp = 139 kJ mol−1 Qm = 6e4 kJ mol−1 R=8.314 J
ZwingerRateFactor::ZwingerRateFactor | ( | Real | a_seconds_per_unit_time | ) |
|
virtual |
Implements RateFactor.
References IceThermodynamics::icepmeltfactor().
Referenced by setParameters().
|
virtual |
Implements RateFactor.
void ZwingerRateFactor::setDefaultParameters | ( | Real | a_seconds_per_unit_time | ) |
References setParameters().
void ZwingerRateFactor::setParameters | ( | Real | a_E, |
Real | a_A0, | ||
Real | a_T0, | ||
Real | a_R, | ||
Real | a_Qm, | ||
Real | a_Qp | ||
) |
References computeA().
Referenced by setDefaultParameters().