Perform inversion

Given the three input files

Parameter file in YAML format 3_input_params/input_params_inv.yaml.
Initial model file in HDF5 format 2_models/model_init_N61_61_61.h5.
Traveltime data file in TEXT format OUTPUT_FILES/OUTPUT_FILES_signal/src_rec_file_forward_noisy.dat.

users can perform inversion by


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

The results are output to the folder OUTPUT_FILES/OUTPUT_FILES_inv.

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

Some remarks about model update

1. Using different types of data

In the default parameter file 3_input_params/input_params_inv.yaml, only absolute traveltimes involoved to update model parameters by setting


model_update:
 
  ...
 
  # -------------- using absolute traveltime data --------------
  abs_time:
    use_abs_time: true # 'true' for using absolute traveltime data to update model parameters; 'false' for not using (no need to set parameters in this section)
    residual_weight: [1, 3, 1, 1] # XXX is the absolute traveltime residual (second) = abs(t^{obs}_{n,i} - t^{syn}_{n,j})
    distance_weight: [100, 200, 1, 1] # XXX is epicenter distance (km) between the source and receiver related to the data
 
  # -------------- using common source differential traveltime data --------------
  cs_dif_time:
    use_cs_time: false # 'true' for using common source differential traveltime data to update model parameters; 'false' for not using (no need to set parameters in this section)
    residual_weight: [1, 3, 1, 1] # XXX is the common source differential traveltime residual (second) = abs(t^{obs}_{n,i} - t^{obs}_{n,j} - t^{syn}_{n,i} + t^{syn}_{n,j}).
    azimuthal_weight: [15, 30, 1, 1] # XXX is the azimuth difference between two separate stations related to the common source.
 
  # -------------- using common receiver differential traveltime data --------------
  cr_dif_time:
    use_cr_time: false # 'true' for using common receiver differential traveltime data to update model parameters; 'false' for not using (no need to set parameters in this section)
    residual_weight: [1, 3, 1, 1] # XXX is the common receiver differential traveltime residual (second) = abs(t^{obs}_{n,i} - t^{obs}_{m,i} - t^{syn}_{n,i} + t^{syn}_{m,i})
    azimuthal_weight: [15, 30, 1, 1] # XXX is the azimuth difference between two separate sources related to the common receiver.

Users can change the options use_abs_time, use_cs_time, use_cr_time to include absolute traveltimes and/or common-source differential arrival times and/or common-receiver differential arrival times.

Perform inversion

Some remarks about model update

1. Using different types of data

2. Important settings about model update