tREADME.md - 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.md (4039B)
---
1 ## "Standard" Greenland spinup examples
2
3 This directory contains scripts for standard Greenland ice dynamics simulations.
4 The first section of the PISM _User's Manual_ documents these examples.
5
6 These runs are motivated by (but not identical to) the "constant-climate" and
7 "paleo-climate" initial states which are evaluated in
8
9 * A. Aschwanden, G. Adalgeirsdottir, and C. Khroulev (2013). _Hindcasting to_
10 _measure ice sheet model sensitivity to initial states._ The Cryosphere 7,
11 doi:10.5194/tc-7-1083-2013, 1083--1093.
12
13 The runs here use 5km SeaRISE data including present-day bedrock
14 topography, ice thickness, and basal heat flux; see references at bottom. Also,
15 in contrast to the above _Hindcasting_ paper, runs here use publicly-available
16 RACMO atmospheric/ocean circulation model outputs, including modeled present-day
17 2m air temperature, precipitation, and surface mass balance, as distributed by
18 the SeaRISE project; see reference at bottom.
19
20 ## Basic usage
21
22 First do
23
24 $ ./preprocess.sh
25
26 This downloads the file `Greenland_5km_v1.1.nc` from the SeaRISE site and it
27 builds PISM-readable NetCDF file `pism_Greenland_5km_v1.1.nc`. It also
28 generates additional climate-forcing files `pism_dT.nc` and `pism_dSL.nc`.
29
30 Then do
31
32 $ ./spinup.sh
33
34 *In fact this does no actual simulation.* Instead it provides a usage/help
35 message. Specifically, it shows two basic usages, one with constant-climate
36 and SIA-only choices, and one with paleo-climate and SIA/SSA-hybrid choices.
37 The first of these examples can be run in a few minutes of computer time,
38 which should help you get started.
39
40 # SeaRISE Greenland background information
41
42 The SeaRISE assessment process included PISM Greenland simulations among many
43 ice sheet model results. These runs are described in the following papers,
44 which evaluate and compare all SeaRISE simulations:
45
46 * S. Nowicki and others (2013). _Insights into spatial sensitivities of ice_
47 _mass response to environmental change from the SeaRISE ice sheet modeling_
48 _project: II. Greenland._ J. Geophys. Res.: Earth Surface 118 (2),
49 doi:10.1002/jgrf.20076, 1025--1044.
50
51 * R. Bindshadler and others (2013). _Ice-sheet model sensitivities to_
52 _environmental forcing and their use in projecting future sea-level_
53 _(The SeaRISE Project)._ J. Glaciol. 59 (214), 195--224.
54
55 The reference for RACMO model outputs is:
56
57 * J. Ettema, M.R. van den Broeke, E. van Meigaard, W.J. van de Berg,
58 J.L. Bamber, J.E. Box, and R.C. Bales (2009). _Higher surface mass balance_
59 _of the Greenland ice sheet revealed by high-resolution climate modeling._
60 Geophys. Res. Lett. 36 (L12501), doi:10.1029/2009GL038110.
61
62 Note that the RACMO surface mass balance is missing where there is no
63 (present-day) ice.
64
65 For further references on the SeaRISE data, run
66
67 $ ncdump -h Greenland_5km_v1.1.nc |less
68
69 after running `preprocess.sh` above, and look at the `source` attributes of
70 various fields.
71
72 Scripts for the actual PISM submitted runs to SeaRISE are more complex than
73 the runs shown here. They can be found in the `stable0.5` branch of PISM. See
74
75 https://github.com/pism/pism/tree/stable0.5/examples/searise-greenland
76
77 # Basic Parameter study
78
79 An example for a basic parameter study is given by running the script `param.sh`
80 and either `runsequential.sh` or `runparallel.sh`. See the User's Manual. Also
81 see scripts in `advanced/`.
82
83 # Visualization
84
85 Results can be visualized using, e.g. `im-plot.py` from
86 [pypismtools](https://github.com/pism/pypismtools):
87
88 im-plot.py -v velsurf_mag --colorbar_label --inner_titles "q=0.1 f=0.01,q=0.25 f=0.01,q=0.8 f=0.01,q=0.1 f=0.02,q=0.25 f=0.02,q=0.1 f=0.02,q=0.8 f=0.0,q=0.1 f=0.05,q=0.25 f=0.05,q=0.8 f=0.05" -o nosgl.pdf g20km_*0.01_*.nc g20km_*0.02_*.nc g20km_*0.05_*.nc
89
90 im-plot.py -v velsurf_mag --colorbar_label --inner_titles "q=0.1 f=0.01,q=0.25 f=0.01,q=0.8 f=0.01,q=0.1 f=0.02,q=0.25 f=0.02,q=0.1 f=0.02,q=0.8 f=0.0,q=0.1 f=0.05,q=0.25 f=0.05,q=0.8 f=0.05" -o sgl.pdf g20km_*0.01.nc g20km_*0.02.nc g20km_*0.05.nc
91