passing the time step size to the function calc_nodal_react_sys.
return the number of equilibrium reactions.
return the number of kinetic reactions.
only when there are kinetic and equilibrium reactions, we perform the solving.
adding countComp function.
fixing bug, getting the right number of chemical components.
if only kinetic reactions, the EqReactSys will not be properly initialized. Therefore, get the nubmer of chemical components from proper class.
adding the type of components in the MCP file.