Takes the raw time-temperature (t-T) output file from the QTQt thermal history modelling software of Gallagher (2012) and replots the output as an image displaying t-T path density.

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