diff --git a/applications/solvers/combustion/coalChemistryFoam/coalChemistryFoam.C b/applications/solvers/combustion/coalChemistryFoam/coalChemistryFoam.C
index 4b28e3b37f5f203f6bb4d0c7d989f0ca96e07500..142eb51a2f9d3dab8b73859321937f78f81a1e18 100644
--- a/applications/solvers/combustion/coalChemistryFoam/coalChemistryFoam.C
+++ b/applications/solvers/combustion/coalChemistryFoam/coalChemistryFoam.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 1991-2007 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -76,6 +76,9 @@ int main(int argc, char *argv[])
 
         Info<< "Time = " << runTime.timeName() << nl << endl;
 
+        rhoEffLagrangian = coalParcels.rhoEff() + limestoneParcels.rhoEff();
+        pDyn = 0.5*rho*magSqr(U);
+
         Info << "Evolving coal cloud" << endl;
 
         coalParcels.evolve();
diff --git a/applications/solvers/combustion/coalChemistryFoam/createFields.H b/applications/solvers/combustion/coalChemistryFoam/createFields.H
index 0cc8632562d06d78195add6f38482b7caa09101d..9c10e89164c0919b7ad3d4d07f7128688da78903 100644
--- a/applications/solvers/combustion/coalChemistryFoam/createFields.H
+++ b/applications/solvers/combustion/coalChemistryFoam/createFields.H
@@ -28,50 +28,36 @@
         thermo->rho()
     );
 
-// lagrangian coal density field
-/*    volScalarField rholc
+    // lagrangian effective density field - used externally (optional)
+    volScalarField rhoEffLagrangian
     (
         IOobject
         (
-            "rholc",
+            "rhoEffLagrangian",
             runTime.timeName(),
             mesh,
             IOobject::NO_READ,
             IOobject::AUTO_WRITE
         ),
         mesh,
-        dimensionedScalar("zero", dimensionSet(1, -3, 0, 0, 0, 0, 0), 0.0)
+        dimensionedScalar("zero", dimDensity, 0.0)
     );
 
-// lagrangian limestone density field
-    volScalarField rhols
+    // dynamic pressure field - used externally (optional)
+    volScalarField pDyn
     (
         IOobject
         (
-            "rhols",
+            "injectionIndicator",
             runTime.timeName(),
             mesh,
             IOobject::NO_READ,
             IOobject::AUTO_WRITE
         ),
         mesh,
-        dimensionedScalar("zero", dimensionSet(1, -3, 0, 0, 0, 0, 0), 0.0)
+        dimensionedScalar("zero", dimPressure, 0.0)
     );
 
-// lagrangian total density field
-    volScalarField rhol
-    (
-        IOobject
-        (
-            "rhol",
-            runTime.timeName(),
-            mesh,
-            IOobject::NO_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh,
-        dimensionedScalar("zero", dimensionSet(1, -3, 0, 0, 0, 0, 0), 0.0)
-    );*/
 
     Info<< "\nReading field U\n" << endl;
     volVectorField U