A clause preventing solution of the density equation has been removed from reactingParcelFoam and chtMultiRegionFoam, so that they are nore consistent with other compressible solvers. In general, the density equation is solved before the pimple loop is entered to make sure that the flux and the density derivative are consistent during the first pimple iteration.