src/engine/ignition/ignite.H

   1     forAll(ign.sites(), i)
   2     {
   3         const ignitionSite& ignSite = ign.sites()[i];
   4
   5         if (ignSite.igniting())
   6         {
   7             forAll(ignSite.cells(), icelli)
   8             {
   9                 label ignCell = ignSite.cells()[icelli];
  10                 Info<< "Igniting cell " << ignCell;
  11
  12                 Info<< " state :"
  13                     << ' ' << b[ignCell]
  14                     << ' ' << Xi[ignCell]
  15                     << ' ' << Su[ignCell]
  16                     << ' ' << mgb[ignCell]
  17                     << endl;
  18
  19                 bEqn.diag()[ignSite.cells()[icelli]] +=
  20                 (
  21                     ignSite.strength()*ignSite.cellVolumes()[icelli]
  22                    *rhou[ignSite.cells()[icelli]]/ignSite.duration()
  23                 )/(b[ignSite.cells()[icelli]] + 0.001);
  24             }
  25         }
  26     }