Some cleaning up.

applications/solvers/cfdemSolverRecurrence/cfdemSolverRecurrence.C

@@ -39,10 +39,11 @@ Rules
 #include "singlePhaseTransportModel.H"
 #include "turbulenceModel.H"
 #include "cfdemCloudRec.H"
+#include "cfdemCloud.H"
 #include "recBase.H"
 #include "recModel.H"
 #include "clockModel.H"
-#include "cfdemCloud.H"
+
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 int main(int argc, char *argv[])

applications/solvers/tfmSolverRecurrence/Make/options

@@ -10,9 +10,12 @@ EXE_INC = \
     -I$(LIB_SRC)/sampling/lnInclude \
     -I$(CFDEM_SRC_DIR)/lagrangian/cfdemParticle/lnInclude \
     -I$(CFDEM_SRC_DIR)/lagrangian/cfdemParticle/cfdTools \
+    -I$(CFDEM_SRC_DIR)/recurrence/lnInclude \
+    -I$(CFDEM_SRC_DIR)/lagrangian/cfdemParticle/derived/cfdemCloudRec
 
 EXE_LIBS = \
     -L$(CFDEM_LIB_DIR)\
+    -lrecurrence \
     -lincompressibleRASModels \
     -lincompressibleLESModels \
     -lincompressibleTransportModels \

applications/solvers/tfmSolverRecurrence/readFields.H

@@ -1,5 +1,5 @@
-alpha1Rec=*(particleCloud.recM().voidfraction());
+alpha1Rec=*(recurrenceBase.recM().voidfraction());
 alpha2Rec=1-alpha1Rec;
-U1Rec=*(particleCloud.recM().U());
-U2Rec=*(particleCloud.recM().Us());
-phi2Rec=*(particleCloud.recM().phiS());
+U1Rec=*(recurrenceBase.recM().U());
+U2Rec=*(recurrenceBase.recM().Us());
+phi2Rec=*(recurrenceBase.recM().phiS());

applications/solvers/tfmSolverRecurrence/tfmSolverRecurrence.C

@@ -40,6 +40,9 @@ Rules
 #include "singlePhaseTransportModel.H"
 #include "turbulenceModel.H"
 #include "cfdemCloudRec.H"
+#include "cfdemCloud.H"
+#include "recBase.H"
+#include "recModel.H"
 #include "clockModel.H"
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -50,14 +53,15 @@ int main(int argc, char *argv[])
     #include "createMesh.H"
     #include "createFields.H"
     #include "createFvOptions.H"
-    cfdemCloudRec particleCloud(mesh);
+    cfdemCloudRec<cfdemCloud> particleCloud(mesh);
+    recBase recurrenceBase(mesh);
 
     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
     Info<< "\nCalculating particle trajectories based on recurrence statistics\n" << endl;
     
     label recTimeIndex(0);
-    scalar recTimeStep_=particleCloud.recM().recTimeStep();
+    scalar recTimeStep_=recurrenceBase.recM().recTimeStep();
     
     while (runTime.run())
     {
@@ -77,7 +81,7 @@ int main(int argc, char *argv[])
         {
 	    Info << "Updating fields at run time " << runTime.timeOutputValue()
 	        << " corresponding to recurrence time " << (recTimeIndex+1)*recTimeStep_ << ".\n" << endl;
-            particleCloud.updateRecFields();
+            recurrenceBase.updateRecFields();
 	    #include "readFields.H"
             recTimeIndex++;
         }

src/lagrangian/cfdemParticle/derived/cfdemCloudRec/cfdemCloudRec.H

@@ -96,7 +96,7 @@ public:
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-#include "cfdemCloudRec.TPP"
+#include "cfdemCloudRecI.H"
 
 #endif