If a SarracenDataFrame contains mixed particle types (say gas and dust), and the particle mass is specified in params rather than per particle, then calc_density() will use the gas mass for all particle types.
This is not an issue if particles have been separated in read_phantom() using separate_types='all', as this assigns the correct particle mass to the global params in each SarracenDataFrame.
If a SarracenDataFrame contains mixed particle types (say gas and dust), and the particle mass is specified in
paramsrather than per particle, thencalc_density()will use the gas mass for all particle types.This is not an issue if particles have been separated in
read_phantom()usingseparate_types='all', as this assigns the correct particle mass to the globalparamsin each SarracenDataFrame.