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