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

◆ tweak()

void pism::bed::LingleClarkSerial::tweak ( petsc::Vec load_thickness,
petsc::Vec U,
int  Nx,
int  Ny,
double  time 
)
private

Modify the plate displacement to correct for the effect of imposing periodic boundary conditions at a finite distance.

See Section 5 in [BuelerLingleBrown].

Parameters
[in]load_thicknessthickness of the load (used to compute the corresponding disc volume)
[in,out]Uviscous plate displacement
[in]Nxgrid size
[in]Nygrid size
[in]timetime, seconds (usually 0 or a large number approximating \infty)

Definition at line 512 of file LingleClarkSerial.cc.

References m_D, m_dx, m_dy, m_eta, m_load_density, m_Lx, m_Ly, m_mantle_density, m_standard_gravity, PISM_CHK, and pism::bed::viscDisc().

Referenced by step(), and uplift_problem().