Subroutines are usually small in size, which means they are much easier to write, test and debug. They are also easy for someone else to understand. As they are written outside of the main program, ...

In the OpenFOAM's native solver, reactingFoam, the original transport equations for species mass fractions were based on the simple asscumptions of unity Schmidt number. And the dynamic viscosity was ...