Processing math: 100%
PISM, A Parallel Ice Sheet Model 2.2.2-d6b3a29ca committed by Constantine Khrulev on 2025-03-28
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages

◆ uplift_problem()

void pism::bed::LingleClarkSerial::uplift_problem ( petsc::Vec load_thickness,
petsc::Vec bed_uplift,
petsc::Vec output 
)
private

Solve

2 \nu |\nabla| \diff{u}{t} + \rho_r g U + D\nabla^4 U = \sigma_{zz}

for U , treating \diff{u}{t} and \sigma_{zz} as known.

Parameters
[in]load_thicknessload thickness, meters
[in]bed_upliftbed uplift, m/second

Here load_thickness is used to compute the load \sigma_{zz} and bed_uplift is \diff{u}{t} itself.

Definition at line 252 of file LingleClarkSerial.cc.

References pism::clear_fftw_array(), pism::copy_fftw_array(), pism::get_real_part(), m_cx, m_cy, m_D, m_dft_forward, m_dft_inverse, m_eta, m_fftw_input, m_fftw_output, m_i0_offset, m_j0_offset, m_load_density, m_loadhat, m_mantle_density, m_Mx, m_My, m_Nx, m_Ny, m_standard_gravity, pism::set_real_part(), and tweak().

Referenced by bootstrap().