ENH: Add functionality to singleStepReactingMixture

src/thermophysicalModels/reactionThermo/mixtures/singleStepReactingMixture/singleStepReactingMixture.C

@@ -108,6 +108,8 @@ void Foam::singleStepReactingMixture<ThermoType>::calculateMaxProducts()
         Yprod0_[specieI] =  this->speciesData()[specieI].W()/Wm*Xi[i];
     }
 
+    Info << "Maximum products mass concentrations :" << Yprod0_<< endl;
+
     // Normalize the stoichiometric coeff to mass
     forAll(specieStoichCoeffs_, i)
     {

src/thermophysicalModels/reactionThermo/mixtures/singleStepReactingMixture/singleStepReactingMixture.H

@@ -152,6 +152,9 @@ public:
             //- Return the list to indicate if specie is produced/consumed
             inline const List<int>& specieProd() const;
 
+            //- Return the list of products mass concentrations
+            inline const scalarList& Yprod0() const;
+
 
         // I-O
 

src/thermophysicalModels/reactionThermo/mixtures/singleStepReactingMixture/singleStepReactingMixtureI.H

@@ -94,4 +94,12 @@ Foam::singleStepReactingMixture<ThermoType>::specieProd() const
 }
 
 
+template<class ThermoType>
+inline const Foam::scalarList&
+Foam::singleStepReactingMixture<ThermoType>::Yprod0() const
+{
+    return Yprod0_;
+}
+
+
 // ************************************************************************* //