pythermo 0.0.0

Tools for performing computational tasks in the field of low-temperature thermochronology.

PyThermo

A set of classes and methods for performing various modeling and computational tasks in the field of low-temperature thermochronology. The current focus is on forward modeling of apatite and zircon (U-Th)/He data using various diffusion and damage annealing kinetic models. Future releases will expand upon the available kinetic models and mineral systems, and introduce additional methods, such as forward modeling of Arrhenius relationships.

The primary objective of this software is to provide an open-source, python-based toolkit for user adaptability and experimentation. The software includes routines for forward modeling and data plotting at higher levels that can be run in a simple fashion, but lower level algorithms are accessible as well. To that end, a secondary objective of this software is as a learning tool to remove some of the black box nature of thermal history modeling routines. Several methods are included (for example, a tridiagonal matrix solver) for instructional purposes, although the main program calls nominally faster scipy routines.

Organization

The source code consists of three separate classes and accompanying methods and/or subclasses. crystal.py contains the class crystal and currently two sub-classes apatite and zircon. Methods are devoted to calculating and parameterizing damage-diffusivity relationships and numerically solving the diffusion equation using a Crank-Nicolson approach. tT_path.py methods interpolate and discretize time-temperature (tT) paths from a handful of tT points, and calculate fission track annealing for apatite and zircon using the equivalent time concept. tT_model.py methods currently allow for one particular approach to forward modeling and plotting (U-Th)/He date-effective Uranium (eU) trends.

Installation

Eventually, PyThermo will be installed as a package by using pip. For now, you can download the 4 python files found in \src and save them to your working directory.

Usage

Once you've downloaded the python files, and if you're just interested in running some forward models, the quickest way to get started is to modify the Jupyter Notebook file that is included in the examples folder. The notebook contains markdown and code that explains and demonstrates forward model date-eU comparisons for the apatite and zircon (U-Th)/He system. The forward modeling method is one particular approach and you can (and should!) modify the forward modeling methods to suite your own needs.The tT_path and crystal classes contain several methods that you may want to call and/or adapt for your own needs. Please read the descriptions for each method in the source code for more details.

Citation

You can find various citation styles for this package here.