tissue-74-test.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 --- tissue-74-test.sh (1342B) --- 1 #!/bin/bash 2 3 # This script uses the circular very similar to the one in run_noshelf.sh to confirm the 4 # issue #74 (https://github.com/pism/pism/issues/74), i.e. that when -calving float_kill 5 # is used the sub-shelf ice flux reported by PISM (the variable basal_mass_flux_floating) 6 # can be non-zero even when the total area of the floating ice (variable 7 # ice_area_glacierized_floating) is zero. 8 9 # This is due to the fact that during time-stepping the sub-shelf ice flux is computed 10 # before calving is applied. This error has an O(dt) character. This will be fixed once 11 # the 2D mass transport code gets an overhaul, but for now basal_mass_flux_floating should 12 # be attributed to calving. (This is consistent with what the model does right now.) 13 14 set -x 15 set -e 16 17 N=4 18 xx=101 19 yy=$xx 20 length=50 21 22 infile="circular_noshelf.nc" 23 24 ./circular_dirichlet.py -Mx $xx -My $yy -o $infile 25 26 grid="-Mx $xx -My $yy -Mz 31 -Mbz 1 -Lz 1500 -Lbz 0" 27 28 pismopts="-i $infile -bootstrap $grid -stress_balance ssa -ssa_dirichlet_bc -o_order zyx -energy none -ssa_method fd -cfbc -part_grid" 29 30 doit="mpiexec -n $N pismr $pismopts" 31 32 extra="-extra_times 1 -extra_vars thk,mask,velbar_mag,ice_area_specific_volume,velbar,usurf,mass_fluxes -extra_file issue-74_ex.nc" 33 ts="-ts_file issue-74_ts.nc -ts_times 1" 34 35 $doit -y $length -o issue-74_o.nc $extra $ts -calving float_kill