troutingtest.sh - 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
       ---
       troutingtest.sh (1837B)
       ---
            1 #!/bin/bash
            2 
            3 #(routingtest.sh) test 'routing' hydrology on Antarctic geometry and constant
            4 #  basal melt rate of 1 cm/year; compare
            5 #  http://www2.gi.alaska.edu/snowice/glaciers/iceflow/bueler-igs-fairbanks-june2012.pdf
            6 
            7 # run preprocess.sh before this script
            8 
            9 NN=4
           10 
           11 DOIT=
           12 #DOIT=echo
           13 
           14 PREFIX=
           15 PISMGO="mpiexec -n $NN ${PREFIX}pismr"
           16 
           17 VERTGRID="-Lz 5000 -Lbz 2000 -Mz 51 -Mbz 21 -z_spacing equal"
           18 
           19 OPTIONS="-sia_e 5.6 -atmosphere given -atmosphere_given_file pism_Antarctica_5km.nc -surface simple -ocean pik -meltfactor_pik 1.5e-2"
           20 
           21 HYDRO="-hydrology routing -hydrology_use_const_bmelt -hydrology_const_bmelt 3.1689e-10 -hydrology_hydraulic_conductivity 1.0e-3 -report_mass_accounting"
           22 
           23 ENDTIME=20000
           24 
           25 dorun () {
           26   GRID=$1
           27   LABEL=$2
           28 
           29   #(from antspinCC.sh)  bootstrapping plus short SIA run for 100 years
           30   cmd="$PISMGO -skip -skip_max 10 -i pism_Antarctica_5km.nc -bootstrap $GRID $VERTGRID $OPTIONS -front_retreat_file pism_Antarctica_5km.nc -y 100 -o pre${LABEL}.nc"
           31   $DOIT $cmd
           32 
           33   EXTRA="-extra_file ex_routing${LABEL}.nc -extra_times 200:100:$ENDTIME -extra_vars bwat,bwp,bwatvel,subglacial_water_input_rate"
           34 
           35   #hydrology only run for $ENDTIME years
           36   cmd="$PISMGO -i pre${LABEL}.nc -bootstrap -Lz 5000 $OPTIONS -front_retreat_file pism_Antarctica_5km.nc $HYDRO -no_mass -energy none -stress_balance ssa -max_dt 10.0 -ys 0 -ye $ENDTIME $EXTRA -o routing${LABEL}.nc"
           37   $DOIT $cmd
           38 }
           39 
           40 HUNDREDKMGRID="-Mx 60 -My 60"
           41 FIFTYKMGRID="-Mx 120 -My 120"
           42 TWENTYFIVEKMGRID="-Mx 240 -My 240"
           43 FIFTEENKMGRID="-Mx 400 -My 400"
           44 TENKMGRID="-Mx 600 -My 600"
           45 FIVEKMGRID="-Mx 1200 -My 1200"
           46 
           47 # these first three regenerate results from IGS talk:
           48 #dorun "$HUNDREDKMGRID" 100km
           49 #dorun "$FIFTYKMGRID" 50km
           50 dorun "$TWENTYFIVEKMGRID" 25km
           51 
           52 # these are more expensive, naturally
           53 #dorun "$FIFTEENKMGRID" 15km
           54 #dorun "$TENKMGRID" 10km
           55 #dorun "$FIVEKMGRID" 5km
           56