applications/solvers/combustion/PDRFoam/hEqn.H

   1 {
   2     solve
   3     (
   4         betav*fvm::ddt(rho, h)
   5       + mvConvection->fvmDiv(phi, h)
   6       - fvm::laplacian(Db, h)
   7      ==
   8         betav*DpDt
   9     );
  10
  11     thermo.correct();
  12 }