Numerical Model of Diffusion-Driven Fe–Mg Isotopic Fractionation

FeMg_diffusion_modeling.py: Python script modeling isotopic fractionation in Mg–Fe interdiffusion systems. The algorithm follows Dauphas et al., (2010) and Sio et al., (2013). MELTS_results.dat: Output from MELTS simulations for the magmatic evolution of CE-5 basalts. FeMg_diffusion_modeling.py reads this data to set initial diffusion temperatures (T0) for olivine grains with different compositions in CE-5 basalts. The four columns represent temperature (T), molar fraction of forsterite in olivine (Fo), mass fraction of olivine (ol_frac, %), mass fraction of melt (melt_frac, %).