DEFINE_NET_REACTION_RATE of species reaction rate
 

Hi everyone,
I am learning how to modify the reaction rate through DEFINE_NET_REACTION_RATE while using the EDC model. I know it is used to "compute the homogeneous net molar reaction rates of all species" and I could program the rate of each reaction, but I have no idea how to program the rate of every species.

For example,
a1A+b1B->c1C;//reaction rate is r1
c2C+d1D->e1E;//reaction rate is r2
Obviously, the rate of C = c1*r1-c2*r2, but how to describe that in the UDF

The codes below are from the header file "material.h", do they help?
real stoich_product[MAX_PDF_SPECIES]; /* stoichiometric coefficients of products */
real stoich_reactant[MAX_PDF_SPECIES]; /* stoichiometric coefficients of reactants */

What's more, after I hook that UDF, should I use a new UDF to modify the source term?

Any suggestion would be appreciated. Thank you.