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

◆ viscDisc()

double pism::bed::viscDisc ( double  t,
double  H0,
double  R0,
double  r,
double  rho,
double  rho_ice,
double  grav,
double  D,
double  eta 
)

Actually compute the response of the viscous half-space model in LingleClark, to a disc load.

Compute the response of the viscous half-space model in [LingleClark] to a disc load.

Parameters
[in]ttime in seconds
[in]H0thickness of the disc load, meters
[in]R0radius of the disc load, meters
[in]rradius, meters
[in]rhomantle density, kg/m3
[in]rho_iceice (load) density, kg/m3
[in]gravacceleration due to gravity, m/s2
[in]Dlithosphere flexural rigidity, N meter
[in]etamantle viscosity, Pascal second

Definition at line 135 of file greens.cc.

References pism::F(), grav, H0, rho, rho_ice, pism::array::sum(), and vd_integrand().

Referenced by pism::bed::LingleClarkSerial::tweak().