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