tiny_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