PISM, A Parallel Ice Sheet Model stable 0.4.1779

src/base/basal_strength/PISMYieldStress.cc File Reference

Process model which computes pseudo-plastic yield stress for the subglacial layer. More...

#include "PISMYieldStress.hh"
#include "Mask.hh"

Go to the source code of this file.


Detailed Description

Process model which computes pseudo-plastic yield stress for the subglacial layer.

The output variable of this submodel is tauc, the pseudo-plastic yield stress field that is used in the ShallowStressBalance (=SSA in 2011) objects.

In the default implementation PISMDefaultYieldStress [BBssasliding], the "dry" strength of the (notional) till is a state variable, but private to the submodel, namely tillphi.

Its initialization is nontrivial: either -topg_to_phi heuristic or inverse modeling (to be implemented ...). Currently tillphi does not evolve during the run.

This submodel uses vHmelt = bwat as an input at each update. Basal water pressure bwp is computed from vHmelt, then that pressure is combined with tillphi to compute an updated tauc by the Mohr-Coulomb criterion.

This submodel is inactive in floating areas.

Definition in file PISMYieldStress.cc.

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines