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 (1691B)
       ---
            1 # Nbreen hydrology example
            2 
            3 This directory contains draft scripts for a hydrology-only simulation using
            4 ice geometry data for the Nordenskioldbreen glacier on Svalbard.  This example
            5 is a basic test of mass-conservative subglacial hydrology models with cavity
            6 evolution (i.e. cavitation by sliding and cavity closure by creep).  The basal
            7 melt rate input is constant and the sliding is artificially set to 50 m/a
            8 everywhere.
            9 
           10 This example is related to a manuscript in preparation by Ward van Pelt and
           11 Ed Bueler.
           12 
           13 In the future there may be an additional ice dynamics simulation coupled to the
           14 hydrology submodel.
           15 
           16 ## Basic usage
           17 
           18 First do
           19 
           20     ln -s ../../util/PISMNC.py
           21     ./preprocess.sh
           22 
           23 This builds PISM-readable NetCDF file `pismnbreen.nc` and it generates a
           24 synthetic summer event input file `fakesummerevent.nc` (which has a large
           25 pulse of water put into the subglacial system in a few weeks in the middle of
           26 the year).  It also converts the `nbreen_config.cdl` to get `nbreen_config.nc`.
           27 
           28 To see how to run examples, type
           29 
           30     ./run.sh
           31 
           32 This will generate a usage message.
           33 
           34 For a quick first run, try
           35 
           36     ./run.sh 2 500 0.1 routing daily
           37 
           38 This uses 2 processors, a 500 m grid, a 0.1 year duration, the simpler "routing"
           39 model, and monthly reporting.  It should finish in under one minute.  View the
           40 time- and space-dependent output file:
           41 
           42     ncview extras_nbreen_y0.1_500m_routing.nc
           43 
           44 Another visualization is with python,
           45 
           46     python showhydrovel.py -t -s extras_nbreen_y0.1_500m_routing.nc
           47 
           48 This shows a quiver plot of the water velocity, among other things.
           49 
           50 ## Changing configuration constants
           51 
           52 Edit `nbreen_config.cdl` to set your preferred value and rerun `preprocess.sh`.