pism

[fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
git clone git://src.adamsgaard.dk/pism # fast
git clone https://src.adamsgaard.dk/pism.git # slow
Log | Files | Refs | README | LICENSE Back to index

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