Difference between revisions of "Kees' assignment"
From Interactive System for Ice sheet Simulation
(→Numerical tips) |
|||
Line 12: | Line 12: | ||
==Model parameters== | ==Model parameters== | ||
*<math>\frac{\partial b}{\partial x} = -0.1</math> | *<math>\frac{\partial b}{\partial x} = -0.1</math> | ||
− | * <math>M(x) = M_0 - x M_1 = 4.0 | + | * <math>M(x) = M_0 - x M_1 = 4.0\times 10^{-3} - 0.2 x</math> km/yr |
* <math>\rho</math> = 920 <math>kg/m^3</math> | * <math>\rho</math> = 920 <math>kg/m^3</math> | ||
*g=9.8 <math>m/s^2</math> | *g=9.8 <math>m/s^2</math> |
Revision as of 13:08, 5 August 2009
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