Go to the documentation of this file.00001 forAll (ign.sites(), i)
00002 {
00003 const ignitionSite& ignSite = ign.sites()[i];
00004
00005 if (ignSite.igniting())
00006 {
00007 forAll(ignSite.cells(), icelli)
00008 {
00009 label ignCell = ignSite.cells()[icelli];
00010 Info<< "Igniting cell " << ignCell;
00011
00012 Info<< " state :"
00013 << ' ' << b[ignCell]
00014 << ' ' << Xi[ignCell]
00015 << ' ' << Su[ignCell]
00016 << ' ' << mgb[ignCell]
00017 << endl;
00018
00019 bEqn.diag()[ignSite.cells()[icelli]] +=
00020 (
00021 ignSite.strength()*ignSite.cellVolumes()[icelli]
00022 *rhou[ignSite.cells()[icelli]]/ignSite.duration()
00023 )/(b[ignSite.cells()[icelli]] + 0.001);
00024 }
00025 }
00026 }
00027
00028