trungum.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
       ---
       trungum.sh (1330B)
       ---
            1 #!/bin/bash
            2 
            3 if [ $# -lt 2 ] ; then
            4   echo "rungum.sh ERROR: needs two arguments ... ENDING NOW"
            5   echo "  format:"
            6   echo "    rungum.sh PROCS MX"
            7   echo "  where"
            8   echo "    PROCS     = 1,2,3,... is number of MPI processes"
            9   echo "    MX        = number of grid points in x,y directions;  MX -> cell width:"
           10   echo "                52 -> 10mm,  104 -> 5mm, 208 -> 2.5mm, 520 -> 1mm"
           11   exit
           12 fi
           13 
           14 NN="$1"
           15 myMx="$2"
           16 
           17 pismexec="pismr"
           18 
           19 oname=lab$myMx.nc
           20 initfile=init$oname
           21 
           22 grid="-Mx $myMx -My $myMx -Mz 101 -Lz 0.1 -Mbz 0 -Lbz 0 -z_spacing equal"
           23 
           24 climate="-surface given -surface_given_file $initfile -surface.given.smb_max 1e10"
           25 
           26 physics="-config_override gumparams.nc -energy none -sia_flow_law isothermal_glen -sia_e 1.0 -gradient mahaffy"
           27 
           28 endtime=2.36555048199e-5   # = 746 / (365*24*60*60); = 746 s in years;
           29                            # Sayag personal communication
           30 
           31 ts_dt=3.1689e-09     # = 0.1 / 31556926 = 0.1 s
           32 timediag="-ts_file ts_$oname -ts_times $ts_dt:$ts_dt:$endtime"
           33 
           34 ex_dt=3.1689e-07     # = 10 / 31556926 = 10 s
           35 exvars="diffusivity,flux_mag,velbar_mag,velsurf_mag,mask,thk,wvelsurf"
           36 exdiag="-extra_file ex_$oname -extra_vars $exvars -extra_times 0:$ex_dt:$endtime"
           37 
           38 mpiexec -n $NN $pismexec -i $initfile -bootstrap $grid $climate $physics \
           39     $timediag $exdiag -ys 0.0 -y $endtime -max_dt $ts_dt -o $oname