trun_withshelf.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_withshelf.sh (982B)
       ---
            1 #!/bin/bash
            2 
            3 N=4
            4 xx=151
            5 yy=151
            6 length=10
            7 
            8 infile=circular_withshelf.nc
            9 
           10 ./circular_ice_sheet.py -Mx $xx -My $yy -o $infile
           11 
           12 grid="-Mx $xx -My $yy -Mz 31 -Mbz 1 -Lz 4500 -Lbz 1000"
           13 
           14 pismopts="-y $length -i $infile -bootstrap $grid -stress_balance ssa+sia -ssa_method fd"
           15 
           16 doit="mpiexec -n $N pismr $pismopts"
           17 
           18 extra="-extra_times 1 -extra_vars thk,mask,velbar_mag,ice_area_specific_volume,velbar -extra_file"
           19 
           20 # run with strength extension and part_grid but no CFBC
           21 # this could be a regression for the option combination "-part_grid"
           22 $doit -part_grid -o ws_part.nc $extra ws_ex_part.nc
           23 
           24 # run with CFBC but no part_grid
           25 # this could be a regression for the option "-cfbc"
           26 $doit -cfbc -o ws_cfbc.nc $extra ws_ex_cfbc.nc
           27 
           28 # run with CFBC and part_grid
           29 # this could be a regression for the option combination "-part_grid -cfbc"
           30 # FIXME: with N=4 processors and xx=yy=301 I observe slight asymmetry in velbase_mag here?
           31 $doit -part_grid -cfbc -o ws_partcfbc.nc $extra ws_ex_partcfbc.nc