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).