PISM, A Parallel Ice Sheet Model  stable v2.0.4 committed by Constantine Khrulev on 2022-05-25 12:02:27 -0800

◆ define_model_state_impl()

void pism::MohrCoulombYieldStress::define_model_state_impl ( const File output) const
protectedvirtual

Define model state variables.

All yield stress models have to write basal yield stress to output files and read it from and input file during initialization because yield stress may be used by PISM's stress balance model. The stress balance code has to be executed early during an update of the model because its output (ice velocity) is used to compute the maximum allowed time step.

Now that PISM's yield stress models are time-dependent YieldStress::update() will be called after the maximum time step is found. This means that during the first time step basal_material_yield_stress() gets called before update().

Reimplemented from pism::YieldStress.

Reimplemented in pism::OptTillphiYieldStress.

Definition at line 224 of file MohrCoulombYieldStress.cc.

References pism::IceModelVec::define(), pism::YieldStress::m_basal_yield_stress, and m_till_phi.

Referenced by pism::OptTillphiYieldStress::define_model_state_impl().