Gear metal diffusion
Compare changes
@@ -52,21 +57,24 @@ ids = np.linspace(1, numPart, numPart)
@@ -100,7 +108,8 @@ grp.create_dataset('Masses', data=m, dtype='f')
In this merge request, I am implementing the metal diffusion from Shen et al. 2010. We still need to test it in a cosmological context, conservation of mass and find the best coefficient / time step (e.g. work on the ideas from Greif et al. 2009), but this will be for a future PhD student.