BUG: density to mass density

Author: aaronjridley

src/calc_neutral_derived.cpp

@@ -585,7 +585,7 @@ void Neutrals::calc_chapman(Grid grid) {
       species[iSpecies].density_scgc.slice(iAlt) %
       species[iSpecies].scale_height_scgc.slice(iAlt);
 
-    species[iSpecies].rho_alt_int_scgc.slice(iAlt) = integral3d.slice(iAlt);
+    species[iSpecies].rho_alt_int_scgc.slice(iAlt) = integral3d.slice(iAlt) * species[iSpecies].mass;
 
     for (iAlt = nAlts - 2; iAlt >= 0; iAlt--) {
       // dr is used here instead of dalt, since we only want the radial integration, while
@@ -596,7 +596,7 @@ void Neutrals::calc_chapman(Grid grid) {
       species[iSpecies].rho_alt_int_scgc.slice(iAlt) =
         species[iSpecies].rho_alt_int_scgc.slice(iAlt + 1) +
         species[iSpecies].density_scgc.slice(iAlt) %
-        grid.dk_edge_m.slice(iAlt + 1);
+        grid.dk_edge_m.slice(iAlt + 1) * species[iSpecies].mass;
     }
 
     erfcy3d = (a + b * y3d) / (c + d * y3d + y3d % y3d);