cr_get_modelling_errors module

Plot deviation between CRMod-derived apparent resistivities and analytical resistivity vs K-factors. Also save the relative deviations to a .dat file.

dev = (R_mod * K - rho0) / rho0

• R_mod - modelled resistance value (using CRMod)

• K - geometric factor over a homogeneous half-space, computed using the analytical formula

• rho0: resistivity of homogeneous half-space

Output files

OUTPUT.png: plot K vs. rel. modelling errors

OUTPUT.dat: Nx2 array, first column: analytical K factors,

second column: rel. modelling errors

Usage example

cr_get_modelling_errors.py –elem grid1/elem.dat –elec grid1/elec.dat –config grid1/config.dat -o grid1_modelling_error.png

cr_get_modelling_errors.compute_K_factors(elem_file, elec_file, config_file)

cr_get_modelling_errors.get_R_mod(options, rho0)

Compute synthetic measurements over a homogeneous half-space

cr_get_modelling_errors.handle_cmd_options()

cr_get_modelling_errors.main()

cr_get_modelling_errors.plot_and_save_deviations(rho0, rho_mod, Kfactors, filename, configs)