Commit 7741d06f authored by Loic Hausammann's avatar Loic Hausammann

GEAR: sphenix does not store div_v at the same place than gadget

parent 30a0eca4
......@@ -298,8 +298,13 @@ __attribute__((always_inline)) INLINE static void star_formation_end_density(
struct part* restrict p, struct xpart* restrict xp,
const struct star_formation* cd, const struct cosmology* cosmo) {
#ifdef SPHENIX_SPH
/* Copy the velocity divergence */
xp->sf_data.div_v = p->viscosity.div_v;
#else
/* Copy the velocity divergence */
xp->sf_data.div_v = p->density.div_v;
#endif
}
/**
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment