A solution of a one-dimensional diffusion equation describing He concentration in diamonds (sphere), affected by production (decay of U and Th) and diffusion. The solution is by finite difference - Crank-Nicolson method. This model was generated for Weiss et al., 2021, Helium in diamonds unravels over a billion years of craton metasomatism, Nature Communications Please cite Weiss et al., 2021 if you use this model.
- An example of an input dataset to demo the software code can be found in supplementary data table 1 at Weiss et al., 2021.
- The code runs in MATLAB
- Instructions to run the code and the expected output are documented in the code file.