GrainSizeTools: a Python script for estimating the dynamically recrystallized grain size from grain sectional areas
Abstract. Paleopiezometry and paleowattometry studies, required to validate models of lithospheric deformation, are increasingly common in structural geology. These studies require a numeric parameter to characterize and compare the dynamically recrystallized grain size of natural mylonites with those obtained in rocks deformed under controlled conditions in the laboratory. We introduce a new tool, a script named GrainSizeTools, to obtain a single numeric value representative of the dynamically recrystallized grain size from the measurement of grain sectional areas (2-D data). For this, it is used an estimate of the most likely grain size of the grain size population, using an alternative tool to the classical histograms and bar plots: the peak of the Gaussian kernel density estimation. The results are comparable to those that can be obtained by other stereological software available, such as the StripStar and CSDCorrections, but with the advantage that the script is specifically developed to produce a single and reproducible value avoiding manual steps in the estimation, which penalizes reproducibility.