trun_diag.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_diag.sh (979B)
       ---
            1 #!/bin/bash
            2 
            3 # '-pik' sets the calving front boundary condition
            4 #     (is equivalent to '-cfbc -kill_icebergs')
            5 # '-ssa_method fem' will not work because it lacks '-cfbc'
            6 
            7 NN=2  # default number of processors
            8 if [ $# -gt 0 ] ; then  # if user says "run_diag.sh 8" then NN = 8
            9   NN="$1"
           10 fi
           11 
           12 M=211   # for 5km;  note M=526 for 2km,  M=421 for 2.5km,  M=351 for 3km
           13 if [ $# -gt 1 ] ; then  # if user says "run_diag.sh 8 351" then NN = 8 and M = 351
           14   M="$2"
           15 fi
           16 
           17 SSAE=0.6
           18 if [ $# -gt 2 ] ; then  # if user says "run_diag.sh 8 351 0.7" then ... and '-ssa_e 0.7'
           19   SSAE="$3"
           20 fi
           21 
           22 PISMPREFIX=""
           23 #PISMPREFIX="../../../bin/"
           24 
           25 cmd="mpiexec -n $NN ${PISMPREFIX}pismr -i ../Ross_combined.nc -bootstrap -Mx $M -My $M \
           26   -Mz 3 -Lz 3000 -z_spacing equal -surface given -stress_balance ssa -energy none -no_mass \
           27   -yield_stress constant -tauc 1e6 -pik -ssa_dirichlet_bc \
           28   -y 1.0 -o diag_Mx${M}.nc -o_order zyx -ssa_e $SSAE -ssafd_ksp_monitor"
           29 
           30 echo "running command:"
           31 echo
           32 echo "$cmd"
           33 $cmd
           34