 |
| Titel |
Numerical ice-sheet model based on the SIA-I algorithm |
| VerfasserIn |
O. Souček, Z. Martinec |
| Konferenz |
EGU General Assembly 2009
|
| Medientyp |
Artikel
|
| Sprache |
Englisch
|
| Digitales Dokument |
PDF |
| Erschienen |
In: GRA - Volume 11 (2009) |
| Datensatznummer |
250028656
|
|
|
|
|
|
| Zusammenfassung |
| We present the numerical results of performance tests of a new numerical model for ice-sheet flow. The geometry evolution of the free surface is modeled by the essentially non-oscillatory (ENO) interpolation schemes with glacier boundaries tracked by the level-set method. The ice velocities and induced viscous stresses are computed by the SIA-I algorithm, developed by the authors, that iteratively improves the shallow-ice approximation. All governing equations are discretized on a non-uniform grid enabling locally increase a model resolution in areas of interest.
We present the results of model output for the ISMIP-HOM benchmark experiment F, that is for a 3D prognostic run. This numerical experiment is additionally extended for a non-linear ice rheology. We compare the results provided by our code with FEM simulations implemented in the Elmer software and show a very good numerical performance of our code with respect to both the accuracy and computational speed. This indicates good applicability of the SIA-I algorithm in complex evolutionary numerical glaciological models. |
|
|
|
|
|
|