fixed MachS calculation bug: MachS = sqrt(density weighted mean vrms^2)/cs_iso, I was incorrectly doing: MachS= density weighted mean vrms/cs_iso