ttiny_tikhonov_lmvm.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_tikhonov_lmvm.sh (1067B)
       ---
            1 #!/bin/bash
            2 # Tests a simple inversion setup.
            3 # Requires PISM's Python bindings
            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 set -x
           11 set -e
           12 
           13 # Create input files
           14 $PYTHONEXEC build_tiny.py -Mx 9 -My 9
           15 
           16 $PYTHONEXEC make_synth_ssa.py -i tiny.nc -o inv_data.nc \
           17               -pseudo_plastic -pseudo_plastic_q 0.25 -regional \
           18               -ssa_dirichlet_bc -generate_ssa_observed -ssa_method fem \
           19               -design_prior_const 70000 -inv_ssa tauc
           20 
           21 # Run the inversion code
           22 $PYTHONEXEC pismi.py \
           23               -i tiny.nc -pseudo_plastic -pseudo_plastic_q 0.25 -inv_data inv_data.nc \
           24               -o tiny_tikhonov_lmvm.nc -regional -ssa_dirichlet_bc -inv_use_tauc_prior \
           25               -inv_design_param trunc -inv_design_cL2 1 -inv_design_cH1 0 \
           26               -inv_method tikhonov_lmvm -tikhonov_penalty 6e-2 -report_coverage
           27 
           28 # Check if we succeeded
           29 $PYTHONEXEC verify_ssa_inv.py tiny_tikhonov_lmvm.nc --desired_misfit 10 --misfit_tolerance .5 --iter_max 100