trun_prog.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 --- trun_prog.sh (2088B) --- 1 #!/bin/bash 2 3 NN=2 # default number of processors 4 if [ $# -gt 0 ] ; then # if user says "run_prog.sh 8" then NN = 8 5 NN="$1" 6 fi 7 8 M=211 # for 5km; note M=526 for 2km, M=421 for 2.5km, M=351 for 3km 9 if [ $# -gt 1 ] ; then # if user says "run_prog.sh 8 211" then NN = 8 and M = 351 10 M="$2" 11 fi 12 13 SSAE=0.6 # ssa enhancement factor 14 if [ $# -gt 2 ] ; then # if user says "run_prog.sh 8 211 0.7" then ... and -ssa_e 0.7 15 SSAE="$3" 16 fi 17 18 YEARS=3000 # integration time 19 if [ $# -gt 3 ] ; then # if user says "run_prog.sh 8 211 0.7 200" then ... and -y 200 20 YEARS="$4" 21 fi 22 23 PISMPREFIX="" 24 #PISMPREFIX="../../../bin/" 25 26 # use these for more robust SSA solves 27 # STRONGKSP="-ssafd_ksp_type gmres -ssafd_ksp_norm_type unpreconditioned -ssafd_ksp_pc_side right -ssafd_pc_type asm -ssafd_sub_pc_type lu" 28 # STRONGKSP="" 29 30 # preliminary bootstrap and diagnostic run: 31 STARTNAME=startfile_Mx${M}.nc 32 cmd_diag="mpiexec -n $NN ${PISMPREFIX}pismr -i ../Ross_combined.nc -bootstrap -Mx $M -My $M \ 33 -Mz 61 -Lz 3000 -z_spacing equal -surface given -stress_balance ssa \ 34 -yield_stress constant -tauc 1e6 -pik -ssa_dirichlet_bc -ssa_e $SSAE \ 35 $STRONGKSP -y 0.1 -ys 0.0 -o $STARTNAME -o_order zyx " 36 echo "running command:" 37 echo 38 echo "$cmd_diag" 39 ${cmd_diag} 40 41 # prognostic run 42 NAME=prog_Mx${M}_yr${YEARS}.nc 43 ECALV=1e18 # constant for eigen_calving parameterization 44 CTHICK=50.0 # constant thickness for thickness_calving 45 exdt=1 46 exvars="thk,mask,velsurf_mag,strain_rates,tendency_of_ice_mass_due_to_flow,tendency_of_ice_mass_due_to_discharge" 47 cmd_prog="mpiexec -n $NN ${PISMPREFIX}pismr -i $STARTNAME \ 48 -surface given -stress_balance ssa -yield_stress constant -tauc 1e6 -pik \ 49 -calving eigen_calving,thickness_calving -eigen_calving_K $ECALV -front_retreat_cfl \ 50 -ssa_dirichlet_bc -ssa_e $SSAE -ys 0 -y $YEARS -o $NAME -o_size big \ 51 -thickness_calving_threshold $CTHICK $STRONGKSP \ 52 -ts_file ts-${NAME} -ts_times 0:monthly:${YEARS} \ 53 -extra_file ex-${NAME} -extra_times 0:${exdt}:${YEARS} -extra_vars ${exvars} \ 54 -options_left" 55 echo "running command:" 56 echo 57 echo "$cmd_prog" 58 ${cmd_prog} 59