[adib_mohammad]

I need to determine the number of carbon and hydrogen atoms of certain species in reactingFoam. I am using hePsiThermo as thermoType.

I looked at the documentation, and I noticed that it can be accessed through chemistry.specieComp() used in DAC.C here.

However, when I use it in my code as

     const List<List<specieElement>>& specieComposition =
         chemistry.specieComp();

compiler throw an erorr saying that

./sootTest.H:44:20: error: ‘class Foam::BasicChemistryModel<Foam::psiReactionThermo>’ has no member named ‘specieComp’
   44 |          chemistry.specieComp();

This is the definition of the chemistry in my solver.

autoPtr<psiReactionThermo> pThermo(psiReactionThermo::New(mesh));

psiReactionThermo& thermo = pThermo();
thermo.validate(args.executable(), "h", "e");

autoPtr<BasicChemistryModel<psiReactionThermo>> pChemistry
(
    BasicChemistryModel<psiReactionThermo>::New(thermo)
);

BasicChemistryModel<psiReactionThermo>& chemistry = pChemistry();