PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800

◆ P_from_W_steady()

void pism::hydrology::Distributed::P_from_W_steady ( const array::Scalar W,
const array::Scalar P_overburden,
const array::Scalar sliding_speed,
array::Scalar result 
)
protected

Compute functional relationship P(W) which applies only in steady state.

In steady state in this model, water pressure is determined by a balance of cavitation (opening) caused by sliding and creep closure.

This will be used in initialization when P is otherwise unknown, and in verification and/or reporting. It is not used during time-dependent model runs. To be more complete, \(P = P(W,P_o,|v_b|)\).

Definition at line 170 of file Distributed.cc.

References pism::Component::m_config, pism::Component::m_grid, pism::array::max(), and pism::Wr.

Referenced by bootstrap_impl().