tantspin-coarse.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 --- tantspin-coarse.sh (4566B) --- 1 #!/bin/bash 2 3 # Copyright (C) 2009-2018 PISM authors 4 ################################################################################## 5 # Coarse grid spinup of Antarctic ice sheet model using data from Anne Le Brocq 6 # (see SeaRISE wiki). Uses PIK physics and enthalpy model 7 # and modified configuration parameters with constant climate. Uses constant 8 # and precip and a parameterization for artm as in Martin et al (2011). 9 # WARNING: at finer resolutions (e.g. under 15 km), output is large! 10 ################################################################################## 11 12 SCRIPTNAME="#(antspin-coarse.sh)" 13 14 echo "$SCRIPTNAME run preprocess.sh before this..." 15 # needed inputs generated by preprocess.sh: pism_Antarctica_5km.nc 16 17 set -e # exit on error 18 19 echo "$SCRIPTNAME Coarse grid, constant-climate spinup script using SeaRISE-Antarctica data" 20 echo "$SCRIPTNAME and -stress_balance ssa+sia and -pik" 21 echo "$SCRIPTNAME Run as './antspinCC.sh NN' for NN procs and 30km grid" 22 23 # naming files, directories, executables 24 RESDIR= 25 BOOTDIR= 26 PISM_EXEC=${PISM_EXEC:=pismr} 27 PISM_MPIDO="mpiexec -n " 28 29 # input data: 30 export PISM_INDATANAME=${BOOTDIR}pism_Antarctica_5km.nc 31 32 source set-physics.sh 33 34 NN=4 # default number of processors 35 if [ $# -gt 0 ] ; then # if user says "antspinCC.sh 8" then NN = 8 36 NN="$1" 37 fi 38 echo "$SCRIPTNAME NN = $NN" 39 set -e # exit on error 40 41 # check if env var PISM_DO was set (i.e. PISM_DO=echo for a 'dry' run) 42 if [ -n "${PISM_DO:+1}" ] ; then # check if env var is already set 43 echo "$SCRIPTNAME PISM_DO = $PISM_DO (already set)" 44 else 45 PISM_DO="" 46 fi 47 DO=$PISM_DO 48 49 # these coarse grid defaults are for development/regression runs, not 50 # "production" or science 51 GRID=$THIRTYKMGRID 52 SKIP=$SKIPTHIRTYKM 53 GRIDNAME=30km 54 55 echo "$SCRIPTNAME PISM = $PISM_EXEC" 56 echo "$SCRIPTNAME FULLPHYS = $FULLPHYS" 57 echo "$SCRIPTNAME PIKPHYS = $PIKPHYS" 58 echo "$SCRIPTNAME PIKPHYS_COUPLING = $PIKPHYS_COUPLING" 59 60 61 # ####################################### 62 # bootstrap and SHORT smoothing run to 100 years 63 # ####################################### 64 stage=earlyone 65 INNAME=$PISM_INDATANAME 66 RESNAMEONE=${RESDIR}${stage}_${GRIDNAME}.nc 67 RUNTIME=1 68 echo 69 echo "$SCRIPTNAME bootstrapping on $GRIDNAME grid plus SIA run for $RUNTIME a" 70 cmd="$PISM_MPIDO $NN $PISM_EXEC -skip -skip_max $SKIP -i ${INNAME} -bootstrap $GRID \ 71 $SIA_ENHANCEMENT $PIKPHYS_COUPLING -front_retreat_file ${INNAME} \ 72 -y $RUNTIME -o $RESNAMEONE" 73 $DO $cmd 74 #exit # <-- uncomment to stop here 75 76 stage=smoothing 77 RESNAME=${RESDIR}${stage}_${GRIDNAME}.nc 78 RUNTIME=100 79 echo 80 echo "$SCRIPTNAME short SIA run for $RUNTIME a" 81 cmd="$PISM_MPIDO $NN $PISM_EXEC -skip -skip_max $SKIP -i $RESNAMEONE \ 82 $SIA_ENHANCEMENT $PIKPHYS_COUPLING -front_retreat_file ${INNAME} \ 83 -y $RUNTIME -o $RESNAME" 84 $DO $cmd 85 86 # ####################################### 87 # run with -no_mass (no surface change) on coarse grid for 200ka 88 # ####################################### 89 stage=nomass 90 INNAME=$RESNAME 91 RESNAME=${RESDIR}${stage}_${GRIDNAME}.nc 92 TSNAME=${RESDIR}ts_${stage}_${GRIDNAME}.nc 93 RUNTIME=200000 94 EXTRANAME=${RESDIR}extra_${stage}_${GRIDNAME}.nc 95 expackage="-extra_times 0:1000:$RUNTIME -extra_vars bmelt,tillwat,velsurf_mag,temppabase,diffusivity,hardav" 96 echo 97 echo "$SCRIPTNAME -no_mass (no surface change) SIA for $RUNTIME a" 98 cmd="$PISM_MPIDO $NN $PISM_EXEC -i $INNAME $PIKPHYS_COUPLING \ 99 $SIA_ENHANCEMENT -no_mass \ 100 -ys 0 -y $RUNTIME \ 101 -extra_file $EXTRANAME $expackage \ 102 -o $RESNAME" 103 $DO $cmd 104 #exit # <-- uncomment to stop here 105 106 107 # ####################################### 108 # run into steady state with constant climate forcing 109 # ####################################### 110 stage=run 111 INNAME=$RESNAME 112 RESNAME=${RESDIR}${stage}_${GRIDNAME}.nc 113 TSNAME=${RESDIR}ts_${stage}_${GRIDNAME}.nc 114 RUNTIME=100000 115 EXTRANAME=${RESDIR}extra_${stage}_${GRIDNAME}.nc 116 exvars="thk,usurf,velbase_mag,velbar_mag,mask,diffusivity,tauc,bmelt,tillwat,temppabase,hardav,cell_grounded_fraction,ice_area_specific_volume,amount_fluxes,basal_mass_flux_grounded,basal_mass_flux_floating" 117 expackage="-extra_times 0:1000:$RUNTIME -extra_vars $exvars" 118 119 echo 120 echo "$SCRIPTNAME run into steady state with constant climate forcing for $RUNTIME a" 121 cmd="$PISM_MPIDO $NN $PISM_EXEC -bootstrap -Lz 5000 -skip -skip_max $SKIP -i $INNAME \ 122 $SIA_ENHANCEMENT $PIKPHYS_COUPLING $PIKPHYS $FULLPHYS \ 123 -ys 0 -y $RUNTIME \ 124 -ts_file $TSNAME -ts_times 0:1:$RUNTIME \ 125 -extra_file $EXTRANAME $expackage \ 126 -o $RESNAME -o_size big" 127 $DO $cmd 128 129 echo 130 echo "$SCRIPTNAME coarse-grid part of spinup done"