# 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