bootstrap compute_elastic_response compute_load_response_matrix elastic_displacement init LingleClarkSerial m_cx m_cy m_D m_dft_forward m_dft_inverse m_dx m_dy m_eta m_fftw_input m_fftw_output m_i0_offset m_include_elastic m_j0_offset m_load_density m_loadhat m_log m_lrm_hat m_Lx m_Ly m_mantle_density m_Mx m_My m_Nx m_Ny m_standard_gravity m_U m_Ue m_Uv precompute_coefficients step total_displacement tweak update_displacement uplift_problem viscous_displacement ~LingleClarkSerial

bootstrap() void pism::bed::LingleClarkSerial::bootstrap ( petsc::Vec &  thickness, petsc::Vec &  uplift  ) Initialize using provided load thickness and the bed uplift rate. Here we solve: rho_r g U + D grad^4 U = -rho g H - 2 eta |grad| uplift Compare equation (16) in Bueler, Lingle, Brown (2007) "Fast computation of a viscoelastic deformable Earth model for ice sheet simulations", Ann. Glaciol. 46, 97–105 Note Probable sign error in eqn (16)?: load "rho g H" should be "- rho g H"] This initialization method is used to "bootstrap" the model. It should not be used to re-start a stopped modeling run. Parameters [in] thickness load thickness, meters [in] uplift initial bed uplift on the PISM grid Sets m_Uv, m_Ue, m_U. Definition at line 317 of file LingleClarkSerial.cc. References compute_elastic_response(), m_include_elastic, m_U, m_Ue, m_Uv, PISM_CHK, update_displacement(), and uplift_problem().