BUG: Corrected construction of inperolation objects in trackData

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

@@ -122,7 +122,7 @@ public:
             ReactingCloud<ParcelType>& cloud_;
 
             //- Interpolator for continuous phase pressure field
-            const interpolation<scalar>& pInterp_;
+            autoPtr<interpolation<scalar> > pInterp_;
 
 
     public:

src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcelI.H

@@ -170,7 +170,7 @@ template<class ParcelType>
 inline const Foam::interpolation<Foam::scalar>&
 Foam::ReactingParcel<ParcelType>::trackData::pInterp() const
 {
-    return pInterp_;
+    return pInterp_();
 }
 
 

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

@@ -138,7 +138,7 @@ public:
             //- Reference to the cloud containing this particle
             ThermoCloud<ParcelType>& cloud_;
 
-            //- Local copy of specific heat field
+            //- Local copy of carrier specific heat field
             //  Cp not stored on acrrier thermo, but returned as tmp<...>
             const volScalarField Cp_;
 
@@ -149,7 +149,7 @@ public:
                 autoPtr<interpolation<scalar> > TInterp_;
 
                 //- Specific heat capacity field interpolator
-                const interpolation<scalar>& CpInterp_;
+                autoPtr<interpolation<scalar> > CpInterp_;
 
 
     public:
@@ -170,6 +170,9 @@ public:
             //- Return access to the owner cloud
             inline ThermoCloud<ParcelType>& cloud();
 
+            //- Return access to the locally stored carrier Cp field
+            inline const volScalarField& Cp() const;
+
             //- Return const access to the interpolator for continuous
             //  phase temperature field
             inline const interpolation<scalar>& TInterp() const;

src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcelI.H

@@ -202,6 +202,14 @@ Foam::ThermoParcel<ParcelType>::trackData::cloud()
 }
 
 
+template<class ParcelType>
+inline const Foam::volScalarField&
+Foam::ThermoParcel<ParcelType>::trackData::Cp() const
+{
+    return Cp_;
+}
+
+
 template<class ParcelType>
 inline const Foam::interpolation<Foam::scalar>&
 Foam::ThermoParcel<ParcelType>::trackData::TInterp() const
@@ -214,7 +222,7 @@ template<class ParcelType>
 inline const Foam::interpolation<Foam::scalar>&
 Foam::ThermoParcel<ParcelType>::trackData::CpInterp() const
 {
-    return CpInterp_;
+    return CpInterp_();
 }