Possible unit error in flux production profile for EmittingEarthAtm?

[pweigel]

The production rate here comes from the interpolator which was created using d(flux)/d(height), where the height is in units of cm:

nuSQuIDS/src/body.cpp

Line 1660 in e9919a2

double prod_rate = flux_interpolators[flv][rho](cz, height_cm, energy_GeV);

This differential flux is added directly to current_external_flux,

nuSQuIDS/src/nuSQuIDS.cpp

Line 288 in e9919a2

body->injected_neutrino_flux(current_external_flux,*track,*this);

Since this is a source term, I think it's being integrated over the track length in natural units. If this is the case, then the production rate should be divided by param.cm.