geoprior.scripts.plot_sm3_identifiability#

Plot synthetic identifiability (SM3).

Panels#

  1. Timescale recovery: log10(tau_true) vs log10(tau_est) + optional tau_cl overlay (closure-implied) + optional tau_prior marker.

  2. Permeability recovery: log10(K_true) vs log10(K_est) OR log10(K_true) vs log10(K_cl(tau_est)).

  3. Degeneracy ridge check: delta_K vs (delta_Ss + 2 delta_Hd).

  4. Error vs identifiability metric: \(|\Delta log10 tau|\) vs chosen identifiability metric.

Inputs#

CSV file produced by SM3 synth runs, typically:

results/sm3_synth_1d/sm3_synth_runs.csv

Outputs#

  • Figure: scripts/figs/<out>.png and <out>.svg

  • Summary table: scripts/out/<out_csv>

  • Summary json: scripts/out/<out_json>

API conventions#

  • Style via scripts.utils.set_paper_style()

  • Figure output via scripts.utils.resolve_fig_out()

  • Tables/JSON via scripts.utils.resolve_out_out()

  • Common plot text toggles via scripts.utils.add_plot_text_args()

Notes

  • Bootstrap CIs are computed row-wise.

  • If –k-from-tau=auto:
    • defaults to false when identify includes “both”

    • else defaults to true.

Functions

main([argv, prog])

plot_sm3_identifiability(df, *, tau_units, ...)

plot_sm3_identifiability_main([argv, prog])

geoprior.scripts.plot_sm3_identifiability.plot_sm3_identifiability(df, *, tau_units, metric, k_from_tau, k_cl_source, show_prior, show_tau_cl, show_legend, show_labels, show_ticklabels, show_panel_titles, show_title, title, show_stats, show_ci, n_boot, ci, boot_seed, out, out_csv, out_json, dpi)[source]#
Parameters:
Return type:

None

geoprior.scripts.plot_sm3_identifiability.plot_sm3_identifiability_main(argv=None, *, prog=None)[source]#
Parameters:
Return type:

None

geoprior.scripts.plot_sm3_identifiability.main(argv=None, *, prog=None)[source]#
Parameters:
Return type:

None