Python Scripts 2, some examples to visualize the output files of TomoATT

TomoATT_path/examples/scripts_of_plotting provides some examples to visualize the output files of TomoATT

Users can directly run this example simply by


./run_this_example.sh

Parameter settings about output files can be found in the YAML file as follows


############################################
#            output file setting           #
############################################
output_setting:
  output_dir: OUTPUT_FILES/OUTPUT_FILES_inv # path to output director (default is ./OUTPUT_FILES/)
  output_source_field:     false # True: output the traveltime field and adjoint field of all sources at each iteration. Default: false. File: 'out_data_sim_group_X'.
  output_kernel:           false # True: output sensitivity kernel and kernel density. Default: false. File: 'out_data_sim_group_X'.
  output_final_model:      true # True: output merged final model. This file can be used as the input model for TomoATT. Default: true. File: 'model_final.h5'.
  output_middle_model:     false # True: output merged intermediate models during inversion. This file can be used as the input model for TomoATT. Default: false. File: 'middle_model_step_XXXX.h5'
  output_in_process:       false # True: output at each inv iteration, otherwise, only output step 0, Niter-1, Niter. Default: true. File: 'out_data_sim_group_0'.
  output_in_process_data:  false # True: output src_rec_file at each inv iteration, otherwise, only output step 0, Niter-2, Niter-1. Default: true. File: 'src_rec_file_step_XXXX.dat'
  single_precision_output: false # True: output results in single precision. Default: false.   
  verbose_output_level:    0 # output internal parameters, (to do).
  output_file_format: 0 # 0: hdf5, 1: ascii
 
# output files:
# File: 'out_data_grid.h5'. Keys: ['Mesh']['elem_conn'], element index; 
#                                 ['Mesh']['node_coords_p'], phi coordinates of nodes; 
#                                 ['Mesh']['node_coords_t'], theta coordinates of nodes; 
#                                 ['Mesh']['node_coords_r'], r coordinates of nodes;
#                                 ['Mesh']['node_coords_x'], phi coordinates of elements; 
#                                 ['Mesh']['node_coords_y'], theta coordinates of elements; 
#                                 ['Mesh']['node_coords_z'], r coordinates of elements; 
# File: 'out_data_sim_group_0'. Keys: ['model']['vel_inv_XXXX'], velocity model at iteration XXXX; 
#                                     ['model']['xi_inv_XXXX'], xi model at iteration XXXX; 
#                                     ['model']['eta_inv_XXXX'], eta model at iteration XXXX
#                                     ['model']['Ks_inv_XXXX'], sensitivity kernel related to slowness at iteration XXXX
#                                     ['model']['Kxi_inv_XXXX'], sensitivity kernel related to xi at iteration XXXX
#                                     ['model']['Keta_inv_XXXX'], sensitivity kernel related to eta at iteration XXXX
#                                     ['model']['Ks_density_inv_XXXX'], kernel density of Ks at iteration XXXX
#                                     ['model']['Kxi_density_inv_XXXX'], kernel density of Kxi at iteration XXXX
#                                     ['model']['Keta_density_inv_XXXX'], kernel density of Keta at iteration XXXX
#                                     ['model']['Ks_over_Kden_inv_XXXX'], slowness kernel over kernel density at iteration XXXX
#                                     ['model']['Kxi_over_Kden_inv_XXXX'], xi kernel over kernel density at iteration XXXX
#                                     ['model']['Keta_over_Kden_inv_XXXX'], eta kernel over kernel density at iteration XXXX
#                                     ['model']['Ks_update_inv_XXXX'], slowness kernel over kernel density at iteration XXXX, smoothed by inversion grid
#                                     ['model']['Kxi_update_inv_XXXX'], xi kernel over kernel density at iteration XXXX, smoothed by inversion grid
#                                     ['model']['Keta_update_inv_XXXX'], eta kernel over kernel density at iteration XXXX, smoothed by inversion grid
#                                     ['1dinv']['vel_1dinv_inv_XXXX'], 2d velocity model at iteration XXXX, in 1d inversion mode
#                                     ['1dinv']['r_1dinv'], r coordinates (depth), in 1d inversion mode
#                                     ['1dinv']['t_1dinv'], t coordinates (epicenter distance), in 1d inversion mode
# File: 'src_rec_file_step_XXXX.dat' or 'src_rec_file_forward.dat'. The synthetic traveltime data file.
# File: 'final_model.h5'. Keys: ['eta'], ['xi'], ['vel'], the final model.
# File: 'middle_model_step_XXXX.h5'. Keys: ['eta'], ['xi'], ['vel'], the model at step XXXX.
# File: 'inversion_grid.txt'. The location of inversion grid nodes
# File: 'objective_function.txt'. The objective function value at each iteration
# File: 'out_data_sim_group_X'. Keys: ['src_YYYY']['time_field_inv_XXXX'], traveltime field of source YYYY at iteration XXXX;
#                                     ['src_YYYY']['adjoint_field_inv_XXXX'], adjoint field of source YYYY at iteration XXXX;
#                                     ['1dinv']['time_field_1dinv_YYYY_inv_XXXX'], 2d traveltime field of source YYYY at iteration XXXX, in 1d inversion mode
#                                     ['1dinv']['adjoint_field_1dinv_YYYY_inv_XXXX'], 2d adjoint field of source YYYY at iteration XXXX, in 1d inversion mode

Here we will demonstrate how to run this example step-by-step.