pythermo 0.0.0__py3-none-any.whl

pythermo/__init__.py

@@ -0,0 +1,23 @@
+"""
+Classes and methods for performing various modeling and computational tasks in the field of low-temperature thermochronology.
+
+This package includes the following modules:
+
+    - .crystal      : Classes and methods for determining radiation damage levels, diffusivities, and thermochronometric dates from different mineral systems
+
+    - .tT_path      : Methods for creating discretized time-temperature (tT) paths and performing crystal-independent radiation damage annealing.
+
+    - .tT_model     : Methods for modeling and plotting thermochronometric date results for comparison to measured data sets. A key focus here is thermal history modeling.
+
+    - .constants    : Important universal constants used by the various modules.
+
+"""
+
+__version__ = '0.1.0'
+
+__all__ = ['constants','crystal','tT_path','tT_model']
+
+from .constants import *
+from .crystal import *
+from .tT_path import *
+from .tT_model import *

pythermo/constants.py

@@ -0,0 +1,46 @@
+""" 
+constants.py
+
+Universal libraries and constants used by files crystals.py, tT_path.py, and tT_model.py 
+
+"""
+import numpy as np
+
+#time constants
+sec_per_myr = 3.15569259747*10**13
+sec_per_yr = 3.15569259747*10**7
+
+#gas constant in kJ/K * mol-1
+gas_constant = 0.008314462618
+
+#jaffey decay constants, 1/s
+lambda_232 = np.log(2)/(1.405*10**10) * (1/sec_per_yr)
+lambda_235 = np.log(2)/(7.0381*10**8) * (1/sec_per_yr)
+lambda_238 = np.log(2)/(4.4683*10**9) * (1/sec_per_yr)
+
+#jaffey decay constants, 1/yr
+lambda_232_yr = np.log(2)/(1.405*10**10) 
+lambda_235_yr = np.log(2)/(7.0381*10**8) 
+lambda_238_yr = np.log(2)/(4.4683*10**9) 
+
+#Sm and FT decay constants, 1/s
+lambda_147 = 6.54*10**-12 * (1/sec_per_yr)
+lambda_f = 8.46*10**-17 * (1/sec_per_yr)
+
+#Sm and FT decay constants, 1/yr
+lambda_147_yr = 6.54*10**-12 
+lambda_f_yr = 8.46*10**-17 
+
+#Avogadro conversion constants
+atom_mol = 6.02214076*10**23
+atom_nmol = 6.02214076*10**14
+
+#ppm to atoms/g conversion factor, atomic weights from 2021 IUPAC report, uses 1/137.818 ratio for U isotopes
+U235_ppm_atom=(1/1000000)*(0.0072559/238.02891)*atom_mol
+U238_ppm_atom=(1/1000000)*(0.99274405/238.02891)*atom_mol
+Th_ppm_atom=(1/1000000)*(1/232.0377)*atom_mol
+Sm_ppm_atom = (1/1000000)*(0.015/150.36)*atom_mol
+
+#mineral densities g/cc
+ap_density = 3.19
+zirc_density = 4.60