PISM, A Parallel Ice Sheet Model
stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800
|
◆ define_model_state_impl()
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. Definition at line 101 of file RegionalYieldStress.cc. References pism::array::Array::define(), pism::YieldStress::m_basal_yield_stress, m_input, and pism::io::PISM_DOUBLE. |