td_correct_temperature module¶
Tool to add or substract the temperature effect on/of resistivity data. The calculation is done after Hayley (2007):
rho_std_i = (m * (T_i - 25°) + 1) / (m * (T_std - 25°) + 1) * rho_i
Options:
m: choose coeffcient m
T_std: choose standard temperatur to which to correct the data
add: add the effect of temperature instead of substracting it (for modelling)
aniso: resistivity data is anisotropic (three columns)
- temp: file with temperature information in 3rd column with a headerline
(mag-format)
filename: file with resistivity information
rho: resistivity input and output files are in rho- instead of mag-format
output: output file
- td_correct_temperature.calc_correction(temp, mag, add=False, T_std=10, m=0.021)[source]¶
Function to add or substract the temperature effect to given data. The function can be called in python scripts. For application via command line in a file system use the script td_correct_temperature.py. The data is taken and given in Ohmm.
rho_std_i = (m * (T_i - 25°) + 1) / (m * (T_std - 25°) + 1) * rho_i rho_i = (m * (T_std - 25°) + 1) / (m * (T_i - 25°) + 1) * rho_std_i
Hayley (2007)
- Parameters:
temp – temperature values corresponding to the individual resistivity values
mag – resistivity values to be corrected
add – switch for adding instead of substracting the effect
T_std – standard temperature t or from which to correct (default=10°)
m – coeffcient (default=0.021)
- td_correct_temperature.main()[source]¶
Function to add or substract the temperature effect to data in a tomodir
- td_correct_temperature.read_iter()[source]¶
Return the path to the final rho*.mag file from the tomodir.
- td_correct_temperature.readin_rho(filename, rhofile=True, aniso=False)[source]¶
Read in the values of the resistivity in Ohmm. The format is variable: rho-file or mag-file.