diff --git a/applications/utilities/miscellaneous/foamListTimes/Make/files b/applications/utilities/miscellaneous/foamListTimes/Make/files
deleted file mode 100644
index fe66ad822487d44cab4d6b56f3fc3013cb55aca7..0000000000000000000000000000000000000000
--- a/applications/utilities/miscellaneous/foamListTimes/Make/files
+++ /dev/null
@@ -1,3 +0,0 @@
-foamListTimes.C
-
-EXE = $(FOAM_APPBIN)/foamListTimes
diff --git a/applications/utilities/miscellaneous/foamListTimes/Make/options b/applications/utilities/miscellaneous/foamListTimes/Make/options
deleted file mode 100644
index 18e6fe47afacb902cddccf82632772447704fd88..0000000000000000000000000000000000000000
--- a/applications/utilities/miscellaneous/foamListTimes/Make/options
+++ /dev/null
@@ -1,2 +0,0 @@
-/* EXE_INC = */
-/* EXE_LIBS = */
diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C
index 61335b315dc1ffb4c9360088ca289e2cb0a591c6..3bd6b97e935b4e6ae1ffd43e20753b8d930e8d8f 100644
--- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C
@@ -285,7 +285,7 @@ bool Foam::KinematicParcel<ParcelType>::move
     const polyMesh& mesh = cloud.pMesh();
     const polyBoundaryMesh& pbMesh = mesh.boundaryMesh();
     const cloudSolution& solution = cloud.solution();
-    const scalarField& cellLengthScale = td.cloud().cellLengthScale();
+    const scalarField& cellLengthScale = cloud.cellLengthScale();
     const scalar maxCo = solution.maxCo();
 
     while (ttd.keepParticle && !ttd.switchProcessor && p.stepFraction() < 1)
@@ -457,7 +457,7 @@ void Foam::KinematicParcel<ParcelType>::hitPatch
 {
     td.keepParticle = false;
 
-    td.cloud().patchInteraction().addToEscapedParcels(nParticle_*mass());
+    cloud.patchInteraction().addToEscapedParcels(nParticle_*mass());
 }
 
 
diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
index 5cf4614cc903b0834f38ae22fc2c25563ecda476..6dcb0aaf1bfa0b4438efe05d1d625857b01ae8b7 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
@@ -58,10 +58,8 @@ void Foam::ReactingParcel<ParcelType>::calcPhaseChange
     scalarField& Cs
 )
 {
-    typedef typename TrackCloudType::reactingCloudType reactingCloudType;
-    const CompositionModel<reactingCloudType>& composition =
-        cloud.composition();
-    PhaseChangeModel<reactingCloudType>& phaseChange = cloud.phaseChange();
+    const auto& composition = cloud.composition();
+    auto& phaseChange = cloud.phaseChange();
 
     if (!phaseChange.active() || (YPhase < SMALL))
     {
@@ -396,9 +394,7 @@ void Foam::ReactingParcel<ParcelType>::calc
     const label celli
 )
 {
-    typedef typename TrackCloudType::reactingCloudType reactingCloudType;
-    const CompositionModel<reactingCloudType>& composition =
-        cloud.composition();
+    const auto& composition = cloud.composition();
 
 
     // Define local properties at beginning of time step
diff --git a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
index ce9bb1b2b7a37b22aa2db45ddea99e126ee5b76d..fa2d780142722414fdc7cc016bda3825f1a57c4a 100644
--- a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
+++ b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
@@ -68,9 +68,8 @@ void Foam::SprayParcel<ParcelType>::calc
     const label celli
 )
 {
-    typedef typename TrackCloudType::reactingCloudType reactingCloudType;
-    const CompositionModel<reactingCloudType>& composition =
-        cloud.composition();
+    const auto& composition = cloud.composition();
+    const auto& liquids = composition.liquids();
 
     // Check if parcel belongs to liquid core
     if (liquidCore() > 0.5)
@@ -80,28 +79,28 @@ void Foam::SprayParcel<ParcelType>::calc
     }
 
     // Get old mixture composition
-    scalarField X0(composition.liquids().X(this->Y()));
+    scalarField X0(liquids.X(this->Y()));
 
     // Check if we have critical or boiling conditions
-    scalar TMax = composition.liquids().Tc(X0);
+    scalar TMax = liquids.Tc(X0);
     const scalar T0 = this->T();
     const scalar pc0 = this->pc_;
-    if (composition.liquids().pv(pc0, T0, X0) >= pc0*0.999)
+    if (liquids.pv(pc0, T0, X0) >= pc0*0.999)
     {
         // Set TMax to boiling temperature
-        TMax = composition.liquids().pvInvert(pc0, X0);
+        TMax = liquids.pvInvert(pc0, X0);
     }
 
     // Set the maximum temperature limit
     cloud.constProps().setTMax(TMax);
 
     // Store the parcel properties
-    this->Cp() = composition.liquids().Cp(pc0, T0, X0);
-    sigma_ = composition.liquids().sigma(pc0, T0, X0);
-    const scalar rho0 = composition.liquids().rho(pc0, T0, X0);
+    this->Cp() = liquids.Cp(pc0, T0, X0);
+    sigma_ = liquids.sigma(pc0, T0, X0);
+    const scalar rho0 = liquids.rho(pc0, T0, X0);
     this->rho() = rho0;
     const scalar mass0 = this->mass();
-    mu_ = composition.liquids().mu(pc0, T0, X0);
+    mu_ = liquids.mu(pc0, T0, X0);
 
     ParcelType::calc(cloud,td, dt, celli);
 
@@ -114,16 +113,16 @@ void Foam::SprayParcel<ParcelType>::calc
         // Update Cp, sigma, density and diameter due to change in temperature
         // and/or composition
         scalar T1 = this->T();
-        scalarField X1(composition.liquids().X(this->Y()));
+        scalarField X1(liquids.X(this->Y()));
 
-        this->Cp() = composition.liquids().Cp(this->pc_, T1, X1);
+        this->Cp() = liquids.Cp(this->pc_, T1, X1);
 
-        sigma_ = composition.liquids().sigma(this->pc_, T1, X1);
+        sigma_ = liquids.sigma(this->pc_, T1, X1);
 
-        scalar rho1 = composition.liquids().rho(this->pc_, T1, X1);
+        scalar rho1 = liquids.rho(this->pc_, T1, X1);
         this->rho() = rho1;
 
-        mu_ = composition.liquids().mu(this->pc_, T1, X1);
+        mu_ = liquids.mu(this->pc_, T1, X1);
 
         scalar d1 = this->d()*cbrt(rho0/rho1);
         this->d() = d1;
@@ -158,29 +157,15 @@ void Foam::SprayParcel<ParcelType>::calcAtomization
     const label celli
 )
 {
-<<<<<<< HEAD
-    typedef typename TrackData::cloudType::sprayCloudType sprayCloudType;
-    const AtomizationModel<sprayCloudType>& atomization =
-        td.cloud().atomization();
+    const auto& atomization = cloud.atomization();
 
     if (!atomization.active())
     {
         return;
     }
 
-    typedef typename TrackData::cloudType::reactingCloudType reactingCloudType;
-=======
-    typedef typename TrackCloudType::reactingCloudType reactingCloudType;
->>>>>>> 4fd4fadab... lagrangian: Un-templated the tracking data
-    const CompositionModel<reactingCloudType>& composition =
-        cloud.composition();
-
-<<<<<<< HEAD
-=======
-    typedef typename TrackCloudType::sprayCloudType sprayCloudType;
-    const AtomizationModel<sprayCloudType>& atomization =
-        cloud.atomization();
->>>>>>> 4fd4fadab... lagrangian: Un-templated the tracking data
+    const auto& composition = cloud.composition();
+    const auto& liquids = composition.liquids();
 
     // Average molecular weight of carrier mix - assumes perfect gas
     scalar Wc = this->rhoc_*RR*this->Tc()/this->pc();
@@ -208,7 +193,7 @@ void Foam::SprayParcel<ParcelType>::calcAtomization
     scalar chi = 0.0;
     if (atomization.calcChi())
     {
-        chi = this->chi(cloud, td, composition.liquids().X(this->Y()));
+        chi = this->chi(cloud, td, liquids.X(this->Y()));
     }
 
     atomization.update
@@ -353,24 +338,23 @@ Foam::scalar Foam::SprayParcel<ParcelType>::chi
 {
     // Modifications to take account of the flash boiling on primary break-up
 
-    typedef typename TrackCloudType::reactingCloudType reactingCloudType;
-    const CompositionModel<reactingCloudType>& composition =
-        cloud.composition();
+    const auto& composition = cloud.composition();
+    const auto& liquids = composition.liquids();
 
     scalar chi = 0.0;
     scalar T0 = this->T();
     scalar p0 = this->pc();
     scalar pAmb = cloud.pAmbient();
 
-    scalar pv = composition.liquids().pv(p0, T0, X);
+    scalar pv = liquids.pv(p0, T0, X);
 
-    forAll(composition.liquids(), i)
+    forAll(liquids, i)
     {
         if (pv >= 0.999*pAmb)
         {
             // Liquid is boiling - calc boiling temperature
 
-            const liquidProperties& liq = composition.liquids().properties()[i];
+            const liquidProperties& liq = liquids.properties()[i];
             scalar TBoil = liq.pvInvert(p0);
 
             scalar hl = liq.hl(pAmb, TBoil);