Modeling melange back-pressure
Equation (33) above, describing the stress boundary condition for ice shelves,
can be written in terms of velocity components:
(38)
Here is the vertically-averaged ice viscosity, is the ice thickness, is the
elevation of the bottom and of the top ice surface, and are pressures
of the column of ice and water, respectively:
(39)
We call the integral on the right hand side of (38) the pressure difference
term.
It can be re-written as
(40)
PISM’s ocean model components provide , the vertically-averaged
pressure of the water column adjacent to an ice margin.
To model the effect of melange [105] on the stress boundary condition
we modify the pressure difference term in (38), adding , the
vertically-averaged melange back pressure:
(41)
By default, is zero, but PISM implements two ocean model components to support
scalar time-dependent melange pressure forcing. Please see the Climate Forcing
Manual for details.