# Difference between revisions of "Kees' assignment"

From Interactive System for Ice sheet Simulation

(added solution section) |
|||

Line 42: | Line 42: | ||

:<math>\frac{\partial H}{\partial t} = -\frac{\partial }{\partial x} H\bar u + M = - \frac{\phi_{i+1/2} - \phi_{i-1/2} }{\Delta x} + M</math> | :<math>\frac{\partial H}{\partial t} = -\frac{\partial }{\partial x} H\bar u + M = - \frac{\phi_{i+1/2} - \phi_{i-1/2} }{\Delta x} + M</math> | ||

+ | |||

+ | ==Solutions== | ||

+ | * [[Team 1 Solution]] | ||

+ | * [[Team 2 Solution]] | ||

+ | * [[Team 3 Solution]] | ||

+ | * [[Team 4 Solution]] | ||

+ | * [[Team 5 Solution]] | ||

+ | * [[Team 6 Solution]] |

## Revision as of 15:04, 5 August 2009

## Contents |

## Model equation

where

and

## Model parameters

- km/yr
- = 920
- g=9.8
- A = 1e-16
- n=3
- dx=1.0 km

## 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