trun_shelfonly.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_shelfonly.sh (1108B)
       ---
            1 #!/bin/bash
            2 
            3 N=4
            4 xx=151
            5 yy=151
            6 length=10
            7 
            8 infile="circular_shelfonly.nc"
            9 
           10 ./circular_dirichlet.py -Mx $xx -My $yy -o $infile -shelf
           11 
           12 grid="-Mx $xx -My $yy -Mz 31 -Mbz 1 -Lz 1500 -Lbz 1000"
           13 
           14 extra="-extra_times 1 -extra_vars thk,mask,velbar_mag,ice_area_specific_volume,velbar -extra_file "
           15 
           16 pismopts="-i $infile -bootstrap $grid -stress_balance ssa+sia -ssa_dirichlet_bc"
           17 
           18 doit="mpiexec -n $N pismr $pismopts"
           19 
           20 # run with strength extension, the old PISM method
           21 #$doit $pismopts -y $length -o so_old.nc
           22 
           23 # check that this result is similar
           24 #$doit -y $length -energy cold -o so_old_cold.nc
           25 
           26 
           27 # run with strength extension and part_grid but no CFBC
           28 # this could be a regression for -part_grid
           29 $doit -y $length -part_grid -o so_part.nc $extra so_ex_part.nc
           30 
           31 # run with CFBC but no part_grid
           32 # this could be a regression for -ssa_method fd_pik only
           33 $doit -y $length -ssa_method fd -cfbc -o so_cfbc.nc $extra so_ex_cfbc.nc
           34 
           35 # run with CFBC and part_grid
           36 # this could be a regression for -ssa_method fd_pik only
           37 $doit $pismopts -y $length -ssa_method fd -cfbc -part_grid -o so_partcfbc.nc $extra so_ex_partcfbc.nc