tREADME.rst - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
(HTM) git clone git://src.adamsgaard.dk/pism
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) LICENSE
---
tREADME.rst (1385B)
---
1 Example: relax Greenland's basal topograpy
2 ======================
3
4 This minimal example shows how we can calculate the isostatically-relaxed
5 basal topography for an ice-free Greenland. The same method can be
6 applied to other areas. It uses the flux correction method
7 "force-to-thickness" with a target thickness of zero to remove all ice
8 combined with the Lingle-Clark bed deformation model to calculate
9 isostatic uplift.
10
11 How to run it
12 -------------------------
13
14 First, download the SeaRISE data set running
15
16 ``$ ./preprocess.sh``
17
18 or, alternatively, provide a data set of your choice. For example, you can
19 make a link to the SeaRISE data set downloaded into
20 ``examples/searise-greenland/``, if the example in Chapter 1 of the
21 User's Manual has already run.
22
23 Then run
24
25 ``$ ./run-relax.sh N M``
26
27 where N is the number of processors and M defines the grid resolution
28 (0: 20km [default], 1: 10km, 2: 5km, 3: 2.5km, 4: 2km, 5: 1km). To run the
29 example on 16 cores and 5km horizontal grid resolution, in the background,
30 and store ``stdout`` output in a text file, type
31
32 ``$ ./run-relax.sh 16 2 &> out.relax &``
33
34 First, the model is run for 1,000 years using flux-correction to
35 remove (most) of the ice. The second simulation continues from the
36 first, but now with the option -no_mass, for another 50,000 years
37 (Note: this number was defined ad-hoc, and may or may not be appropriate).