|
stars_data[:,7] = data.stars.smoothing_lengths[mask_stars].value * siminfo.a * 1e3 #kpc |
This is not so physically relevant for the stellar mass distributions, could instead use a fixed smoothing, e.g:
baryon_maxsoft = float(data.metadata.parameters["Gravity:max_physical_baryon_softening"])*1000. #kpc
stars_data[:,7] = 2*baryon_maxsoft*np.ones(stars_mass.size)
morpholopy/morpholopy/particles.py
Line 66 in 0d46597
This is not so physically relevant for the stellar mass distributions, could instead use a fixed smoothing, e.g: