Perform seismic tomography

The entire inversion process consists of three steps:

perform 1-D model inversion to obtain a better initial model for subsequent relocation and model update;
relocate earthquake in the obtained initial model;
update model parameters and relocate earthquakes simultaneously;

Step 1. 1-D model inversion

Given the three input files

users can perform inversion by


mpirun -n 8 TOMOATT -i 3_input_params/input_params_step1_1D_inv.yaml

The results are output to the folder OUTPUT_FILES/OUTPUT_FILES_step1_1D_inv.

The updated model can be found as OUTPUT_FILES/OUTPUT_FILES_step1_1D_inv/final_model.h5.

Given the three input files

Parameter file in YAML format 3_input_params/input_params_step2_reloc.yaml.
Initial model file in HDF5 format OUTPUT_FILES/OUTPUT_FILES_step1_1D_inv/final_model.h5.
Traveltime data file in TEXT format 1_src_rec_files/src_rec_file.dat.

users can perform inversion by


mpirun -n 8 TOMOATT -i 3_input_params/input_params_step2_reloc.yaml

The results are output to the folder OUTPUT_FILES/OUTPUT_FILES_step2_reloc.

The updated catalog can be found as OUTPUT_FILES/OUTPUT_FILES_step2_reloc/src_rec_file_reloc_0030_obs.dat.

Given the three input files

Parameter file in YAML format 3_input_params/OUTPUT_FILES_step3_inv_reloc.yaml.
Initial model file in HDF5 format OUTPUT_FILES/OUTPUT_FILES_step1_1D_inv/final_model.h5.
Traveltime data file in TEXT format OUTPUT_FILES/OUTPUT_FILES_step2_reloc/src_rec_file_reloc_0030_obs.dat.

users can perform inversion by


mpirun -n 8 TOMOATT -i 3_input_params/input_params_step3_inv_reloc.yaml

The results are output to the folder OUTPUT_FILES/OUTPUT_FILES_step3_inv_reloc.

The updated model can be found as OUTPUT_FILES/OUTPUT_FILES_step3_inv_reloc/final_model.h5.