BUG: Corrected cloud radiation model handling

src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -477,7 +477,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::calc
             const scalar T4 = pow4(this->T_);
             td.cloud().radAreaP()[cellI] += dt*np0*ap;
             td.cloud().radT4()[cellI] += dt*np0*T4;
-            td.cloud().radAreaP()[cellI] += dt*np0*ap*T4;
+            td.cloud().radAreaPT4()[cellI] += dt*np0*ap*T4;
         }
     }
 }

src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C

@@ -466,7 +466,7 @@ void Foam::ReactingParcel<ParcelType>::calc
             const scalar T4 = pow4(this->T_);
             td.cloud().radAreaP()[cellI] += dt*np0*ap;
             td.cloud().radT4()[cellI] += dt*np0*T4;
-            td.cloud().radAreaP()[cellI] += dt*np0*ap*T4;
+            td.cloud().radAreaPT4()[cellI] += dt*np0*ap*T4;
         }
     }
 }

src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -260,7 +260,7 @@ void Foam::ThermoParcel<ParcelType>::calc
             const scalar T4 = pow4(this->T_);
             td.cloud().radAreaP()[cellI] += dt*np0*ap;
             td.cloud().radT4()[cellI] += dt*np0*T4;
-            td.cloud().radAreaP()[cellI] += dt*np0*ap*T4;
+            td.cloud().radAreaPT4()[cellI] += dt*np0*ap*T4;
         }
     }
 }