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