How to compute consumption/production rate of species
 

Dears,

How can I compute consumption or production rate of a specific species using UDF? I guess I should use the DEFINE_VR_RATE and I know the algorithm of calculation of O2 consumption rate for instance:
forAll Reactions
{
O2_consumption +=reaction->reaction rate * stoichiometric coefficient;
}