tREADME.md - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
HTML git clone git://src.adamsgaard.dk/pism
DIR Log
DIR Files
DIR Refs
DIR LICENSE
---
tREADME.md (1152B)
---
1 ### A simplified setup showing that `-kill_icebergs` should always be "on".
2
3 We start with three blobs of ice on a flat bed. As the run goes on,
4 the sea level rises, and blobs start floating, starting with the
5 smallest one.
6
7 Once the whole blob is afloat the SSA solver will fail because of a
8 zero pivot. The option `-kill_icebergs` removes such disconnected floating ice.
9
10 To run this example, first put `util/PISMNC.py` on your `PYTHONPATH`, or do
11
12 $ ln -s ../../../util/PISMNC.py
13
14 Now run all examples
15
16 $ make
17
18 First look at variable `thk` in `flood.nc` to see initial geometry.
19 See evolving `mask` and `thk` variables in output files `ex_*.nc` to see
20 the effect of various options including `[none]`, `-cfbc`,
21 `-cfbc -kill_icebergs`.
22
23 Note that without the calving front boundary condition (`make no_cfbc`
24 which uses `[none]`), the strength extension "ties" floating blobs to the
25 biggest blob, which remains grounded. In this (not-recommended) case the
26 SSA solver does not see a zero pivot.
27
28 As a result of this example and others, the option combination
29 `-cfbc -kill_icebergs` is generally recommended for marine ice sheet
30 simulations.