ttiny_nlcg.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
       ---
       ttiny_nlcg.sh (1186B)
       ---
            1 #!/bin/bash
            2 # Tests a simple inversion setup.
            3 # Requires PISM's Python bindings and siple.
            4 PYTHONEXEC=$5
            5 PISM_BUILD_DIR=$1
            6 
            7 # make sure that Python imports the right modules
            8 export PYTHONPATH=${PISM_BUILD_DIR}/site-packages:$PYTHONPATH
            9 
           10 # check if siple is installed
           11 $PYTHONEXEC -c 'import siple'
           12 if [ $? != 0 ];
           13 then
           14     echo "Please make sure that siple is installed!"
           15     exit 1
           16 fi
           17 
           18 set -x
           19 set -e
           20 
           21 # Create input files
           22 $PYTHONEXEC build_tiny.py -Mx 9 -My 9
           23 
           24 $PYTHONEXEC make_synth_ssa.py -i tiny.nc -o inv_data.nc \
           25               -pseudo_plastic -pseudo_plastic_q 0.25 -regional \
           26               -ssa_dirichlet_bc -generate_ssa_observed -ssa_method fem \
           27               -design_prior_const 70000 -inv_ssa tauc
           28 
           29 # Run the inversion code
           30 $PYTHONEXEC pismi.py \
           31               -i tiny.nc -pseudo_plastic -pseudo_plastic_q 0.25 -inv_data inv_data.nc \
           32               -o tiny_inv.nc -regional -ssa_dirichlet_bc -inv_use_tauc_prior \
           33               -inv_design_param trunc -inv_design_cL2 1 -inv_design_cH1 0 \
           34               -inv_method nlcg -inv_target_misfit 100 -report_coverage
           35 
           36 # Check if we succeeded
           37 $PYTHONEXEC verify_ssa_inv.py tiny_inv.nc --desired_misfit 100 --iter_max 40 --morozov