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