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). The ZRDAAM is the radiation damage model for zircon published in Guenthner et al. (2013).

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.

0 = apatite or 1 = zircon
U ppm
Th ppm
Sm ppm
grain radius in microns

number of time-Temperature pairs, i.e., 2
list t-T points in descending order

i.e.,
2
500 150
0 0

./euage <#of eU values>. <step between eU values>