Output Files - Geronimorz/Test GitHub Wiki

Output Directory Structure

After running an inversion, QtomoMCMC generates an output directory named after the input YAML file (e.g., TongaLau_Spherical/). Below is an overview of its structure and contents:


Top-Level Contents

Item Description
inp.yml A copy of the input configuration file used for this run.
data/ Preprocessed t*, raypath, and trace data used internally during inversion.
figures/ All model visualization outputs, separated into subfolders by type.
models/ Raw MCMC chain outputs saved in .jld format.
TongaAttenData/ (Optional) Directory for storing interpolated models for plotting in GMT.

data/

Only needs to run prepare_data.jl once to pre-process the data.

Contains internal representations of the inputs used for inversion:

  • p_tstar.dat, raypaths.dat – input files reformatted for internal use
  • raypaths.jld, traces.jld – binary storage of ray and data arrays in Julia's JLD format

figures/

⚠️ Only used to visualize figures directly in Julia for quick inspection and diagnostics. For publication-quality figures, it is recommended to export the results and use GMT (Generic Mapping Tools).

Folder Description
llh/ Likelihood evolution per iteration across chains
nCells/ Number of Voronoi cells per model over iterations
phi/ Model misfit evolution per iteration across chains
xyMean/, xzMean/ Mean model in horizontal (xy) and vertical (xz) cross sections
xyUncertainty/, xzUncertainty/ Standard deviation (σ) maps representing model uncertainty
xyMasked/, xzMasked/ Mean models masked by standard deviation
xyContour/, xzContour/ Contour plots of attenuation strength
xyVoronoi_*/ Individual Voronoi realizations plotted at selected depths (e.g., 50, 100, 300, 500 km)
xzVoronoi_-20/ Voronoi realizations at a specific y-slice (e.g., y = -20 km)

models/

MCMC model chain outputs saved as .jld files.

e.g.,
chain1_iter1000_100%.jld = 1000 iterations of chain 1, all models saved.