Kees' assignment
From Interactive System for Ice sheet Simulation
Contents |
Model equation
where
and
Model parameters
-
km/yr
-
= 920
- g=9.8
- A = 1e-16
- n=3
Boundary conditions
-
(left boundary)
-
(right boundary)
Numerical tips
Use a staggered grid such that the are computed at the centers of the grid (as opposed to the vertices, as we have been doing), so
From the diffusivity, the flux is computed
,
where
and then the flux () can be used to compute the rate of change of the surface from