https://github.com/openthermochronology/euage
Runs a loop across a range of effective uranium (eU=U+0.238*Th) values (having the same U/Th/Sm proportions as the input sample, i.e., apatite or zircon) for an input thermal history using the RDAAM of Flowers et al. 2009 or ZRDAAM of Guenthner et al. 2013
https://github.com/openthermochronology/euage
effective-uranium eu radiation-damage thermal-history thermochronology
Last synced: 4 months ago
JSON representation
Runs a loop across a range of effective uranium (eU=U+0.238*Th) values (having the same U/Th/Sm proportions as the input sample, i.e., apatite or zircon) for an input thermal history using the RDAAM of Flowers et al. 2009 or ZRDAAM of Guenthner et al. 2013
- Host: GitHub
- URL: https://github.com/openthermochronology/euage
- Owner: OpenThermochronology
- License: gpl-3.0
- Created: 2023-09-19T20:41:36.000Z (over 1 year ago)
- Default Branch: main
- Last Pushed: 2023-12-07T01:14:42.000Z (over 1 year ago)
- Last Synced: 2023-12-07T02:50:24.158Z (over 1 year ago)
- Topics: effective-uranium, eu, radiation-damage, thermal-history, thermochronology
- Language: C++
- Homepage:
- Size: 64.5 KB
- Stars: 0
- Watchers: 1
- Forks: 0
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
- License: LICENSE
Awesome Lists containing this project
README
# eUage
Program eUage 1.0
Written in 2014 by Peter Zeitler (Lehigh University) in C++ to be more easily compatible with the Radiation Damage Accumulation & Annealing Model (RDAAM) code written by Richard Ketcham (UT Austin). This program was put together by grabbing the wrapper from the crankfunc code (i.e., Crank-Nicholson diffusion) and using that to call RDAAM/ZRDAAM. The RDAAM for the mineral apatite is published in [Flowers et al. (2009)](https://doi.org/10.1016/j.gca.2009.01.015). The ZRDAAM is the radiation damage model for zircon published in [Guenthner et al. (2013)](https://doi.org/10.2475/03.2013.01).
The wrapper pulls in one input file with required sample data, and another with the desired thermal history. Written only for apatite and zircon. It runs in a loop across a range of effective uranium (eU=U+0.238*Th) values, having the same U/Th/Sm proportions as the sample and provides output of the predicted date-eU relationship across the range of eU from the damage models. Minerals are modeled using spherical geometry and output 'raw' uncorrected ages and alpha-ejection corrected ages. Zoning can't be modeled.
## Input Files
### euage.in
0 = apatite or 1 = zircon
U ppm
Th ppm
Sm ppm
grain radius in microns
### euage_history.in
number of time-Temperature pairs, i.e., 2
list t-T points in descending order
i.e.,
2
500 150
0 0
## Usage
```
./euage <#of eU values>.
```