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

◆ 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 251 of file LingleClarkSerial.cc.

References C, 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().