release on 2013-03-19_21-11-25

2025-12-08 06:37:44 +00:00 · 2013-03-19 21:11:25 +01:00
parent 2be5ead220
commit c803f5a8f1
78 changed files with 1819 additions and 835 deletions
--- a/81
+++ b/81
@ -0,0 +1,81 @@
+/*---------------------------------------------------------------------------*\
+    CFDEMcoupling - Open Source CFD-DEM coupling
+
+    CFDEMcoupling is part of the CFDEMproject
+    www.cfdem.com
+                                Christoph Goniva, christoph.goniva@cfdem.com
+                                Copyright 2009-2012 JKU Linz
+                                Copyright 2012-     DCS Computing GmbH, Linz
+-------------------------------------------------------------------------------
+License
+    This file is part of CFDEMcoupling.
+
+    CFDEMcoupling is free software; you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation; either version 3 of the License, or (at your
+    option) any later version.
+
+    CFDEMcoupling is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with CFDEMcoupling; if not, write to the Free Software Foundation,
+    Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
+
+Description
+    This code is designed to realize coupled CFD-DEM simulations using LIGGGHTS
+    and OpenFOAM. Note: this code is not part of OpenFOAM (see DISCLAIMER).
+\*---------------------------------------------------------------------------*/
+
+
+CFDEM coupling provides an open source parallel coupled CFD-DEM framework 
+combining the strengths of LIGGGHTS DEM code and the Open Source 
+CFD package OpenFOAM(R)(*). The CFDEMcoupling toolbox allows to expand 
+standard CFD solvers of OpenFOAM(R)(*) to include a coupling to the DEM 
+code LIGGGHTS. In this toolbox the particle representation within the 
+CFD solver is organized by "cloud" classes. Key functionalities are organised 
+in sub-models (e.g. force models, data exchange models, etc.) which can easily 
+be selected and combined by dictionary settings.
+
+The coupled solvers run fully parallel on distributed-memory clusters. 
+
+Features are:
+
+- its modular approach allows users to easily implement new models
+- its MPI parallelization enables to use it for large scale problems
+- the "forum"_lws on CFD-DEM gives the possibility to exchange with other 
+  users / developers
+- the use of GIT allows to easily update to the latest version
+- basic documentation is provided
+
+The file structure:
+
+- "src" directory including the source files of the coupling toolbox and models
+- "applications" directory including the solver files for coupled CFD-DEM simulations
+- "doc" directory including the documentation of CFDEMcoupling
+- "tutorials" directory including basic tutorial cases showing the functionality
+
+
+
+Details on installation are given on the "www.cfdem.com"
+
+The functionality of this CFD-DEM framwork is described via "tutorial cases" showing 
+how to use different solvers and models.
+
+CFDEMcoupling stands for Computational Fluid Dynamics (CFD) -
+Discrete Element Method (DEM) coupling.
+
+CFDEMcoupling is an open-source code, distributed freely under the terms of the 
+GNU Public License (GPL).
+
+Core development of CFDEMcoupling is done by 
+Christoph Goniva and Christoph Kloss, both at DCS Computing GmbH, 2012
+
+
+\*---------------------------------------------------------------------------*/
+(*) "OpenFOAM(R)"_of is a registered trade mark of the ESI Group.
+This offering is not affiliated, approved or endorsed by ESI Group, 
+the producer of the OpenFOAM® software and owner of the OpenFOAM® trade mark.
+\*---------------------------------------------------------------------------*/
--- a/doc/CFDEMcoupling_Manual.pdf
+++ b/doc/CFDEMcoupling_Manual.pdf
--- a/doc/CFDEMcoupling_Manual.txt
+++ b/doc/CFDEMcoupling_Manual.txt
@ -170,6 +170,38 @@ Reasonable example settings for the "liggghtsCommands" dictionary are given in t



+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
 1.6 Models/Solvers :h4,link(cmd_5),link(comm)

 This section lists all CFDEMcoupling sub-models and solvers alphabetically, with a separate
--- a/doc/IOModel_sophIO.html
+++ b/doc/IOModel_sophIO.html
@ -21,8 +21,7 @@
 </PRE>
 <P><B>Description:</B>
 </P>
-<P>The soph IO-model writes particle positions, velocities, voidfraction, implicit forces, explicit forces , DEM forces and radii to files. the same as for basicIO model. Data is written every write time of the CFD simulation.
-Currently paraview can only interpret position, velocity and radius (a port to VTK might help).
+<P>The sophIO-model is based on basicIO model and additionally writes voidfraction, implicit forces, explicit forces. Data is written every write time of the CFD simulation.
 </P>
 <P><B>Restrictions:</B> None.
 </P>
--- a/doc/IOModel_sophIO.txt
+++ b/doc/IOModel_sophIO.txt
@ -19,8 +19,7 @@ IOModel "sophIO"; :pre

 [Description:]

-The soph IO-model writes particle positions, velocities, voidfraction, implicit forces, explicit forces , DEM forces and radii to files. the same as for basicIO model. Data is written every write time of the CFD simulation.
-Currently paraview can only interpret position, velocity and radius (a port to VTK might help).
+The sophIO-model is based on basicIO model and additionally writes voidfraction, implicit forces, explicit forces. Data is written every write time of the CFD simulation.

 [Restrictions:] None.

--- a/doc/IOModel_trackIO.html
+++ b/doc/IOModel_trackIO.html
@ -21,7 +21,7 @@
 </PRE>
 <P><B>Description:</B>
 </P>
-<P>The basic IO-model writes particle positions, velocities, index and radii to files. The output directories ($casePath/CFD/$timestep) are created automatically. Data is written every write time of the CFD simulation.
+<P>The trackIO-model is based on sophIO model and additionally writes fields necessary to use the particleTracks utility (which needs a particleTrackProperties file in the constant dir). The particleTracks generats tracks of the particles and writes them to a vtk file.
 </P>
 <P><B>Restrictions:</B> None.
 </P>
--- a/doc/IOModel_trackIO.txt
+++ b/doc/IOModel_trackIO.txt
@ -19,7 +19,7 @@ IOModel "trackIO"; :pre

 [Description:]

-The basic IO-model writes particle positions, velocities, index and radii to files. The output directories ($casePath/CFD/$timestep) are created automatically. Data is written every write time of the CFD simulation.
+The trackIO-model is based on sophIO model and additionally writes fields necessary to use the particleTracks utility (which needs a particleTrackProperties file in the constant dir). The particleTracks generats tracks of the particles and writes them to a vtk file.

 [Restrictions:] None.

--- a/doc/forceModel_KochHillDrag.html
+++ b/doc/forceModel_KochHillDrag.html
@ -33,6 +33,8 @@ KochHillDragProps

 <LI><I>interpolation</I> = flag to use interpolated voidfraction and fluid velocity values (normally off) 

+<LI><I>implDEM</I> = flag to use implicit formulation of drag on DEM side (normally off) 
+

 </UL>
 <P><B>Examples:</B>
--- a/doc/forceModel_KochHillDrag.txt
+++ b/doc/forceModel_KochHillDrag.txt
@ -27,6 +27,7 @@ KochHillDragProps
 {density} = name of the finite volume gravity field :l
 {voidfraction} = name of the finite volume voidfraction field :l
 {interpolation} = flag to use interpolated voidfraction and fluid velocity values (normally off) :l
+{implDEM} = flag to use implicit formulation of drag on DEM side (normally off) :l
 :ule

 [Examples:]
--- a/doc/githubAccess_public.pdf
+++ b/doc/githubAccess_public.pdf
--- a/src/lagrangian/cfdemParticle/Make/files
+++ b/src/lagrangian/cfdemParticle/Make/files
@ -6,7 +6,6 @@ voidFractionModels = subModels/voidFractionModel
 locateModels = subModels/locateModel
 meshMotionModels = subModels/meshMotionModel
 momCoupleModels = subModels/momCoupleModel
-regionModels = subModels/regionModel
 dataExchangeModels = subModels/dataExchangeModel
 averagingModels = subModels/averagingModel
 clockModels = subModels/clockModel
@ -21,18 +20,33 @@ $(forceModels)/forceModel/newForceModel.C
 $(forceModels)/noDrag/noDrag.C
 $(forceModels)/checkCouplingInterval/checkCouplingInterval.C
 $(forceModels)/DiFeliceDrag/DiFeliceDrag.C
+$(forceModels)/DiFeliceDragNLift/DiFeliceDragNLift.C
 $(forceModels)/GidaspowDrag/GidaspowDrag.C
 $(forceModels)/SchillerNaumannDrag/SchillerNaumannDrag.C
 $(forceModels)/Archimedes/Archimedes.C
 $(forceModels)/ArchimedesIB/ArchimedesIB.C
 $(forceModels)/interface/interface.C
 $(forceModels)/ShirgaonkarIB/ShirgaonkarIB.C
+$(forceModels)/fieldTimeAverage/fieldTimeAverage.C
+$(forceModels)/fieldBound/fieldBound.C
+$(forceModels)/volWeightedAverage/volWeightedAverage.C
+$(forceModels)/totalMomentumExchange/totalMomentumExchange.C
 $(forceModels)/KochHillDrag/KochHillDrag.C
+$(forceModels)/BeetstraDrag/multiphaseFlowBasic/multiphaseFlowBasic.C
+$(forceModels)/BeetstraDrag/BeetstraDrag.C
+$(forceModels)/LaEuScalarLiquid/LaEuScalarLiquid.C
 $(forceModels)/LaEuScalarTemp/LaEuScalarTemp.C
+$(forceModels)/LaEuScalarDust/LaEuScalarDust.C
 $(forceModels)/virtualMassForce/virtualMassForce.C
 $(forceModels)/gradPForce/gradPForce.C
+$(forceModels)/gradULiftForce/gradULiftForce.C
 $(forceModels)/viscForce/viscForce.C
 $(forceModels)/MeiLift/MeiLift.C
+$(forceModels)/KochHillDragNLift/KochHillDragNLift.C
+$(forceModels)/solidsPressureForce/solidsPressureForce.C
+$(forceModels)/periodicPressure/periodicPressure.C
+$(forceModels)/periodicPressureControl/periodicPressureControl.C
+$(forceModels)/averageSlipVel/averageSlipVel.C

 $(forceModelsMS)/forceModelMS/forceModelMS.C
 $(forceModelsMS)/forceModelMS/newForceModelMS.C
@ -43,6 +57,7 @@ $(IOModels)/IOModel/newIOModel.C
 $(IOModels)/noIO/noIO.C
 $(IOModels)/basicIO/basicIO.C
 $(IOModels)/trackIO/trackIO.C
+$(IOModels)/sophIO/sophIO.C

 $(voidFractionModels)/voidFractionModel/voidFractionModel.C
 $(voidFractionModels)/voidFractionModel/newVoidFractionModel.C
@ -61,21 +76,18 @@ $(locateModels)/turboEngineSearch/turboEngineSearch.C
 $(locateModels)/turboEngineSearchM2M/turboEngineSearchM2M.C
 $(locateModels)/engineSearchIB/engineSearchIB.C

-
 $(meshMotionModels)/meshMotionModel/meshMotionModel.C
 $(meshMotionModels)/meshMotionModel/newMeshMotionModel.C
 $(meshMotionModels)/noMeshMotion/noMeshMotion.C
+$(meshMotionModels)/DEMdrivenMeshMotion/DEMdrivenMeshMotion.C

 $(momCoupleModels)/momCoupleModel/momCoupleModel.C
 $(momCoupleModels)/momCoupleModel/newMomCoupleModel.C
 $(momCoupleModels)/explicitCouple/explicitCouple.C
+$(momCoupleModels)/explicitCoupleSource/explicitCoupleSource.C
 $(momCoupleModels)/implicitCouple/implicitCouple.C
 $(momCoupleModels)/noCouple/noCouple.C

-$(regionModels)/regionModel/regionModel.C
-$(regionModels)/regionModel/newRegionModel.C
-$(regionModels)/allRegion/allRegion.C
-
 $(dataExchangeModels)/dataExchangeModel/dataExchangeModel.C
 $(dataExchangeModels)/dataExchangeModel/newDataExchangeModel.C
 $(dataExchangeModels)/oneWayVTK/oneWayVTK.C
@ -101,4 +113,4 @@ $(liggghtsCommandModels)/runLiggghts/runLiggghts.C
 $(liggghtsCommandModels)/writeLiggghts/writeLiggghts.C
 $(liggghtsCommandModels)/readLiggghtsData/readLiggghtsData.C

-LIB = $(FOAM_USER_LIBBIN)/lib$(CFDEM_LIB_NAME)
+LIB = $(CFDEM_LIB_DIR)/lib$(CFDEM_LIB_NAME)
--- a/src/lagrangian/cfdemParticle/cfdemCloud/cfdemCloud.C
+++ b/src/lagrangian/cfdemParticle/cfdemCloud/cfdemCloud.C
@ -76,9 +76,11 @@ Foam::cfdemCloud::cfdemCloud
    modelType_(couplingProperties_.lookup("modelType")),
    positions_(NULL),
    velocities_(NULL),
+    fluidVel_(NULL),
    impForces_(NULL),
    expForces_(NULL),
    DEMForces_(NULL),
+    Cds_(NULL),
    radii_(NULL),
    voidfractions_(NULL),
    cellIDs_(NULL),
@ -92,6 +94,7 @@ Foam::cfdemCloud::cfdemCloud
    liggghtsCommandModelList_(liggghtsCommandDict_.lookup("liggghtsCommandModels")),
    turbulenceModelType_(couplingProperties_.lookup("turbulenceModelType")),
    cgOK_(true),
+    impDEMdrag_(false),
    turbulence_
    (
        #if defined(version21) || defined(version16ext)
@ -227,9 +230,11 @@ Foam::cfdemCloud::~cfdemCloud()
    clockM().normHist();
    dataExchangeM().destroy(positions_,3);
    dataExchangeM().destroy(velocities_,3);
+    dataExchangeM().destroy(fluidVel_,3);
    dataExchangeM().destroy(impForces_,3);
    dataExchangeM().destroy(expForces_,3);
    dataExchangeM().destroy(DEMForces_,3);
+    dataExchangeM().destroy(Cds_,1);
    dataExchangeM().destroy(radii_,1);
    dataExchangeM().destroy(voidfractions_,1);
    dataExchangeM().destroy(cellIDs_,1);
@ -246,12 +251,16 @@ void Foam::cfdemCloud::getDEMdata()

 void Foam::cfdemCloud::giveDEMdata()
 {
-    for(int index = 0;index <  numberOfParticles(); ++index){
-        for(int i=0;i<3;i++){
-            impForces_[index][i] += expForces_[index][i] + DEMForces_[index][i];
+    if(forceM(0).coupleForce())
+    {
+        dataExchangeM().giveData("dragforce","vector-atom",DEMForces_);
+
+        if(impDEMdrag_)
+        {
+            dataExchangeM().giveData("Ksl","scalar-atom",Cds_);
+            dataExchangeM().giveData("uf","vector-atom",fluidVel_);
        }
    }
-    if(forceM(0).coupleForce()) dataExchangeM().giveData("dragforce","vector-atom",impForces_);
    if(verbose_) Info << "giveDEMdata done." << endl;
 }

@ -275,10 +284,12 @@ void Foam::cfdemCloud::findCells()

 void Foam::cfdemCloud::setForces()
 {
+    resetArray(fluidVel_,numberOfParticles(),3);
    resetArray(impForces_,numberOfParticles(),3);
    resetArray(expForces_,numberOfParticles(),3);
    resetArray(DEMForces_,numberOfParticles(),3);
-    for (int i=0;i<cfdemCloud::nrForceModels();i++) cfdemCloud::forceM(i).setForce(NULL,impForces_,expForces_,DEMForces_);
+    resetArray(Cds_,numberOfParticles(),1);
+    for (int i=0;i<cfdemCloud::nrForceModels();i++) cfdemCloud::forceM(i).setForce();
 }

 // * * * * * * * * * * * * * * * public Member Functions  * * * * * * * * * * * * * //
@ -323,6 +334,13 @@ vector Foam::cfdemCloud::velocity(int index)
    return vel;
 }

+vector Foam::cfdemCloud::fluidVel(int index)
+{
+    vector vel;
+    for(int i=0;i<3;i++) vel[i] = fluidVels()[index][i];
+    return vel;
+}
+
 const forceModel& Foam::cfdemCloud::forceM(int i)
 {
    return forceModel_[i];
@ -491,9 +509,11 @@ bool Foam::cfdemCloud::reAllocArrays() const
        // get arrays of new length
        dataExchangeM().allocateArray(positions_,0.,3);
        dataExchangeM().allocateArray(velocities_,0.,3);
+        dataExchangeM().allocateArray(fluidVel_,0.,3);
        dataExchangeM().allocateArray(impForces_,0.,3);
        dataExchangeM().allocateArray(expForces_,0.,3);
        dataExchangeM().allocateArray(DEMForces_,0.,3);
+        dataExchangeM().allocateArray(Cds_,0.,1);
        dataExchangeM().allocateArray(radii_,0.,1);
        dataExchangeM().allocateArray(voidfractions_,1.,voidFractionM().maxCellsPerParticle());
        dataExchangeM().allocateArray(cellIDs_,0.,voidFractionM().maxCellsPerParticle());
--- a/src/lagrangian/cfdemParticle/cfdemCloud/cfdemCloud.H
+++ b/src/lagrangian/cfdemParticle/cfdemCloud/cfdemCloud.H
@ -101,12 +101,16 @@ protected:

    mutable double **velocities_;

+    mutable double **fluidVel_;
+
    mutable double **impForces_;

    mutable double **expForces_;

    mutable double **DEMForces_;

+    mutable double **Cds_;
+
    mutable double **radii_;

    mutable double **voidfractions_;
@ -133,6 +137,8 @@ protected:

    bool cgOK_;

+    bool impDEMdrag_;
+
    #if defined(version21) || defined(version16ext)
        #ifdef comp
            const compressible::turbulenceModel& turbulence_;
@ -178,6 +184,7 @@ public:

    friend class dataExchangeModel;
    friend class voidFractionModel;
+    friend class forceModel;

 // Constructors

@ -207,6 +214,8 @@ public:

       vector velocity(int);

+       vector fluidVel(int);
+
       virtual const forceModel& forceM(int);

       virtual int nrForceModels();
@ -229,12 +238,16 @@ public:

       inline double ** velocities() const;

+       inline double ** fluidVels() const;
+
       inline double ** impForces() const;

       inline double ** expForces() const;

       inline double ** DEMForces() const;

+       inline double ** Cds() const;
+
       inline double ** radii() const;

       inline double ** voidfractions() const;
--- a/src/lagrangian/cfdemParticle/cfdemCloud/cfdemCloudI.H
+++ b/src/lagrangian/cfdemParticle/cfdemCloud/cfdemCloudI.H
@ -64,6 +64,11 @@ inline double ** cfdemCloud::velocities() const
    return velocities_;
 }

+inline double ** cfdemCloud::fluidVels() const
+{
+    return fluidVel_;
+}
+
 inline double ** cfdemCloud::impForces() const
 {
    return impForces_;
@ -79,6 +84,11 @@ inline double ** cfdemCloud::DEMForces() const
    return DEMForces_;
 }

+inline double ** cfdemCloud::Cds() const
+{
+    return Cds_;
+}
+
 inline double ** cfdemCloud::radii() const
 {
    return radii_;
--- a/src/lagrangian/cfdemParticle/derived/cfdemCloudIB/cfdemCloudIB.C
+++ b/src/lagrangian/cfdemParticle/derived/cfdemCloudIB/cfdemCloudIB.C
@ -113,9 +113,10 @@ bool Foam::cfdemCloudIB::evolve()
            for(int i=0;i<3;i++){
                impForces_[index][i] = 0;
                expForces_[index][i] = 0;
+                DEMForces_[index][i] = 0;
            }
        }
-        for (int i=0;i<nrForceModels();i++) forceM(i).setForce(NULL,impForces_,expForces_,expForces_);
+        for (int i=0;i<nrForceModels();i++) forceM(i).setForce();
        if(verbose_) Info << "setForce done." << endl;

        // write DEM data
--- a/src/lagrangian/cfdemParticle/derived/cfdemCloudMS/cfdemCloudMS.C
+++ b/src/lagrangian/cfdemParticle/derived/cfdemCloudMS/cfdemCloudMS.C
@ -221,7 +221,7 @@ void Foam::cfdemCloudMS::setForces()
    resetArray(impForces_,numberOfParticles(),3);
    resetArray(expForces_,numberOfParticles(),3);
    resetArray(DEMForces_,numberOfParticles(),3);
-    for (int i=0;i<cfdemCloudMS::nrForceModels();i++) cfdemCloudMS::forceM(i).setForce(NULL,impForces_,expForces_,DEMForces_);
+    for (int i=0;i<cfdemCloudMS::nrForceModels();i++) cfdemCloudMS::forceM(i).setForce();
 }

 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
--- a/src/lagrangian/cfdemParticle/etc/controlDict_cgs_2.1.x
+++ b/src/lagrangian/cfdemParticle/etc/controlDict_cgs_2.1.x
@ -0,0 +1,999 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.0.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      controlDict;
+}
+
+// NB: the #functions do not work here
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+Documentation
+{
+    docBrowser      "firefox %f";
+    doxyDocDirs
+    (
+        "$WM_PROJECT_USER_DIR/html"
+        "~OpenFOAM/html"
+        "$WM_PROJECT_DIR/doc/Doxygen/html"
+    );
+    doxySourceFileExts
+    (
+        "App_8C.html"
+        "_8C.html"
+    );
+}
+
+
+InfoSwitches
+{
+    writePrecision  6;
+    writeJobInfo    0;
+    writeDictionaries 0;
+
+    // Allow case-supplied C++ code (#codeStream, codedFixedValue)
+    allowSystemOperations   0;
+}
+
+
+OptimisationSwitches
+{
+    fileModificationSkew 10;
+
+    //- Modification checking:
+    //  - timeStamp         : use modification time on file
+    //  - inotify           : use inotify framework
+    //  - timeStampMaster   : do time stamp (and file reading) only on master.
+    //  - inotifyMaster     : do inotify (and file reading) only on master.
+    fileModificationChecking timeStampMaster;//inotify;timeStamp;inotifyMaster;
+
+    commsType       nonBlocking; //scheduled; //blocking;
+    floatTransfer   0;
+    nProcsSimpleSum 0;
+
+    // Force dumping (at next timestep) upon signal (-1 to disable)
+    writeNowSignal              -1; //10;
+    // Force dumping (at next timestep) upon signal (-1 to disable) and exit
+    stopAtWriteNowSignal        -1;
+}
+
+
+DebugSwitches
+{
+    Analytical          0;
+    APIdiffCoefFunc     0;
+    Ar                  0;
+    BICCG               0;
+    BirdCarreau         0;
+    C10H22              0;
+    C12H26              0;
+    C13H28              0;
+    C14H30              0;
+    C16H34              0;
+    C2H5OH              0;
+    C2H6                0;
+    C2H6O               0;
+    C3H6O               0;
+    C3H8                0;
+    C4H10O              0;
+    C6H14               0;
+    C6H6                0;
+    C7H16               0;
+    C7H8                0;
+    C8H10               0;
+    C8H18               0;
+    C9H20               0;
+    CH3OH               0;
+    CH4N2O              0;
+    CarnahanStarling    0;
+    CallbackRegistry    0;
+    CentredFitData<biLinearFitPolynomial>           0;
+    CentredFitData<linearFitPolynomial>             0;
+    CentredFitData<quadraticFitPolynomial>          0;
+    CentredFitData<quadraticLinearFitPolynomial>    0;
+    ChomiakInjector     0;
+    Cloud<basicKinematicParcel> 0;
+    Cloud<basicReactingParcel>  0;
+    Cloud<basicThermoParcel>    0;
+    Cloud<indexedParticle>  0;
+    Cloud<parcel>           0;
+    Cloud<passiveParticle>  0;
+    Cloud<trackedParticle>  0;
+    CoEuler             0;
+    CompositionModel    0;
+    ConeInjection       0;
+    Constant            0;
+    ConstantRateDevolatilisation    0;
+    CrankNicholson      0;
+    CrossPowerLaw       0;
+    Cs                  0;
+    DIC                 0;
+    DICGaussSeidel      0;
+    DILU                0;
+    DILUGaussSeidel     0;
+    DataEntry           0;
+    DeardorffDiffStress 0;
+    DispersionModel     0;
+    DispersionRASModel  0;
+    DragModel           0;
+    ETAB                0;
+    Ergun               0;
+    Euler               0;
+    EulerImplicit       0;
+    EulerRotation       0;
+    extendedCellToFaceStencil 0;
+    FDIC                0;
+    FaceCellWave        0;
+    GAMG                0;
+    GAMGAgglomeration   0;
+    GAMGInterface       0;
+    GAMGInterfaceField  0;
+    Gamma               0;
+    Gamma01             0;
+    GammaV              0;
+    Gauss               0;
+    GaussSeidel         0;
+    Gibilaro            0;
+    Gidaspow            0;
+    GidaspowErgunWenYu  0;
+    GidaspowSchillerNaumann 0;
+    GradientDispersionRAS   0;
+    Gulders             0;
+    GuldersEGR          0;
+    H2O                 0;
+    HashTable           0;
+    HeatTransferModel   0;
+    HerschelBulkley     0;
+    HrenyaSinclair      0;
+    IC8H18              0;
+    ICCG                0;
+    IDDESDelta          0;
+    IDEA                0;
+    IFstream            0;
+    IOMap<dictionary>   0;
+    IOPtrList<MRFZone>  0;
+    IOPtrList<coordinateSystem> 0;
+    IOPtrList<injector> 0;
+    IOPtrList<porousZone>   0;
+    IOobject            0;
+    InjectionModel      0;
+    IntegrationScheme   0;
+    JohnsonJackson      0;
+    KRR4                0;
+    KinematicCloud<basicKinematicParcel> 0;
+    KinematicCloud<basicReactingParcel> 0;
+    KinematicCloud<basicThermoParcel>   0;
+    KinematicParcel<basicReactingParcel> 0;
+    KinematicParcel<basicThermoParcel>  0;
+    LESModel            0;
+    LESdelta            0;
+    LESfilter           0;
+    LISA                0;
+    LRR                 0;
+    LRRDiffStress       0;
+    LamBremhorstKE      0;
+    LaunderGibsonRSTM   0;
+    LaunderSharmaKE     0;
+    LienCubicKE         0;
+    LienCubicKELowRe    0;
+    LienLeschzinerLowRe 0;
+    MB                  0;
+    MC                  0;
+    MCV                 0;
+    MUSCL               0;
+    MUSCL01             0;
+    MUSCLV              0;
+    ManualInjection     0;
+    MarshakRadiation    0;
+    MarshakRadiationFixedT  0;
+    MassTransferModel   0;
+    MeshWave            0;
+    Minmod              0;
+    MinmodV             0;
+    N2                  0;
+    NSRDSfunc0          0;
+    NSRDSfunc1          0;
+    NSRDSfunc14         0;
+    NSRDSfunc2          0;
+    NSRDSfunc3          0;
+    NSRDSfunc4          0;
+    NSRDSfunc5          0;
+    NSRDSfunc6          0;
+    NSRDSfunc7          0;
+    Newtonian           0;
+    NoDispersion        0;
+    NoDrag              0;
+    NoHeatTransfer      0;
+    NoInjection         0;
+    NoMassTransfer      0;
+    NoSurfaceReaction   0;
+    NonlinearKEShih     0;
+    ODE                 0;
+    ODESolver           0;
+    OFstream            0;
+    ORourke             0;
+    OSPRE               0;
+    OSPREV              0;
+    P1                  0;
+    PBiCG               0;
+    PCG                 0;
+    PackedList          0;
+    ParSortableList     0;
+    PatchToPatchInterpolation 0;
+    Phi                 0;
+    PointEdgeWave       0;
+    POSIX               0;
+    Prandtl             0;
+    PrimitivePatch      0;
+    Pstream             0;
+    QUICK               0;
+    QUICKV              0;
+    QZeta               0;
+    RASModel            0;
+    RK                  0;
+    RNGkEpsilon         0;
+    RanzMarshall        0;
+    ReactingCloud<basicReactingParcel>  0;
+    ReactingParcel<basicReactingParcel> 0;
+    Rebound             0;
+    ReitzDiwakar        0;
+    ReitzKHRT           0;
+    RosinRammler        0;
+    RutlandFlashBoil    0;
+    SDA                 0;
+    SFCD                0;
+    SFCDV               0;
+    SHF                 0;
+    SIBS                0;
+    SKA                 0;
+    SLTS                0;
+    SRFModel            0;
+    SRFVelocity         0;
+    STARCDRotation      0;
+    Schaeffer           0;
+    SchillerNaumann     0;
+    SinclairJackson     0;
+    SingleKineticRateDevolatilisation 0;
+    SingleMixtureFraction   0;
+    Smagorinsky         0;
+    SpalartAllmaras     0;
+    SpalartAllmarasDDES 0;
+    SpalartAllmarasIDDES    0;
+    SphereDrag          0;
+    StandardWallInteraction 0;
+    StaticHashTable     0;
+    StochasticDispersionRAS 0;
+    SuperBee            0;
+    SuperBeeV           0;
+    SurfaceReactionModel 0;
+    Syamlal             0;
+    SyamlalOBrien       0;
+    SyamlalRogersOBrien 0;
+    TAB                 0;
+    Table               0;
+    ThermoCloud<basicReactingParcel>    0;
+    ThermoCloud<basicThermoParcel>      0;
+    ThermoParcel<basicReactingParcel>   0;
+    ThermoParcel<basicThermoParcel>     0;
+    UMIST               0;
+    UMISTV              0;
+    UpwindFitData<cubicUpwindFitPolynomial> 0;
+    UpwindFitData<quadraticLinearUpwindFitPolynomial> 0;
+    UpwindFitData<quadraticUpwindFitPolynomial>     0;
+    WallInteractionModel 0;
+    WenYu               0;
+    aC11H10             0;
+    absorptionEmissionModel 0;
+    addCell             0;
+    addFace             0;
+    addPatchCellLayer   0;
+    addPoint            0;
+    advective           0;
+    algebraicPair       0;
+    alphaContactAngle   0;
+    alphaFixedPressure  0;
+    alphatWallFunction  0;
+    angularOscillatingDisplacement  0;
+    angularOscillatingVelocity  0;
+    anisotropic         0;
+    ash                 0;
+    atomizationModel    0;
+    attachDetach        0;
+    autoDensity         0;
+    autoHexMeshDriver   0;
+    autoLayerDriver     0;
+    autoRefineDriver    0;
+    autoSnapDriver      0;
+    bC11H10             0;
+    backgroundMeshDecomposition 0;
+    backward            0;
+    basePatch           0;
+    basicKinematicCloud 0;
+    basicKinematicParcel 0;
+    basicMixture        0;
+    basicReactingCloud  0;
+    basicReactingParcel 0;
+    basicThermo         0;
+    basicThermoCloud    0;
+    basicThermoParcel   0;
+    biLinearFit         0;
+    binaryAbsorptionEmission 0;
+    blended             0;
+    blobsSheetAtomization   0;
+    blobsSwirlInjector  0;
+    booleanSurface      0;
+    boundaryCutter      0;
+    boundaryMesh        0;
+    boundaryToFace      0;
+    boundedBackward     0;
+    boxToCell           0;
+    boxToFace           0;
+    boxToPoint          0;
+    breakupModel        0;
+    calculated          0;
+    cell                0;
+    cellClassification  0;
+    cellCuts            0;
+    cellDistFuncs       0;
+    cellLimited         0;
+    cellList            0;
+    cellLooper          0;
+    cellMDLimited       0;
+    cellMotion          0;
+    cellPoint           0;
+    cellPointFace       0;
+    cellPointWeight     0;
+    cellSet             0;
+    cellSizeControlSurfaces 0;
+    cellToCell          0;
+    cellToFace          0;
+    cellToPoint         0;
+    cellZone            0;
+    centredCECStencil   0;
+    centredCFCStencil   0;
+    chemistryReader     0;
+    chemistrySolver     0;
+    chemkinReader       0;
+    clippedLinear       0;
+    cloud               0;
+    cloudAbsorptionEmission 0;
+    cloudScatter        0;
+    collisionModel      0;
+    combineFaces        0;
+    commSchedule        0;
+    commonRailInjector  0;
+    compound            0;
+    constInjector       0;
+    constant            0;
+    constantAbsorptionEmission  0;
+    constantAlphaContactAngle   0;
+    constantScatter     0;
+    coordinateRotation  0;
+    coordinateSystem    0;
+    coordinateSystems   0;
+    corrected           0;
+    coupled             0;
+    cubeRootVol         0;
+    cubic               0;
+    cubicUpwindFit      0;
+    curve               0;
+    cyclic              0;
+    cyclicLduInterface  0;
+    cyclicLduInterfaceField 0;
+    cylinderToCell      0;
+    cylindrical         0;
+    decompositionMethod 0;
+    definedHollowConeInjector 0;
+    definedInjector     0;
+    definedPressureSwirlInjector 0;
+    diagTensorField     0;
+    diagonal            0;
+    dictionary          0;
+    dimensionSet        1;
+    mappedBase    0;
+    mappedPatch   0;
+    mappedVelocityFlux 0;
+    directionMixed      0;
+    directional         0;
+    disallowGenericFvPatchField 0;
+    disallowGenericPointPatchField 0;
+    disallowGenericPolyPatch    0;
+    dispersionLESModel  0;
+    dispersionModel     0;
+    dispersionRASModel  0;
+    displacementComponentLaplacian 0;
+    displacementInterpolation 0;
+    displacementLaplacian 0;
+    displacementSBRStress 0;
+    distanceSurface     0;
+    Distribution        0;
+    downwind            0;
+    dragModel           0;
+    duplicatePoints     0;
+    dx                  0;
+    dynMixedSmagorinsky 0;
+    dynOneEqEddy        0;
+    dynSmagorinsky      0;
+    dynamicAlphaContactAngle 0;
+    dynamicFvMesh       0;
+    dynamicInkJetFvMesh 0;
+    dynamicMotionSolverFvMesh 0;
+    dynamicRefineFvMesh 0;
+    edgeIntersections   0;
+    edgeList            0;
+    edgeSurface         0;
+    empty               0;
+    engineMesh          0;
+    enrichedPatch       0;
+    epsilonWallFunction 0;
+    errorDrivenRefinement   0;
+    evaporationModel    0;
+    exponential         0;
+    extendedLeastSquares    0;
+    extendedLeastSquaresVectors 0;
+    face                0;
+    faceAreaPair        0;
+    faceCoupleInfo      0;
+    faceLimited         0;
+    faceList            0;
+    faceMDLimited       0;
+    faceSet             0;
+    faceToCell          0;
+    faceToFace          0;
+    faceToPoint         0;
+    faceZone            0;
+    fan                 0;
+    featureEdgeMesh     0;
+    fieldToCell         0;
+    file                0;
+    fileName            2;
+    filteredLinear      0;
+    filteredLinear2     0;
+    filteredLinear2V    0;
+    filteredLinear3     0;
+    filteredLinear3V    0;
+    fixedEnthalpy       0;
+    buoyantPressure 0;
+    fixedFluxBoussinesqBuoyantPressure 0;
+    fixedFluxPressure   0;
+    fixedGradient       0;
+    fixedInternalEnergy 0;
+    fixedInternalValue  0;
+    fixedNormalSlip     0;
+    fixedPressureCompressibleDensity 0;
+    fixedUnburntEnthalpy 0;
+    fixedValue          0;
+    flowRateInletVelocity   0;
+    fluxCorrectedVelocity   0;
+    foamChemistryReader 0;
+    foamFile            0;
+    forceCoeffs         0;
+    forces              0;
+    fourth              0;
+    freestream          0;
+    freestreamPressure  0;
+    frictionalStressModel   0;
+    functionObject      0;
+    fv                  0;
+    fvMesh              0;
+    fvMeshDistribute    0;
+    fvMotionSolver      0;
+    fvPatchField        0;
+    fvScalarMatrix      0;
+    fvSchemes           0;
+    fvSphericalTensorMatrix 0;
+    fvSymmTensorMatrix  0;
+    fvTensorMatrix      0;
+    fvVectorMatrix      0;
+    fvsPatchField       0;
+    general             0;
+    generic             0;
+    genericPatch        0;
+    geomCellLooper      0;
+    geometricSurfacePatch   0;
+    global              0;
+    globalIndexAndTransform 0;
+    globalMeshData      0;
+    globalPoints        0;
+    gnuplot             0;
+    gradientDispersionRAS   0;
+    gradientEnthalpy        0;
+    gradientInternalEnergy  0;
+    gradientUnburntEnthalpy 0;
+    granularPressureModel   0;
+    hCombustionThermo   0;
+    hMixtureThermo<dieselMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>> 0;
+    hMixtureThermo<homogeneousMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>> 0;
+    hMixtureThermo<homogeneousMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>> 0;
+    hMixtureThermo<inhomogeneousMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>> 0;
+    hMixtureThermo<inhomogeneousMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>> 0;
+    hMixtureThermo<multiComponentMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>> 0;
+    hMixtureThermo<reactingMixture> 0;
+    hMixtureThermo<veryInhomogeneousMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>> 0;
+    hMixtureThermo<veryInhomogeneousMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>> 0;
+    hThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>> 0;
+    hThermo<pureMixture<sutherlandTransport<specieThermo<hConstThermo<perfectGas>>>>> 0;
+    hThermo<pureMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>> 0;
+    harmonic            0;
+    heatTransferModel   0;
+    hexCellLooper       0;
+    hexRef8             0;
+    hhuCombustionThermo 0;
+    hhuMixtureThermo<egrMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>> 0;
+    hhuMixtureThermo<egrMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>> 0;
+    hhuMixtureThermo<homogeneousMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>> 0;
+    hhuMixtureThermo<homogeneousMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>> 0;
+    hhuMixtureThermo<inhomogeneousMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>> 0;
+    hhuMixtureThermo<inhomogeneousMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>> 0;
+    hhuMixtureThermo<veryInhomogeneousMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>> 0;
+    hhuMixtureThermo<veryInhomogeneousMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>> 0;
+    hierarchical        0;
+    hollowConeInjector  0;
+    iC3H8O              0;
+    indexedOctree       0;
+    indexedParticle     0;
+    injectorModel       0;
+    injectorType        0;
+    inletOutlet         0;
+    inletOutletTotalTemperature 0;
+    interfaceCompression 0;
+    intersectedSurface  0;
+    inverseDistance     0;
+    inverseFaceDistance 0;
+    inversePointDistance 0;
+    inverseVolume       0;
+    isoSurface          0;
+    isoSurfaceCell      0;
+    jplot               0;
+    jumpCyclic          0;
+    kEpsilon            0;
+    kOmega              0;
+    kOmegaSST           0;
+    kOmegaSSTSAS        0;
+    kqRWallFunction     0;
+    kinematicCloud      0;
+    labelField          0;
+    labelList           0;
+    labelListList       0;
+    labelToCell         0;
+    labelToFace         0;
+    labelToPoint        0;
+    laminar             0;
+    laminarFlameSpeed   0;
+    laplace             0;
+    layerAdditionRemoval 0;
+    layered             0;
+    lduInterface        0;
+    lduInterfaceField   0;
+    lduMatrix           1;
+    lduMesh             0;
+    leastSquares        0;
+    leastSquaresVectors 0;
+    level               2;
+    limitWith           0;
+    limited             0;
+    limitedCubic        0;
+    limitedCubic01      0;
+    limitedCubicV       0;
+    limitedGamma        0;
+    limitedLimitedCubic 0;
+    limitedLimitedLinear 0;
+    limitedLinear       0;
+    limitedLinear01     0;
+    limitedLinearV      0;
+    limitedMUSCL        0;
+    limitedSurfaceInterpolationScheme 0;
+    limitedVanLeer      0;
+    linear              0;
+    linearFit           0;
+    linearUpwind        0;
+    linearUpwindV       0;
+    liquid              0;
+    locDynOneEqEddy     0;
+    localBlended        0;
+    localMax            0;
+    localMin            0;
+    localPointRegion    0;
+    lowReOneEqEddy      0;
+    manual              0;
+    meshCutAndRemove    0;
+    meshCutter          0;
+    meshModifier        0;
+    meshRefinement      0;
+    meshSearch          0;
+    meshToMesh          0;
+    metis               0;
+    midPoint            0;
+    midPointAndFace     0;
+    mixed               0;
+    mixedEnthalpy       0;
+    mixedInternalEnergy 0;
+    mixedSmagorinsky    0;
+    mixedUnburntEnthalpy 0;
+    mixerFvMesh         0;
+    modifyCell          0;
+    modifyFace          0;
+    modifyPoint         0;
+    motionDiffusivity   0;
+    motionDirectional   0;
+    motionSmoother      0;
+    motionSolver        0;
+    movingConeTopoFvMesh 0;
+    movingWallVelocity  0;
+    muSgsSpalartAllmarasWallFunction 0;
+    multiDirRefinement  0;
+    multiHoleInjector   0;
+    multiLevel          1;
+    multivariateSelection 0;
+    mutRoughWallFunction 0;
+    mutSpalartAllmarasStandardRoughWallFunction  0;
+    mutSpalartAllmarasStandardWallFunction  0;
+    mutSpalartAllmarasWallFunction  0;
+    mutWallFunction     0;
+    nC3H8O              0;
+    nbrToCell           0;
+    nearestToCell       0;
+    nearestToPoint      0;
+    noAbsorptionEmission 0;
+    noDragModel         0;
+    noRadiation         0;
+    none                0;
+    normal              0;
+    normalToFace        0;
+    nuSgsSpalartAllmarasWallFunction 0;
+    nutRoughWallFunction 0;
+    nutSpalartAllmarasStandardRoughWallFunction  0;
+    nutSpalartAllmarasStandardWallFunction  0;
+    nutSpalartAllmarasWallFunction  0;
+    nutWallFunction     0;
+    obj                 0;
+    objectRegistry      0;
+    off                 0;
+    omegaWallFunction   0;
+    oneEqEddy           0;
+    orientedSurface     0;
+    oscillatingDisplacement 0;
+    oscillatingFixedValue 0;
+    oscillatingVelocity 0;
+    outletInlet         0;
+    outletStabilised    0;
+    pair                0;
+    parabolicCylindrical 0;
+    parcel              0;
+    partialSlip         0;
+    passiveParticle     0;
+    patch               0;
+    patchToFace         0;
+    patchZones          0;
+    pdf                 0;
+    perfectInterface    0;
+    pointIndexHitList   0;
+    pointPatchField     0;
+    pointScalarField    0;
+    pointScalarField::DimensionedInternalField 0;
+    pointSet            0;
+    pointSphericalTensorField 0;
+    pointSphericalTensorField::DimensionedInternalField 0;
+    pointSymmTensorField 0;
+    pointSymmTensorField::DimensionedInternalField 0;
+    pointTensorField    0;
+    pointTensorField::DimensionedInternalField 0;
+    pointToCell         0;
+    pointToFace         0;
+    pointToPoint        0;
+    pointVectorField    0;
+    pointVectorField::DimensionedInternalField 0;
+    pointZone           0;
+    polyBoundaryMesh    0;
+    polyMesh            0;
+    polyMeshGeometry    0;
+    polyMeshInfo        0;
+    polyTopoChange      0;
+    polyTopoChanger     0;
+    powerLaw            0;
+    pressureDirectedInletOutletVelocity 0;
+    pressureDirectedInletVelocity 0;
+    pressureInletOutletVelocity 0;
+    pressureInletUniformVelocity 0;
+    pressureInletVelocity 0;
+    pressureNormalInletOutletVelocity 0;
+    pressureSwirlInjector 0;
+    primitiveMesh       0;
+    primitiveMeshGeometry 0;
+    probes              0;
+    processor           0;
+    processorLduInterface 0;
+    processorLduInterfaceField 0;
+    pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>> 0;
+    pureMixture<sutherlandTransport<specieThermo<hConstThermo<perfectGas>>>> 0;
+    pureMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>> 0;
+    quadratic           0;
+    quadraticFit        0;
+    quadraticLinearFit  0;
+    quadraticLinearUpwindFit 0;
+    quadraticUpwindFit  0;
+    radiationModel      0;
+    raw                 0;
+    reactingCloud       0;
+    reaction            0;
+    realizableKE        0;
+    refinementHistory   0;
+    refinementIterator  0;
+    reflect             0;
+    regIOobject         0;
+    regionSplit         0;
+    regionToCell        0;
+    remove              0;
+    removeCell          0;
+    removeCells         0;
+    removeFace          0;
+    removeFaces         0;
+    removePoint         0;
+    removePoints        0;
+    reverseLinear       0;
+    rotatedBoxToCell    0;
+    rotatingPressureInletOutletVelocity 0;
+    rotatingTotalPressure 0;
+    sampledPatch        0;
+    sampledPlane        0;
+    sampledSet          0;
+    sampledSurface      0;
+    saturateEvaporationModel 0;
+    scalarAverageField  0;
+    scalarField         0;
+    scalarRange         0;
+    scaleSimilarity     0;
+    scatterModel        0;
+    searchableBox       0;
+    searchableSurface   0;
+    sequential          0;
+    setUpdater          0;
+    sets                0;
+    shapeList           0;
+    shapeToCell         0;
+    simple              0;
+    sixDoFRigidBodyMotionConstraint 0;
+    skewCorrected       0;
+    skewCorrectionVectors 0;
+    sliced              0;
+    slidingInterface    0;
+    slip                0;
+    smooth              0;
+    smoothSolver        0;
+    solid               0;
+    solidBodyMotionFunction 0;
+    solidBodyMotionFvMesh 0;
+    solution            0;
+    spectEddyVisc       0;
+    sphereToCell        0;
+    spherical           0;
+    sphericalTensorAverageField 0;
+    sphericalTensorField 0;
+    standardDragModel   0;
+    standardEvaporationModel 0;
+    staticFvMesh        0;
+    steadyState         0;
+    stl                 0;
+    string              0;
+    stochasticDispersionRAS 0;
+    supersonicFreestream 0;
+    surfaceFeatures     0;
+    surfaceInterpolation 0;
+    surfaceInterpolationScheme 0;
+    surfaceIntersection 0;
+    surfaceNormalFixedValue 0;
+    surfacePatch        0;
+    surfacePatchIOList  0;
+    surfaceScalarField  0;
+    surfaceScalarField::DimensionedInternalField 0;
+    surfaceSlipDisplacement 0;
+    surfaceSphericalTensorField 0;
+    surfaceSphericalTensorField::DimensionedInternalField 0;
+    surfaceSymmTensorField 0;
+    surfaceSymmTensorField::DimensionedInternalField 0;
+    surfaceTensorField 0;
+    surfaceTensorField::DimensionedInternalField 0;
+    surfaceToCell       0;
+    surfaceToPoint      0;
+    surfaceVectorField  0;
+    surfaceVectorField::DimensionedInternalField 0;
+    surfaceWriter       0;
+    surfaces            0;
+    swirlInjector       0;
+    symmTensorAverageField 0;
+    symmTensorField     0;
+    symmetryPlane       0;
+    syringePressure     0;
+    tensorAverageField  0;
+    tensorField         0;
+    tetDecomposedPolyMesh 0;
+    thermoCloud         0;
+    thermophysicalFunction 0;
+    time                0;
+    timeVaryingAlphaContactAngle 0;
+    timeVaryingFlowRateInletVelocity 0;
+    timeVaryingMappedFixedValue 0;
+    timeVaryingTotalPressure 0;
+    timeVaryingUniformFixedValue 0;
+    timer               0;
+    topoAction          0;
+    topoCellLooper      0;
+    topoChangerFvMesh   0;
+    topoSet             0;
+    topoSetSource       0;
+    toroidal            0;
+    totalPressure       0;
+    totalTemperature    0;
+    trackedParticle     0;
+    trajectory          0;
+    transform           0;
+    treeDataCell        0;
+    treeDataFace        0;
+    treeDataTriSurface  0;
+    treeLeaf            0;
+    treeNode            0;
+    triSurface          0;
+    triSurfaceMesh      0;
+    turbulenceModel     0;
+    turbulentHeatFluxTemperature 0;
+    turbulentInlet      0;
+    turbulentIntensityKineticEnergyInlet 0;
+    turbulentMixingLengthDissipationRateInlet 0;
+    turbulentMixingLengthFrequencyInlet 0;
+    uncorrected         0;
+    undoableMeshCutter  0;
+    uniform             0;
+    uniformFixedValue   0;
+    unitInjector        0;
+    upwind              0;
+    upwindCFCStencil    0;
+    value               0;
+    vanAlbada           0;
+    vanAlbadaV          0;
+    vanDriest           0;
+    vanLeer             0;
+    vanLeer01           0;
+    vanLeerV            0;
+    vector2DField       0;
+    vectorAverageField  0;
+    vectorField         0;
+    velocityComponentLaplacian 0;
+    velocityLaplacian   0;
+    viscosityModel      0;
+    volPointInterpolation 0;
+    volScalarField      0;
+    volScalarField::DimensionedInternalField 0;
+    volSphericalTensorField 0;
+    volSphericalTensorField::DimensionedInternalField 0;
+    volSymmTensorField  0;
+    volSymmTensorField::DimensionedInternalField 0;
+    volTensorField      0;
+    volTensorField::DimensionedInternalField 0;
+    volVectorField      0;
+    volVectorField::DimensionedInternalField 0;
+    vtk                 0;
+    walkPatch           0;
+    wall                0;
+    wallHeatTransfer    0;
+    wallLayerCells      0;
+    wallModel           0;
+    waveTransmissive    0;
+    wedge               0;
+    weighted            0;
+    word                2;
+    writer              0;
+    xmgr                0;
+    zeroGradient        0;
+    zoneToCell          0;
+    zoneToFace          0;
+    zoneToPoint         0;
+}
+
+
+DimensionedConstants
+{
+    unitSet             CGS; // SI // USCS
+
+    CGSCoeffs
+    {
+        universal
+        {
+            c               c [ 0 1 -1 0 0 0 0 ] 2.99792e+10; // speed of light in vacuum (cm/s)
+            G               G [ -1 3 -2 0 0 0 0 ] 6.67429e-8; // gravitational constant (cm^3/(gs^2))
+            h               h [ 1 2 -1 0 0 0 0 ] 6.62607e-27; // Planck's constant (erg.s)
+        }
+        electromagnetic
+        {
+            e               e [ 0 0 1 0 0 1 0 ] 4.803204e-10; // elementary charge (statcoulomb)
+        }
+        atomic
+        {
+            me              me [ 1 0 0 0 0 0 0 ] 9.10938e-28; // electron mass (g)
+            mp              mp [ 1 0 0 0 0 0 0 ] 1.67262e-24; // proton mass (g)
+        }
+        physicoChemical
+        {
+            mu              mu [ 1 0 0 0 0 0 0 ] 1.66054e-24; // atomic mass unit (g)
+            k               k [ 1 2 -2 -1 0 0 0 ] 1.38065e-16; // Boltzman constant (erg/K)
+        }
+        standard
+        {
+            //- Standard pressure [Pa]
+            Pstd            Pstd [ 1 -1 -2 0 0 0 0 ] 1000000; // 1 bar (barye)
+            //- Standard temperature [degK]
+            Tstd            Tstd [ 0 0 0 1 0 0 0 ] 298.15; // should be same as in SI unit system
+        }
+    }
+
+    SICoeffs
+    {
+        universal
+        {
+            c               c [ 0 1 -1 0 0 0 0 ] 2.99792e+08;
+            G               G [ -1 3 -2 0 0 0 0 ] 6.67429e-11;
+            h               h [ 1 2 -1 0 0 0 0 ] 6.62607e-34;
+        }
+        electromagnetic
+        {
+            e               e [ 0 0 1 0 0 1 0 ] 1.60218e-19;
+        }
+        atomic
+        {
+            me              me [ 1 0 0 0 0 0 0 ] 9.10938e-31;
+            mp              mp [ 1 0 0 0 0 0 0 ] 1.67262e-27;
+        }
+        physicoChemical
+        {
+            mu              mu [ 1 0 0 0 0 0 0 ] 1.66054e-27;
+            k               k [ 1 2 -2 -1 0 0 0 ] 1.38065e-23;
+        }
+        standard
+        {
+            //- Standard pressure [Pa]
+            Pstd            Pstd [ 1 -1 -2 0 0 0 0 ] 100000;
+            //- Standard temperature [degK]
+            Tstd            Tstd [ 0 0 0 1 0 0 0 ] 298.15;
+        }
+    }
+    USCSCoeffs
+    {
+        universal
+        {
+            c               c [ 0 1 -1 0 0 0 0 ] 9.83558e+08;
+            G               G [ -1 3 -2 0 0 0 0 ] 1.06909e-09;
+            h               h [ 1 2 -1 0 0 0 0 ] 1.57234e-32;
+        }
+        electromagnetic
+        {
+            e               e [ 0 0 1 0 0 1 0 ] 1.60218e-19;
+        }
+        atomic
+        {
+            me              me [ 1 0 0 0 0 0 0 ] 2.00825e-30;
+            mp              mp [ 1 0 0 0 0 0 0 ] 3.68746e-27;
+        }
+        physicoChemical
+        {
+            mu              mu [ 1 0 0 0 0 0 0 ] 3.66083e-27;
+            k               k [ 1 2 -2 -1 0 0 0 ] 1.82012e-22;
+        }
+        standard
+        {
+            //- Standard pressure [lbm/ft^2]
+            Pstd            Pstd [ 1 -1 -2 0 0 0 0 ] 2088.6;
+            //- Standard temperature [degR]
+            Tstd            Tstd [ 0 0 0 1 0 0 0 ] 536.67;
+        }
+    }
+}
+
+
+// ************************************************************************* //
--- a/src/lagrangian/cfdemParticle/subModels/IOModel/basicIO/basicIO.C
+++ b/src/lagrangian/cfdemParticle/subModels/IOModel/basicIO/basicIO.C
@ -146,7 +146,14 @@ void basicIO::streamDataToPath(fileName path, double** array,word name,word type
    *fileStream << "FoamFile\n";
    *fileStream << "{version 2.0; format ascii;class "<< className << "; location 0;object  "<< name <<";}\n";
    *fileStream << nPProc_ <<"\n";
-    *fileStream << "(\n";
+    //*fileStream << "(\n";
+
+    if(type!="origProcId")*fileStream << "(\n";
+    else if(type=="origProcId")
+    {
+        if(nPProc_>0) *fileStream <<"{0}"<< "\n";
+        else *fileStream <<"{}"<< "\n";
+    }

    for(int index = 0;index < particleCloud_.numberOfParticles(); ++index)
    {
@ -154,13 +161,16 @@ void basicIO::streamDataToPath(fileName path, double** array,word name,word type
        {
            if (type=="scalar"){
                *fileStream << array[index][0] << " \n";
-            }else {
+            }else if (type=="position" || type=="vector"){
                for(int i=0;i<3;i++) vec[i] = array[index][i];
                *fileStream <<"( "<< vec[0] <<" "<<vec[1]<<" "<<vec[2]<<" ) "<< finaliser << " \n";
+            }else if (type=="label"){
+                *fileStream << index << finaliser << " \n";
            }
        }
    }
-    *fileStream << ")\n";
+    //*fileStream << ")\n";
+    if(type!="origProcId")*fileStream << ")\n";
    delete fileStream;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
--- a/src/lagrangian/cfdemParticle/subModels/IOModel/trackIO/trackIO.C
+++ b/src/lagrangian/cfdemParticle/subModels/IOModel/trackIO/trackIO.C
@ -60,16 +60,8 @@ trackIO::trackIO
    cfdemCloud& sm
 )
 :
-    IOModel(dict,sm),
-    //propsDict_(dict.subDict(typeName + "Props")),
-    dirName_(""),
-    path_("dev/null"),
-    lagPath_("dev/null")
-{
-    //if (propsDict_.found("dirName")) dirName_=word(propsDict_.lookup("dirName"));
-    path_ = buildFilePath(dirName_);
-
-}
+    sophIO(dict,sm)
+{}

 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //

@ -85,60 +77,16 @@ void trackIO::dumpDEMdata() const
 {
    if (time_.outputTime())
    {
-        // make time directory
-        lagPath_=createTimeDir(path_);
-        lagPath_=createLagrangianDir(fileName(lagPath_));
+        sophIO::dumpDEMdata();

        // stream data to file
-        streamDataToPath(lagPath_, particleCloud_.positions(), particleCloud_.numberOfParticles(), "positions","position","Cloud<passiveParticle>","0");
-        streamDataToPath(lagPath_, particleCloud_.velocities(), particleCloud_.numberOfParticles(), "v","vector","vectorField","");
-        streamDataToPath(lagPath_, particleCloud_.velocities(), particleCloud_.numberOfParticles(), "origId","label","labelField","");
-        streamDataToPath(lagPath_, particleCloud_.velocities(), particleCloud_.numberOfParticles(), "origProcId","origProcId","labelField","");
-        streamDataToPath(lagPath_, particleCloud_.radii(), particleCloud_.numberOfParticles(), "r","scalar","scalarField","");
+        streamDataToPath(lagPath_, particleCloud_.velocities(), "origId","label","labelField","");
+        streamDataToPath(lagPath_, particleCloud_.velocities(), "origProcId","origProcId","labelField","");
    }
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 // Private Member Functions

-fileName trackIO::buildFilePath(word dirName) const
-{
-    // create file structure
-    fileName path("."/dirName);
-    OFstream* stubFile = new OFstream(fileName(path/"particles.foam"));
-    delete stubFile;
-    return path;
-}
-
-void trackIO::streamDataToPath(fileName path, double** array,int n,word name,word type,word className,word finaliser) const
-{
-    vector vec;
-    OFstream* fileStream = new OFstream(fileName(path/name));
-    *fileStream << "FoamFile\n";
-    *fileStream << "{version 2.0; format ascii;class "<< className << "; location 0;object  "<< name <<";}\n";
-    *fileStream << n <<"\n";
-    if(type!="origProcId")*fileStream << "(\n";
-    else if(type=="origProcId")*fileStream <<"{0}"<< "\n";
-
-    for(int index = 0;index < n; ++index)
-    {
-        if (type=="scalar"){
-            *fileStream << array[index][0] << finaliser << " \n";
-        }else if (type=="position"){
-            for(int i=0;i<3;i++) vec[i] = array[index][i];
-//          You may need to use these two lines if you have cyclics
-//		    if(vec[0]<0)vec[0]+=0.12;if(vec[0]>0.12)vec[0]-=0.12;
-//		    if(vec[1]<0)vec[1]+=0.06;if(vec[1]>0.06)vec[1]-=0.06;
-            *fileStream <<"( "<< vec[0] <<" "<<vec[1]<<" "<<vec[2]<<" ) "<< finaliser << " \n";
-        }else if (type=="label"){
-            *fileStream << index << finaliser << " \n";
-        }else  if (type=="vector"){
-            for(int i=0;i<3;i++)  vec[i] = array[index][i]; 
-            *fileStream <<"( "<< vec[0] <<" "<<vec[1]<<" "<<vec[2]<<" ) " << finaliser << " \n";
-        }
-    }
-    if(type!="origProcId")*fileStream << ")\n";
-    delete fileStream;
-}
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

 } // End namespace Foam
--- a/src/lagrangian/cfdemParticle/subModels/IOModel/trackIO/trackIO.H
+++ b/src/lagrangian/cfdemParticle/subModels/IOModel/trackIO/trackIO.H
@ -39,7 +39,7 @@ SourceFiles
 #ifndef trackIO_H
 #define trackIO_H

-#include "IOModel.H"
+#include "sophIO.H"
 #include "interpolationCellPoint.H"
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

@ -52,24 +52,12 @@ namespace Foam

 class trackIO
 :
-    public IOModel
+    public sophIO
 {

 private:

-    //dictionary propsDict_;
-
-    word dirName_;
-
-    fileName path_;
-
    // Member Functions
-    fileName buildFilePath(word) const;
-
-protected:
-    mutable fileName lagPath_;
-
-    void streamDataToPath(fileName,double**,int,word,word type,word className,word finaliser) const;

 public:

--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/error.d
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/error.d
@ -1,54 +1,54 @@
 error.o: error.cpp /usr/lib/openmpi/include/mpi.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stddef.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stddef.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/mpicxx.h \
- /usr/include/c++/4.4/map /usr/include/c++/4.4/bits/stl_tree.h \
- /usr/include/c++/4.4/bits/stl_algobase.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++config.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/os_defines.h \
+ /usr/include/c++/4.6/map /usr/include/c++/4.6/bits/stl_tree.h \
+ /usr/include/c++/4.6/bits/stl_algobase.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++config.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/os_defines.h \
 /usr/include/features.h /usr/include/bits/predefs.h \
 /usr/include/sys/cdefs.h /usr/include/bits/wordsize.h \
 /usr/include/gnu/stubs.h /usr/include/gnu/stubs-64.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/cpu_defines.h \
- /usr/include/c++/4.4/cstddef /usr/include/c++/4.4/bits/functexcept.h \
- /usr/include/c++/4.4/exception_defines.h \
- /usr/include/c++/4.4/bits/cpp_type_traits.h \
- /usr/include/c++/4.4/ext/type_traits.h \
- /usr/include/c++/4.4/ext/numeric_traits.h \
- /usr/include/c++/4.4/bits/stl_pair.h /usr/include/c++/4.4/bits/move.h \
- /usr/include/c++/4.4/bits/concept_check.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_types.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_funcs.h \
- /usr/include/c++/4.4/bits/stl_iterator.h \
- /usr/include/c++/4.4/debug/debug.h /usr/include/c++/4.4/bits/allocator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++allocator.h \
- /usr/include/c++/4.4/ext/new_allocator.h /usr/include/c++/4.4/new \
- /usr/include/c++/4.4/exception /usr/include/c++/4.4/bits/stl_function.h \
- /usr/include/c++/4.4/backward/binders.h \
- /usr/include/c++/4.4/bits/stl_map.h \
- /usr/include/c++/4.4/initializer_list \
- /usr/include/c++/4.4/bits/stl_multimap.h /usr/include/c++/4.4/utility \
- /usr/include/c++/4.4/bits/stl_relops.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stdarg.h \
- /usr/include/stdio.h /usr/include/bits/types.h \
- /usr/include/bits/typesizes.h /usr/include/libio.h \
- /usr/include/_G_config.h /usr/include/wchar.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/cpu_defines.h \
+ /usr/include/c++/4.6/bits/functexcept.h \
+ /usr/include/c++/4.6/bits/exception_defines.h \
+ /usr/include/c++/4.6/bits/cpp_type_traits.h \
+ /usr/include/c++/4.6/ext/type_traits.h \
+ /usr/include/c++/4.6/ext/numeric_traits.h \
+ /usr/include/c++/4.6/bits/stl_pair.h /usr/include/c++/4.6/bits/move.h \
+ /usr/include/c++/4.6/bits/concept_check.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_types.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_funcs.h \
+ /usr/include/c++/4.6/bits/stl_iterator.h \
+ /usr/include/c++/4.6/debug/debug.h /usr/include/c++/4.6/bits/allocator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++allocator.h \
+ /usr/include/c++/4.6/ext/new_allocator.h /usr/include/c++/4.6/new \
+ /usr/include/c++/4.6/exception /usr/include/c++/4.6/bits/stl_function.h \
+ /usr/include/c++/4.6/backward/binders.h \
+ /usr/include/c++/4.6/bits/stl_map.h \
+ /usr/include/c++/4.6/initializer_list \
+ /usr/include/c++/4.6/bits/stl_multimap.h \
+ /usr/include/c++/4.6/bits/range_access.h /usr/include/c++/4.6/utility \
+ /usr/include/c++/4.6/bits/stl_relops.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stdarg.h /usr/include/stdio.h \
+ /usr/include/bits/types.h /usr/include/bits/typesizes.h \
+ /usr/include/libio.h /usr/include/_G_config.h /usr/include/wchar.h \
 /usr/include/bits/stdio_lim.h /usr/include/bits/sys_errlist.h \
 /usr/include/bits/stdio.h /usr/include/bits/stdio2.h \
- /usr/include/c++/4.4/iostream /usr/include/c++/4.4/ostream \
- /usr/include/c++/4.4/ios /usr/include/c++/4.4/iosfwd \
- /usr/include/c++/4.4/bits/stringfwd.h \
- /usr/include/c++/4.4/bits/postypes.h /usr/include/c++/4.4/cwchar \
+ /usr/include/c++/4.6/iostream /usr/include/c++/4.6/ostream \
+ /usr/include/c++/4.6/ios /usr/include/c++/4.6/iosfwd \
+ /usr/include/c++/4.6/bits/stringfwd.h \
+ /usr/include/c++/4.6/bits/postypes.h /usr/include/c++/4.6/cwchar \
 /usr/include/bits/wchar.h /usr/include/xlocale.h \
- /usr/include/bits/wchar2.h /usr/include/c++/4.4/bits/char_traits.h \
- /usr/include/c++/4.4/bits/localefwd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++locale.h \
- /usr/include/c++/4.4/clocale /usr/include/locale.h \
- /usr/include/bits/locale.h /usr/include/c++/4.4/cctype \
+ /usr/include/bits/wchar2.h /usr/include/c++/4.6/bits/char_traits.h \
+ /usr/include/c++/4.6/bits/localefwd.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++locale.h \
+ /usr/include/c++/4.6/clocale /usr/include/locale.h \
+ /usr/include/bits/locale.h /usr/include/c++/4.6/cctype \
 /usr/include/ctype.h /usr/include/endian.h /usr/include/bits/endian.h \
- /usr/include/bits/byteswap.h /usr/include/c++/4.4/bits/ios_base.h \
- /usr/include/c++/4.4/ext/atomicity.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr-default.h \
+ /usr/include/bits/byteswap.h /usr/include/c++/4.6/bits/ios_base.h \
+ /usr/include/c++/4.6/ext/atomicity.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr-default.h \
 /usr/include/pthread.h /usr/include/sched.h /usr/include/time.h \
 /usr/include/bits/sched.h /usr/include/bits/time.h /usr/include/signal.h \
 /usr/include/bits/sigset.h /usr/include/bits/pthreadtypes.h \
@ -56,24 +56,24 @@ error.o: error.cpp /usr/lib/openmpi/include/mpi.h \
 /usr/include/bits/posix_opt.h /usr/include/bits/environments.h \
 /usr/include/bits/confname.h /usr/include/getopt.h \
 /usr/include/bits/unistd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/atomic_word.h \
- /usr/include/c++/4.4/bits/locale_classes.h /usr/include/c++/4.4/string \
- /usr/include/c++/4.4/bits/ostream_insert.h \
- /usr/include/c++/4.4/cxxabi-forced.h \
- /usr/include/c++/4.4/bits/basic_string.h \
- /usr/include/c++/4.4/bits/basic_string.tcc \
- /usr/include/c++/4.4/bits/locale_classes.tcc \
- /usr/include/c++/4.4/streambuf /usr/include/c++/4.4/bits/streambuf.tcc \
- /usr/include/c++/4.4/bits/basic_ios.h \
- /usr/include/c++/4.4/bits/locale_facets.h /usr/include/c++/4.4/cwctype \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/atomic_word.h \
+ /usr/include/c++/4.6/bits/locale_classes.h /usr/include/c++/4.6/string \
+ /usr/include/c++/4.6/bits/ostream_insert.h \
+ /usr/include/c++/4.6/bits/cxxabi_forced.h \
+ /usr/include/c++/4.6/bits/basic_string.h \
+ /usr/include/c++/4.6/bits/basic_string.tcc \
+ /usr/include/c++/4.6/bits/locale_classes.tcc \
+ /usr/include/c++/4.6/streambuf /usr/include/c++/4.6/bits/streambuf.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.h \
+ /usr/include/c++/4.6/bits/locale_facets.h /usr/include/c++/4.6/cwctype \
 /usr/include/wctype.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_base.h \
- /usr/include/c++/4.4/bits/streambuf_iterator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_inline.h \
- /usr/include/c++/4.4/bits/locale_facets.tcc \
- /usr/include/c++/4.4/bits/basic_ios.tcc \
- /usr/include/c++/4.4/bits/ostream.tcc /usr/include/c++/4.4/istream \
- /usr/include/c++/4.4/bits/istream.tcc \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_base.h \
+ /usr/include/c++/4.6/bits/streambuf_iterator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_inline.h \
+ /usr/include/c++/4.6/bits/locale_facets.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.tcc \
+ /usr/include/c++/4.6/bits/ostream.tcc /usr/include/c++/4.6/istream \
+ /usr/include/c++/4.6/bits/istream.tcc \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/constants.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/functions.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/datatype.h \
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/error.o
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/error.o
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/files.d
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/files.d
@ -2,10 +2,10 @@ files.o: files.cpp /usr/include/stdio.h /usr/include/features.h \
 /usr/include/bits/predefs.h /usr/include/sys/cdefs.h \
 /usr/include/bits/wordsize.h /usr/include/gnu/stubs.h \
 /usr/include/gnu/stubs-64.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stddef.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stddef.h \
 /usr/include/bits/types.h /usr/include/bits/typesizes.h \
 /usr/include/libio.h /usr/include/_G_config.h /usr/include/wchar.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stdarg.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stdarg.h \
 /usr/include/bits/stdio_lim.h /usr/include/bits/sys_errlist.h \
 /usr/include/bits/stdio.h /usr/include/bits/stdio2.h \
 /usr/include/string.h /usr/include/xlocale.h /usr/include/bits/string3.h \
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/files.o
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/files.o
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/irregular.d
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/irregular.d
@ -2,10 +2,10 @@ irregular.o: irregular.cpp /usr/include/stdio.h /usr/include/features.h \
 /usr/include/bits/predefs.h /usr/include/sys/cdefs.h \
 /usr/include/bits/wordsize.h /usr/include/gnu/stubs.h \
 /usr/include/gnu/stubs-64.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stddef.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stddef.h \
 /usr/include/bits/types.h /usr/include/bits/typesizes.h \
 /usr/include/libio.h /usr/include/_G_config.h /usr/include/wchar.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stdarg.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stdarg.h \
 /usr/include/bits/stdio_lim.h /usr/include/bits/sys_errlist.h \
 /usr/include/bits/stdio.h /usr/include/bits/stdio2.h \
 /usr/include/stdlib.h /usr/include/bits/waitflags.h \
@ -18,66 +18,67 @@ irregular.o: irregular.cpp /usr/include/stdio.h /usr/include/features.h \
 /usr/include/alloca.h /usr/include/bits/stdlib.h /usr/include/string.h \
 /usr/include/bits/string3.h irregular.h /usr/lib/openmpi/include/mpi.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/mpicxx.h \
- /usr/include/c++/4.4/map /usr/include/c++/4.4/bits/stl_tree.h \
- /usr/include/c++/4.4/bits/stl_algobase.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++config.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/os_defines.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/cpu_defines.h \
- /usr/include/c++/4.4/cstddef /usr/include/c++/4.4/bits/functexcept.h \
- /usr/include/c++/4.4/exception_defines.h \
- /usr/include/c++/4.4/bits/cpp_type_traits.h \
- /usr/include/c++/4.4/ext/type_traits.h \
- /usr/include/c++/4.4/ext/numeric_traits.h \
- /usr/include/c++/4.4/bits/stl_pair.h /usr/include/c++/4.4/bits/move.h \
- /usr/include/c++/4.4/bits/concept_check.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_types.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_funcs.h \
- /usr/include/c++/4.4/bits/stl_iterator.h \
- /usr/include/c++/4.4/debug/debug.h /usr/include/c++/4.4/bits/allocator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++allocator.h \
- /usr/include/c++/4.4/ext/new_allocator.h /usr/include/c++/4.4/new \
- /usr/include/c++/4.4/exception /usr/include/c++/4.4/bits/stl_function.h \
- /usr/include/c++/4.4/backward/binders.h \
- /usr/include/c++/4.4/bits/stl_map.h \
- /usr/include/c++/4.4/initializer_list \
- /usr/include/c++/4.4/bits/stl_multimap.h /usr/include/c++/4.4/utility \
- /usr/include/c++/4.4/bits/stl_relops.h /usr/include/c++/4.4/iostream \
- /usr/include/c++/4.4/ostream /usr/include/c++/4.4/ios \
- /usr/include/c++/4.4/iosfwd /usr/include/c++/4.4/bits/stringfwd.h \
- /usr/include/c++/4.4/bits/postypes.h /usr/include/c++/4.4/cwchar \
+ /usr/include/c++/4.6/map /usr/include/c++/4.6/bits/stl_tree.h \
+ /usr/include/c++/4.6/bits/stl_algobase.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++config.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/os_defines.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/cpu_defines.h \
+ /usr/include/c++/4.6/bits/functexcept.h \
+ /usr/include/c++/4.6/bits/exception_defines.h \
+ /usr/include/c++/4.6/bits/cpp_type_traits.h \
+ /usr/include/c++/4.6/ext/type_traits.h \
+ /usr/include/c++/4.6/ext/numeric_traits.h \
+ /usr/include/c++/4.6/bits/stl_pair.h /usr/include/c++/4.6/bits/move.h \
+ /usr/include/c++/4.6/bits/concept_check.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_types.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_funcs.h \
+ /usr/include/c++/4.6/bits/stl_iterator.h \
+ /usr/include/c++/4.6/debug/debug.h /usr/include/c++/4.6/bits/allocator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++allocator.h \
+ /usr/include/c++/4.6/ext/new_allocator.h /usr/include/c++/4.6/new \
+ /usr/include/c++/4.6/exception /usr/include/c++/4.6/bits/stl_function.h \
+ /usr/include/c++/4.6/backward/binders.h \
+ /usr/include/c++/4.6/bits/stl_map.h \
+ /usr/include/c++/4.6/initializer_list \
+ /usr/include/c++/4.6/bits/stl_multimap.h \
+ /usr/include/c++/4.6/bits/range_access.h /usr/include/c++/4.6/utility \
+ /usr/include/c++/4.6/bits/stl_relops.h /usr/include/c++/4.6/iostream \
+ /usr/include/c++/4.6/ostream /usr/include/c++/4.6/ios \
+ /usr/include/c++/4.6/iosfwd /usr/include/c++/4.6/bits/stringfwd.h \
+ /usr/include/c++/4.6/bits/postypes.h /usr/include/c++/4.6/cwchar \
 /usr/include/bits/wchar.h /usr/include/bits/wchar2.h \
- /usr/include/c++/4.4/bits/char_traits.h \
- /usr/include/c++/4.4/bits/localefwd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++locale.h \
- /usr/include/c++/4.4/clocale /usr/include/locale.h \
- /usr/include/bits/locale.h /usr/include/c++/4.4/cctype \
- /usr/include/ctype.h /usr/include/c++/4.4/bits/ios_base.h \
- /usr/include/c++/4.4/ext/atomicity.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr-default.h \
+ /usr/include/c++/4.6/bits/char_traits.h \
+ /usr/include/c++/4.6/bits/localefwd.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++locale.h \
+ /usr/include/c++/4.6/clocale /usr/include/locale.h \
+ /usr/include/bits/locale.h /usr/include/c++/4.6/cctype \
+ /usr/include/ctype.h /usr/include/c++/4.6/bits/ios_base.h \
+ /usr/include/c++/4.6/ext/atomicity.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr-default.h \
 /usr/include/pthread.h /usr/include/sched.h /usr/include/bits/sched.h \
 /usr/include/signal.h /usr/include/bits/setjmp.h /usr/include/unistd.h \
 /usr/include/bits/posix_opt.h /usr/include/bits/environments.h \
 /usr/include/bits/confname.h /usr/include/getopt.h \
 /usr/include/bits/unistd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/atomic_word.h \
- /usr/include/c++/4.4/bits/locale_classes.h /usr/include/c++/4.4/string \
- /usr/include/c++/4.4/bits/ostream_insert.h \
- /usr/include/c++/4.4/cxxabi-forced.h \
- /usr/include/c++/4.4/bits/basic_string.h \
- /usr/include/c++/4.4/bits/basic_string.tcc \
- /usr/include/c++/4.4/bits/locale_classes.tcc \
- /usr/include/c++/4.4/streambuf /usr/include/c++/4.4/bits/streambuf.tcc \
- /usr/include/c++/4.4/bits/basic_ios.h \
- /usr/include/c++/4.4/bits/locale_facets.h /usr/include/c++/4.4/cwctype \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/atomic_word.h \
+ /usr/include/c++/4.6/bits/locale_classes.h /usr/include/c++/4.6/string \
+ /usr/include/c++/4.6/bits/ostream_insert.h \
+ /usr/include/c++/4.6/bits/cxxabi_forced.h \
+ /usr/include/c++/4.6/bits/basic_string.h \
+ /usr/include/c++/4.6/bits/basic_string.tcc \
+ /usr/include/c++/4.6/bits/locale_classes.tcc \
+ /usr/include/c++/4.6/streambuf /usr/include/c++/4.6/bits/streambuf.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.h \
+ /usr/include/c++/4.6/bits/locale_facets.h /usr/include/c++/4.6/cwctype \
 /usr/include/wctype.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_base.h \
- /usr/include/c++/4.4/bits/streambuf_iterator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_inline.h \
- /usr/include/c++/4.4/bits/locale_facets.tcc \
- /usr/include/c++/4.4/bits/basic_ios.tcc \
- /usr/include/c++/4.4/bits/ostream.tcc /usr/include/c++/4.4/istream \
- /usr/include/c++/4.4/bits/istream.tcc \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_base.h \
+ /usr/include/c++/4.6/bits/streambuf_iterator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_inline.h \
+ /usr/include/c++/4.6/bits/locale_facets.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.tcc \
+ /usr/include/c++/4.6/bits/ostream.tcc /usr/include/c++/4.6/istream \
+ /usr/include/c++/4.6/bits/istream.tcc \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/constants.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/functions.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/datatype.h \
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/irregular.o
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/irregular.o
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/lammps_data_write.d
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/lammps_data_write.d
@ -2,7 +2,7 @@ lammps_data_write.o: lammps_data_write.cpp /usr/include/stdlib.h \
 /usr/include/features.h /usr/include/bits/predefs.h \
 /usr/include/sys/cdefs.h /usr/include/bits/wordsize.h \
 /usr/include/gnu/stubs.h /usr/include/gnu/stubs-64.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stddef.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stddef.h \
 /usr/include/bits/waitflags.h /usr/include/bits/waitstatus.h \
 /usr/include/endian.h /usr/include/bits/endian.h \
 /usr/include/bits/byteswap.h /usr/include/xlocale.h \
@ -15,71 +15,72 @@ lammps_data_write.o: lammps_data_write.cpp /usr/include/stdlib.h \
 /usr/include/bits/string3.h lammps_data_write.h send2one.h \
 /usr/lib/openmpi/include/mpi.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/mpicxx.h \
- /usr/include/c++/4.4/map /usr/include/c++/4.4/bits/stl_tree.h \
- /usr/include/c++/4.4/bits/stl_algobase.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++config.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/os_defines.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/cpu_defines.h \
- /usr/include/c++/4.4/cstddef /usr/include/c++/4.4/bits/functexcept.h \
- /usr/include/c++/4.4/exception_defines.h \
- /usr/include/c++/4.4/bits/cpp_type_traits.h \
- /usr/include/c++/4.4/ext/type_traits.h \
- /usr/include/c++/4.4/ext/numeric_traits.h \
- /usr/include/c++/4.4/bits/stl_pair.h /usr/include/c++/4.4/bits/move.h \
- /usr/include/c++/4.4/bits/concept_check.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_types.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_funcs.h \
- /usr/include/c++/4.4/bits/stl_iterator.h \
- /usr/include/c++/4.4/debug/debug.h /usr/include/c++/4.4/bits/allocator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++allocator.h \
- /usr/include/c++/4.4/ext/new_allocator.h /usr/include/c++/4.4/new \
- /usr/include/c++/4.4/exception /usr/include/c++/4.4/bits/stl_function.h \
- /usr/include/c++/4.4/backward/binders.h \
- /usr/include/c++/4.4/bits/stl_map.h \
- /usr/include/c++/4.4/initializer_list \
- /usr/include/c++/4.4/bits/stl_multimap.h /usr/include/c++/4.4/utility \
- /usr/include/c++/4.4/bits/stl_relops.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stdarg.h \
- /usr/include/stdio.h /usr/include/libio.h /usr/include/_G_config.h \
- /usr/include/wchar.h /usr/include/bits/stdio_lim.h \
- /usr/include/bits/sys_errlist.h /usr/include/bits/stdio.h \
- /usr/include/bits/stdio2.h /usr/include/c++/4.4/iostream \
- /usr/include/c++/4.4/ostream /usr/include/c++/4.4/ios \
- /usr/include/c++/4.4/iosfwd /usr/include/c++/4.4/bits/stringfwd.h \
- /usr/include/c++/4.4/bits/postypes.h /usr/include/c++/4.4/cwchar \
+ /usr/include/c++/4.6/map /usr/include/c++/4.6/bits/stl_tree.h \
+ /usr/include/c++/4.6/bits/stl_algobase.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++config.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/os_defines.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/cpu_defines.h \
+ /usr/include/c++/4.6/bits/functexcept.h \
+ /usr/include/c++/4.6/bits/exception_defines.h \
+ /usr/include/c++/4.6/bits/cpp_type_traits.h \
+ /usr/include/c++/4.6/ext/type_traits.h \
+ /usr/include/c++/4.6/ext/numeric_traits.h \
+ /usr/include/c++/4.6/bits/stl_pair.h /usr/include/c++/4.6/bits/move.h \
+ /usr/include/c++/4.6/bits/concept_check.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_types.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_funcs.h \
+ /usr/include/c++/4.6/bits/stl_iterator.h \
+ /usr/include/c++/4.6/debug/debug.h /usr/include/c++/4.6/bits/allocator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++allocator.h \
+ /usr/include/c++/4.6/ext/new_allocator.h /usr/include/c++/4.6/new \
+ /usr/include/c++/4.6/exception /usr/include/c++/4.6/bits/stl_function.h \
+ /usr/include/c++/4.6/backward/binders.h \
+ /usr/include/c++/4.6/bits/stl_map.h \
+ /usr/include/c++/4.6/initializer_list \
+ /usr/include/c++/4.6/bits/stl_multimap.h \
+ /usr/include/c++/4.6/bits/range_access.h /usr/include/c++/4.6/utility \
+ /usr/include/c++/4.6/bits/stl_relops.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stdarg.h /usr/include/stdio.h \
+ /usr/include/libio.h /usr/include/_G_config.h /usr/include/wchar.h \
+ /usr/include/bits/stdio_lim.h /usr/include/bits/sys_errlist.h \
+ /usr/include/bits/stdio.h /usr/include/bits/stdio2.h \
+ /usr/include/c++/4.6/iostream /usr/include/c++/4.6/ostream \
+ /usr/include/c++/4.6/ios /usr/include/c++/4.6/iosfwd \
+ /usr/include/c++/4.6/bits/stringfwd.h \
+ /usr/include/c++/4.6/bits/postypes.h /usr/include/c++/4.6/cwchar \
 /usr/include/bits/wchar.h /usr/include/bits/wchar2.h \
- /usr/include/c++/4.4/bits/char_traits.h \
- /usr/include/c++/4.4/bits/localefwd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++locale.h \
- /usr/include/c++/4.4/clocale /usr/include/locale.h \
- /usr/include/bits/locale.h /usr/include/c++/4.4/cctype \
- /usr/include/ctype.h /usr/include/c++/4.4/bits/ios_base.h \
- /usr/include/c++/4.4/ext/atomicity.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr-default.h \
+ /usr/include/c++/4.6/bits/char_traits.h \
+ /usr/include/c++/4.6/bits/localefwd.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++locale.h \
+ /usr/include/c++/4.6/clocale /usr/include/locale.h \
+ /usr/include/bits/locale.h /usr/include/c++/4.6/cctype \
+ /usr/include/ctype.h /usr/include/c++/4.6/bits/ios_base.h \
+ /usr/include/c++/4.6/ext/atomicity.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr-default.h \
 /usr/include/pthread.h /usr/include/sched.h /usr/include/bits/sched.h \
 /usr/include/signal.h /usr/include/bits/setjmp.h /usr/include/unistd.h \
 /usr/include/bits/posix_opt.h /usr/include/bits/environments.h \
 /usr/include/bits/confname.h /usr/include/getopt.h \
 /usr/include/bits/unistd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/atomic_word.h \
- /usr/include/c++/4.4/bits/locale_classes.h /usr/include/c++/4.4/string \
- /usr/include/c++/4.4/bits/ostream_insert.h \
- /usr/include/c++/4.4/cxxabi-forced.h \
- /usr/include/c++/4.4/bits/basic_string.h \
- /usr/include/c++/4.4/bits/basic_string.tcc \
- /usr/include/c++/4.4/bits/locale_classes.tcc \
- /usr/include/c++/4.4/streambuf /usr/include/c++/4.4/bits/streambuf.tcc \
- /usr/include/c++/4.4/bits/basic_ios.h \
- /usr/include/c++/4.4/bits/locale_facets.h /usr/include/c++/4.4/cwctype \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/atomic_word.h \
+ /usr/include/c++/4.6/bits/locale_classes.h /usr/include/c++/4.6/string \
+ /usr/include/c++/4.6/bits/ostream_insert.h \
+ /usr/include/c++/4.6/bits/cxxabi_forced.h \
+ /usr/include/c++/4.6/bits/basic_string.h \
+ /usr/include/c++/4.6/bits/basic_string.tcc \
+ /usr/include/c++/4.6/bits/locale_classes.tcc \
+ /usr/include/c++/4.6/streambuf /usr/include/c++/4.6/bits/streambuf.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.h \
+ /usr/include/c++/4.6/bits/locale_facets.h /usr/include/c++/4.6/cwctype \
 /usr/include/wctype.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_base.h \
- /usr/include/c++/4.4/bits/streambuf_iterator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_inline.h \
- /usr/include/c++/4.4/bits/locale_facets.tcc \
- /usr/include/c++/4.4/bits/basic_ios.tcc \
- /usr/include/c++/4.4/bits/ostream.tcc /usr/include/c++/4.4/istream \
- /usr/include/c++/4.4/bits/istream.tcc \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_base.h \
+ /usr/include/c++/4.6/bits/streambuf_iterator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_inline.h \
+ /usr/include/c++/4.6/bits/locale_facets.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.tcc \
+ /usr/include/c++/4.6/bits/ostream.tcc /usr/include/c++/4.6/istream \
+ /usr/include/c++/4.6/bits/istream.tcc \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/constants.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/functions.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/datatype.h \
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/lammps_data_write.o
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/lammps_data_write.o
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/libcouple.a
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/libcouple.a
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/many2many.d
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/many2many.d
@ -1,54 +1,54 @@
 many2many.o: many2many.cpp /usr/lib/openmpi/include/mpi.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stddef.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stddef.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/mpicxx.h \
- /usr/include/c++/4.4/map /usr/include/c++/4.4/bits/stl_tree.h \
- /usr/include/c++/4.4/bits/stl_algobase.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++config.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/os_defines.h \
+ /usr/include/c++/4.6/map /usr/include/c++/4.6/bits/stl_tree.h \
+ /usr/include/c++/4.6/bits/stl_algobase.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++config.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/os_defines.h \
 /usr/include/features.h /usr/include/bits/predefs.h \
 /usr/include/sys/cdefs.h /usr/include/bits/wordsize.h \
 /usr/include/gnu/stubs.h /usr/include/gnu/stubs-64.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/cpu_defines.h \
- /usr/include/c++/4.4/cstddef /usr/include/c++/4.4/bits/functexcept.h \
- /usr/include/c++/4.4/exception_defines.h \
- /usr/include/c++/4.4/bits/cpp_type_traits.h \
- /usr/include/c++/4.4/ext/type_traits.h \
- /usr/include/c++/4.4/ext/numeric_traits.h \
- /usr/include/c++/4.4/bits/stl_pair.h /usr/include/c++/4.4/bits/move.h \
- /usr/include/c++/4.4/bits/concept_check.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_types.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_funcs.h \
- /usr/include/c++/4.4/bits/stl_iterator.h \
- /usr/include/c++/4.4/debug/debug.h /usr/include/c++/4.4/bits/allocator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++allocator.h \
- /usr/include/c++/4.4/ext/new_allocator.h /usr/include/c++/4.4/new \
- /usr/include/c++/4.4/exception /usr/include/c++/4.4/bits/stl_function.h \
- /usr/include/c++/4.4/backward/binders.h \
- /usr/include/c++/4.4/bits/stl_map.h \
- /usr/include/c++/4.4/initializer_list \
- /usr/include/c++/4.4/bits/stl_multimap.h /usr/include/c++/4.4/utility \
- /usr/include/c++/4.4/bits/stl_relops.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stdarg.h \
- /usr/include/stdio.h /usr/include/bits/types.h \
- /usr/include/bits/typesizes.h /usr/include/libio.h \
- /usr/include/_G_config.h /usr/include/wchar.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/cpu_defines.h \
+ /usr/include/c++/4.6/bits/functexcept.h \
+ /usr/include/c++/4.6/bits/exception_defines.h \
+ /usr/include/c++/4.6/bits/cpp_type_traits.h \
+ /usr/include/c++/4.6/ext/type_traits.h \
+ /usr/include/c++/4.6/ext/numeric_traits.h \
+ /usr/include/c++/4.6/bits/stl_pair.h /usr/include/c++/4.6/bits/move.h \
+ /usr/include/c++/4.6/bits/concept_check.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_types.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_funcs.h \
+ /usr/include/c++/4.6/bits/stl_iterator.h \
+ /usr/include/c++/4.6/debug/debug.h /usr/include/c++/4.6/bits/allocator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++allocator.h \
+ /usr/include/c++/4.6/ext/new_allocator.h /usr/include/c++/4.6/new \
+ /usr/include/c++/4.6/exception /usr/include/c++/4.6/bits/stl_function.h \
+ /usr/include/c++/4.6/backward/binders.h \
+ /usr/include/c++/4.6/bits/stl_map.h \
+ /usr/include/c++/4.6/initializer_list \
+ /usr/include/c++/4.6/bits/stl_multimap.h \
+ /usr/include/c++/4.6/bits/range_access.h /usr/include/c++/4.6/utility \
+ /usr/include/c++/4.6/bits/stl_relops.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stdarg.h /usr/include/stdio.h \
+ /usr/include/bits/types.h /usr/include/bits/typesizes.h \
+ /usr/include/libio.h /usr/include/_G_config.h /usr/include/wchar.h \
 /usr/include/bits/stdio_lim.h /usr/include/bits/sys_errlist.h \
 /usr/include/bits/stdio.h /usr/include/bits/stdio2.h \
- /usr/include/c++/4.4/iostream /usr/include/c++/4.4/ostream \
- /usr/include/c++/4.4/ios /usr/include/c++/4.4/iosfwd \
- /usr/include/c++/4.4/bits/stringfwd.h \
- /usr/include/c++/4.4/bits/postypes.h /usr/include/c++/4.4/cwchar \
+ /usr/include/c++/4.6/iostream /usr/include/c++/4.6/ostream \
+ /usr/include/c++/4.6/ios /usr/include/c++/4.6/iosfwd \
+ /usr/include/c++/4.6/bits/stringfwd.h \
+ /usr/include/c++/4.6/bits/postypes.h /usr/include/c++/4.6/cwchar \
 /usr/include/bits/wchar.h /usr/include/xlocale.h \
- /usr/include/bits/wchar2.h /usr/include/c++/4.4/bits/char_traits.h \
- /usr/include/c++/4.4/bits/localefwd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++locale.h \
- /usr/include/c++/4.4/clocale /usr/include/locale.h \
- /usr/include/bits/locale.h /usr/include/c++/4.4/cctype \
+ /usr/include/bits/wchar2.h /usr/include/c++/4.6/bits/char_traits.h \
+ /usr/include/c++/4.6/bits/localefwd.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++locale.h \
+ /usr/include/c++/4.6/clocale /usr/include/locale.h \
+ /usr/include/bits/locale.h /usr/include/c++/4.6/cctype \
 /usr/include/ctype.h /usr/include/endian.h /usr/include/bits/endian.h \
- /usr/include/bits/byteswap.h /usr/include/c++/4.4/bits/ios_base.h \
- /usr/include/c++/4.4/ext/atomicity.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr-default.h \
+ /usr/include/bits/byteswap.h /usr/include/c++/4.6/bits/ios_base.h \
+ /usr/include/c++/4.6/ext/atomicity.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr-default.h \
 /usr/include/pthread.h /usr/include/sched.h /usr/include/time.h \
 /usr/include/bits/sched.h /usr/include/bits/time.h /usr/include/signal.h \
 /usr/include/bits/sigset.h /usr/include/bits/pthreadtypes.h \
@ -56,24 +56,24 @@ many2many.o: many2many.cpp /usr/lib/openmpi/include/mpi.h \
 /usr/include/bits/posix_opt.h /usr/include/bits/environments.h \
 /usr/include/bits/confname.h /usr/include/getopt.h \
 /usr/include/bits/unistd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/atomic_word.h \
- /usr/include/c++/4.4/bits/locale_classes.h /usr/include/c++/4.4/string \
- /usr/include/c++/4.4/bits/ostream_insert.h \
- /usr/include/c++/4.4/cxxabi-forced.h \
- /usr/include/c++/4.4/bits/basic_string.h \
- /usr/include/c++/4.4/bits/basic_string.tcc \
- /usr/include/c++/4.4/bits/locale_classes.tcc \
- /usr/include/c++/4.4/streambuf /usr/include/c++/4.4/bits/streambuf.tcc \
- /usr/include/c++/4.4/bits/basic_ios.h \
- /usr/include/c++/4.4/bits/locale_facets.h /usr/include/c++/4.4/cwctype \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/atomic_word.h \
+ /usr/include/c++/4.6/bits/locale_classes.h /usr/include/c++/4.6/string \
+ /usr/include/c++/4.6/bits/ostream_insert.h \
+ /usr/include/c++/4.6/bits/cxxabi_forced.h \
+ /usr/include/c++/4.6/bits/basic_string.h \
+ /usr/include/c++/4.6/bits/basic_string.tcc \
+ /usr/include/c++/4.6/bits/locale_classes.tcc \
+ /usr/include/c++/4.6/streambuf /usr/include/c++/4.6/bits/streambuf.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.h \
+ /usr/include/c++/4.6/bits/locale_facets.h /usr/include/c++/4.6/cwctype \
 /usr/include/wctype.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_base.h \
- /usr/include/c++/4.4/bits/streambuf_iterator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_inline.h \
- /usr/include/c++/4.4/bits/locale_facets.tcc \
- /usr/include/c++/4.4/bits/basic_ios.tcc \
- /usr/include/c++/4.4/bits/ostream.tcc /usr/include/c++/4.4/istream \
- /usr/include/c++/4.4/bits/istream.tcc \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_base.h \
+ /usr/include/c++/4.6/bits/streambuf_iterator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_inline.h \
+ /usr/include/c++/4.6/bits/locale_facets.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.tcc \
+ /usr/include/c++/4.6/bits/ostream.tcc /usr/include/c++/4.6/istream \
+ /usr/include/c++/4.6/bits/istream.tcc \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/constants.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/functions.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/datatype.h \
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/many2many.o
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/many2many.o
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/many2one.d
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/many2one.d
@ -1,54 +1,54 @@
 many2one.o: many2one.cpp /usr/lib/openmpi/include/mpi.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stddef.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stddef.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/mpicxx.h \
- /usr/include/c++/4.4/map /usr/include/c++/4.4/bits/stl_tree.h \
- /usr/include/c++/4.4/bits/stl_algobase.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++config.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/os_defines.h \
+ /usr/include/c++/4.6/map /usr/include/c++/4.6/bits/stl_tree.h \
+ /usr/include/c++/4.6/bits/stl_algobase.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++config.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/os_defines.h \
 /usr/include/features.h /usr/include/bits/predefs.h \
 /usr/include/sys/cdefs.h /usr/include/bits/wordsize.h \
 /usr/include/gnu/stubs.h /usr/include/gnu/stubs-64.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/cpu_defines.h \
- /usr/include/c++/4.4/cstddef /usr/include/c++/4.4/bits/functexcept.h \
- /usr/include/c++/4.4/exception_defines.h \
- /usr/include/c++/4.4/bits/cpp_type_traits.h \
- /usr/include/c++/4.4/ext/type_traits.h \
- /usr/include/c++/4.4/ext/numeric_traits.h \
- /usr/include/c++/4.4/bits/stl_pair.h /usr/include/c++/4.4/bits/move.h \
- /usr/include/c++/4.4/bits/concept_check.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_types.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_funcs.h \
- /usr/include/c++/4.4/bits/stl_iterator.h \
- /usr/include/c++/4.4/debug/debug.h /usr/include/c++/4.4/bits/allocator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++allocator.h \
- /usr/include/c++/4.4/ext/new_allocator.h /usr/include/c++/4.4/new \
- /usr/include/c++/4.4/exception /usr/include/c++/4.4/bits/stl_function.h \
- /usr/include/c++/4.4/backward/binders.h \
- /usr/include/c++/4.4/bits/stl_map.h \
- /usr/include/c++/4.4/initializer_list \
- /usr/include/c++/4.4/bits/stl_multimap.h /usr/include/c++/4.4/utility \
- /usr/include/c++/4.4/bits/stl_relops.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stdarg.h \
- /usr/include/stdio.h /usr/include/bits/types.h \
- /usr/include/bits/typesizes.h /usr/include/libio.h \
- /usr/include/_G_config.h /usr/include/wchar.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/cpu_defines.h \
+ /usr/include/c++/4.6/bits/functexcept.h \
+ /usr/include/c++/4.6/bits/exception_defines.h \
+ /usr/include/c++/4.6/bits/cpp_type_traits.h \
+ /usr/include/c++/4.6/ext/type_traits.h \
+ /usr/include/c++/4.6/ext/numeric_traits.h \
+ /usr/include/c++/4.6/bits/stl_pair.h /usr/include/c++/4.6/bits/move.h \
+ /usr/include/c++/4.6/bits/concept_check.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_types.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_funcs.h \
+ /usr/include/c++/4.6/bits/stl_iterator.h \
+ /usr/include/c++/4.6/debug/debug.h /usr/include/c++/4.6/bits/allocator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++allocator.h \
+ /usr/include/c++/4.6/ext/new_allocator.h /usr/include/c++/4.6/new \
+ /usr/include/c++/4.6/exception /usr/include/c++/4.6/bits/stl_function.h \
+ /usr/include/c++/4.6/backward/binders.h \
+ /usr/include/c++/4.6/bits/stl_map.h \
+ /usr/include/c++/4.6/initializer_list \
+ /usr/include/c++/4.6/bits/stl_multimap.h \
+ /usr/include/c++/4.6/bits/range_access.h /usr/include/c++/4.6/utility \
+ /usr/include/c++/4.6/bits/stl_relops.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stdarg.h /usr/include/stdio.h \
+ /usr/include/bits/types.h /usr/include/bits/typesizes.h \
+ /usr/include/libio.h /usr/include/_G_config.h /usr/include/wchar.h \
 /usr/include/bits/stdio_lim.h /usr/include/bits/sys_errlist.h \
 /usr/include/bits/stdio.h /usr/include/bits/stdio2.h \
- /usr/include/c++/4.4/iostream /usr/include/c++/4.4/ostream \
- /usr/include/c++/4.4/ios /usr/include/c++/4.4/iosfwd \
- /usr/include/c++/4.4/bits/stringfwd.h \
- /usr/include/c++/4.4/bits/postypes.h /usr/include/c++/4.4/cwchar \
+ /usr/include/c++/4.6/iostream /usr/include/c++/4.6/ostream \
+ /usr/include/c++/4.6/ios /usr/include/c++/4.6/iosfwd \
+ /usr/include/c++/4.6/bits/stringfwd.h \
+ /usr/include/c++/4.6/bits/postypes.h /usr/include/c++/4.6/cwchar \
 /usr/include/bits/wchar.h /usr/include/xlocale.h \
- /usr/include/bits/wchar2.h /usr/include/c++/4.4/bits/char_traits.h \
- /usr/include/c++/4.4/bits/localefwd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++locale.h \
- /usr/include/c++/4.4/clocale /usr/include/locale.h \
- /usr/include/bits/locale.h /usr/include/c++/4.4/cctype \
+ /usr/include/bits/wchar2.h /usr/include/c++/4.6/bits/char_traits.h \
+ /usr/include/c++/4.6/bits/localefwd.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++locale.h \
+ /usr/include/c++/4.6/clocale /usr/include/locale.h \
+ /usr/include/bits/locale.h /usr/include/c++/4.6/cctype \
 /usr/include/ctype.h /usr/include/endian.h /usr/include/bits/endian.h \
- /usr/include/bits/byteswap.h /usr/include/c++/4.4/bits/ios_base.h \
- /usr/include/c++/4.4/ext/atomicity.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr-default.h \
+ /usr/include/bits/byteswap.h /usr/include/c++/4.6/bits/ios_base.h \
+ /usr/include/c++/4.6/ext/atomicity.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr-default.h \
 /usr/include/pthread.h /usr/include/sched.h /usr/include/time.h \
 /usr/include/bits/sched.h /usr/include/bits/time.h /usr/include/signal.h \
 /usr/include/bits/sigset.h /usr/include/bits/pthreadtypes.h \
@ -56,24 +56,24 @@ many2one.o: many2one.cpp /usr/lib/openmpi/include/mpi.h \
 /usr/include/bits/posix_opt.h /usr/include/bits/environments.h \
 /usr/include/bits/confname.h /usr/include/getopt.h \
 /usr/include/bits/unistd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/atomic_word.h \
- /usr/include/c++/4.4/bits/locale_classes.h /usr/include/c++/4.4/string \
- /usr/include/c++/4.4/bits/ostream_insert.h \
- /usr/include/c++/4.4/cxxabi-forced.h \
- /usr/include/c++/4.4/bits/basic_string.h \
- /usr/include/c++/4.4/bits/basic_string.tcc \
- /usr/include/c++/4.4/bits/locale_classes.tcc \
- /usr/include/c++/4.4/streambuf /usr/include/c++/4.4/bits/streambuf.tcc \
- /usr/include/c++/4.4/bits/basic_ios.h \
- /usr/include/c++/4.4/bits/locale_facets.h /usr/include/c++/4.4/cwctype \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/atomic_word.h \
+ /usr/include/c++/4.6/bits/locale_classes.h /usr/include/c++/4.6/string \
+ /usr/include/c++/4.6/bits/ostream_insert.h \
+ /usr/include/c++/4.6/bits/cxxabi_forced.h \
+ /usr/include/c++/4.6/bits/basic_string.h \
+ /usr/include/c++/4.6/bits/basic_string.tcc \
+ /usr/include/c++/4.6/bits/locale_classes.tcc \
+ /usr/include/c++/4.6/streambuf /usr/include/c++/4.6/bits/streambuf.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.h \
+ /usr/include/c++/4.6/bits/locale_facets.h /usr/include/c++/4.6/cwctype \
 /usr/include/wctype.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_base.h \
- /usr/include/c++/4.4/bits/streambuf_iterator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_inline.h \
- /usr/include/c++/4.4/bits/locale_facets.tcc \
- /usr/include/c++/4.4/bits/basic_ios.tcc \
- /usr/include/c++/4.4/bits/ostream.tcc /usr/include/c++/4.4/istream \
- /usr/include/c++/4.4/bits/istream.tcc \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_base.h \
+ /usr/include/c++/4.6/bits/streambuf_iterator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_inline.h \
+ /usr/include/c++/4.6/bits/locale_facets.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.tcc \
+ /usr/include/c++/4.6/bits/ostream.tcc /usr/include/c++/4.6/istream \
+ /usr/include/c++/4.6/bits/istream.tcc \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/constants.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/functions.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/datatype.h \
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/many2one.o
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/many2one.o
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/memory.d
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/memory.d
@ -1,54 +1,54 @@
 memory.o: memory.cpp /usr/lib/openmpi/include/mpi.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stddef.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stddef.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/mpicxx.h \
- /usr/include/c++/4.4/map /usr/include/c++/4.4/bits/stl_tree.h \
- /usr/include/c++/4.4/bits/stl_algobase.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++config.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/os_defines.h \
+ /usr/include/c++/4.6/map /usr/include/c++/4.6/bits/stl_tree.h \
+ /usr/include/c++/4.6/bits/stl_algobase.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++config.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/os_defines.h \
 /usr/include/features.h /usr/include/bits/predefs.h \
 /usr/include/sys/cdefs.h /usr/include/bits/wordsize.h \
 /usr/include/gnu/stubs.h /usr/include/gnu/stubs-64.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/cpu_defines.h \
- /usr/include/c++/4.4/cstddef /usr/include/c++/4.4/bits/functexcept.h \
- /usr/include/c++/4.4/exception_defines.h \
- /usr/include/c++/4.4/bits/cpp_type_traits.h \
- /usr/include/c++/4.4/ext/type_traits.h \
- /usr/include/c++/4.4/ext/numeric_traits.h \
- /usr/include/c++/4.4/bits/stl_pair.h /usr/include/c++/4.4/bits/move.h \
- /usr/include/c++/4.4/bits/concept_check.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_types.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_funcs.h \
- /usr/include/c++/4.4/bits/stl_iterator.h \
- /usr/include/c++/4.4/debug/debug.h /usr/include/c++/4.4/bits/allocator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++allocator.h \
- /usr/include/c++/4.4/ext/new_allocator.h /usr/include/c++/4.4/new \
- /usr/include/c++/4.4/exception /usr/include/c++/4.4/bits/stl_function.h \
- /usr/include/c++/4.4/backward/binders.h \
- /usr/include/c++/4.4/bits/stl_map.h \
- /usr/include/c++/4.4/initializer_list \
- /usr/include/c++/4.4/bits/stl_multimap.h /usr/include/c++/4.4/utility \
- /usr/include/c++/4.4/bits/stl_relops.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stdarg.h \
- /usr/include/stdio.h /usr/include/bits/types.h \
- /usr/include/bits/typesizes.h /usr/include/libio.h \
- /usr/include/_G_config.h /usr/include/wchar.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/cpu_defines.h \
+ /usr/include/c++/4.6/bits/functexcept.h \
+ /usr/include/c++/4.6/bits/exception_defines.h \
+ /usr/include/c++/4.6/bits/cpp_type_traits.h \
+ /usr/include/c++/4.6/ext/type_traits.h \
+ /usr/include/c++/4.6/ext/numeric_traits.h \
+ /usr/include/c++/4.6/bits/stl_pair.h /usr/include/c++/4.6/bits/move.h \
+ /usr/include/c++/4.6/bits/concept_check.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_types.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_funcs.h \
+ /usr/include/c++/4.6/bits/stl_iterator.h \
+ /usr/include/c++/4.6/debug/debug.h /usr/include/c++/4.6/bits/allocator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++allocator.h \
+ /usr/include/c++/4.6/ext/new_allocator.h /usr/include/c++/4.6/new \
+ /usr/include/c++/4.6/exception /usr/include/c++/4.6/bits/stl_function.h \
+ /usr/include/c++/4.6/backward/binders.h \
+ /usr/include/c++/4.6/bits/stl_map.h \
+ /usr/include/c++/4.6/initializer_list \
+ /usr/include/c++/4.6/bits/stl_multimap.h \
+ /usr/include/c++/4.6/bits/range_access.h /usr/include/c++/4.6/utility \
+ /usr/include/c++/4.6/bits/stl_relops.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stdarg.h /usr/include/stdio.h \
+ /usr/include/bits/types.h /usr/include/bits/typesizes.h \
+ /usr/include/libio.h /usr/include/_G_config.h /usr/include/wchar.h \
 /usr/include/bits/stdio_lim.h /usr/include/bits/sys_errlist.h \
 /usr/include/bits/stdio.h /usr/include/bits/stdio2.h \
- /usr/include/c++/4.4/iostream /usr/include/c++/4.4/ostream \
- /usr/include/c++/4.4/ios /usr/include/c++/4.4/iosfwd \
- /usr/include/c++/4.4/bits/stringfwd.h \
- /usr/include/c++/4.4/bits/postypes.h /usr/include/c++/4.4/cwchar \
+ /usr/include/c++/4.6/iostream /usr/include/c++/4.6/ostream \
+ /usr/include/c++/4.6/ios /usr/include/c++/4.6/iosfwd \
+ /usr/include/c++/4.6/bits/stringfwd.h \
+ /usr/include/c++/4.6/bits/postypes.h /usr/include/c++/4.6/cwchar \
 /usr/include/bits/wchar.h /usr/include/xlocale.h \
- /usr/include/bits/wchar2.h /usr/include/c++/4.4/bits/char_traits.h \
- /usr/include/c++/4.4/bits/localefwd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++locale.h \
- /usr/include/c++/4.4/clocale /usr/include/locale.h \
- /usr/include/bits/locale.h /usr/include/c++/4.4/cctype \
+ /usr/include/bits/wchar2.h /usr/include/c++/4.6/bits/char_traits.h \
+ /usr/include/c++/4.6/bits/localefwd.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++locale.h \
+ /usr/include/c++/4.6/clocale /usr/include/locale.h \
+ /usr/include/bits/locale.h /usr/include/c++/4.6/cctype \
 /usr/include/ctype.h /usr/include/endian.h /usr/include/bits/endian.h \
- /usr/include/bits/byteswap.h /usr/include/c++/4.4/bits/ios_base.h \
- /usr/include/c++/4.4/ext/atomicity.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr-default.h \
+ /usr/include/bits/byteswap.h /usr/include/c++/4.6/bits/ios_base.h \
+ /usr/include/c++/4.6/ext/atomicity.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr-default.h \
 /usr/include/pthread.h /usr/include/sched.h /usr/include/time.h \
 /usr/include/bits/sched.h /usr/include/bits/time.h /usr/include/signal.h \
 /usr/include/bits/sigset.h /usr/include/bits/pthreadtypes.h \
@ -56,24 +56,24 @@ memory.o: memory.cpp /usr/lib/openmpi/include/mpi.h \
 /usr/include/bits/posix_opt.h /usr/include/bits/environments.h \
 /usr/include/bits/confname.h /usr/include/getopt.h \
 /usr/include/bits/unistd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/atomic_word.h \
- /usr/include/c++/4.4/bits/locale_classes.h /usr/include/c++/4.4/string \
- /usr/include/c++/4.4/bits/ostream_insert.h \
- /usr/include/c++/4.4/cxxabi-forced.h \
- /usr/include/c++/4.4/bits/basic_string.h \
- /usr/include/c++/4.4/bits/basic_string.tcc \
- /usr/include/c++/4.4/bits/locale_classes.tcc \
- /usr/include/c++/4.4/streambuf /usr/include/c++/4.4/bits/streambuf.tcc \
- /usr/include/c++/4.4/bits/basic_ios.h \
- /usr/include/c++/4.4/bits/locale_facets.h /usr/include/c++/4.4/cwctype \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/atomic_word.h \
+ /usr/include/c++/4.6/bits/locale_classes.h /usr/include/c++/4.6/string \
+ /usr/include/c++/4.6/bits/ostream_insert.h \
+ /usr/include/c++/4.6/bits/cxxabi_forced.h \
+ /usr/include/c++/4.6/bits/basic_string.h \
+ /usr/include/c++/4.6/bits/basic_string.tcc \
+ /usr/include/c++/4.6/bits/locale_classes.tcc \
+ /usr/include/c++/4.6/streambuf /usr/include/c++/4.6/bits/streambuf.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.h \
+ /usr/include/c++/4.6/bits/locale_facets.h /usr/include/c++/4.6/cwctype \
 /usr/include/wctype.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_base.h \
- /usr/include/c++/4.4/bits/streambuf_iterator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_inline.h \
- /usr/include/c++/4.4/bits/locale_facets.tcc \
- /usr/include/c++/4.4/bits/basic_ios.tcc \
- /usr/include/c++/4.4/bits/ostream.tcc /usr/include/c++/4.4/istream \
- /usr/include/c++/4.4/bits/istream.tcc \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_base.h \
+ /usr/include/c++/4.6/bits/streambuf_iterator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_inline.h \
+ /usr/include/c++/4.6/bits/locale_facets.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.tcc \
+ /usr/include/c++/4.6/bits/ostream.tcc /usr/include/c++/4.6/istream \
+ /usr/include/c++/4.6/bits/istream.tcc \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/constants.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/functions.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/datatype.h \
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/memory.o
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/memory.o
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/one2many.d
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/one2many.d
@ -1,54 +1,54 @@
 one2many.o: one2many.cpp /usr/lib/openmpi/include/mpi.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stddef.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stddef.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/mpicxx.h \
- /usr/include/c++/4.4/map /usr/include/c++/4.4/bits/stl_tree.h \
- /usr/include/c++/4.4/bits/stl_algobase.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++config.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/os_defines.h \
+ /usr/include/c++/4.6/map /usr/include/c++/4.6/bits/stl_tree.h \
+ /usr/include/c++/4.6/bits/stl_algobase.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++config.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/os_defines.h \
 /usr/include/features.h /usr/include/bits/predefs.h \
 /usr/include/sys/cdefs.h /usr/include/bits/wordsize.h \
 /usr/include/gnu/stubs.h /usr/include/gnu/stubs-64.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/cpu_defines.h \
- /usr/include/c++/4.4/cstddef /usr/include/c++/4.4/bits/functexcept.h \
- /usr/include/c++/4.4/exception_defines.h \
- /usr/include/c++/4.4/bits/cpp_type_traits.h \
- /usr/include/c++/4.4/ext/type_traits.h \
- /usr/include/c++/4.4/ext/numeric_traits.h \
- /usr/include/c++/4.4/bits/stl_pair.h /usr/include/c++/4.4/bits/move.h \
- /usr/include/c++/4.4/bits/concept_check.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_types.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_funcs.h \
- /usr/include/c++/4.4/bits/stl_iterator.h \
- /usr/include/c++/4.4/debug/debug.h /usr/include/c++/4.4/bits/allocator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++allocator.h \
- /usr/include/c++/4.4/ext/new_allocator.h /usr/include/c++/4.4/new \
- /usr/include/c++/4.4/exception /usr/include/c++/4.4/bits/stl_function.h \
- /usr/include/c++/4.4/backward/binders.h \
- /usr/include/c++/4.4/bits/stl_map.h \
- /usr/include/c++/4.4/initializer_list \
- /usr/include/c++/4.4/bits/stl_multimap.h /usr/include/c++/4.4/utility \
- /usr/include/c++/4.4/bits/stl_relops.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stdarg.h \
- /usr/include/stdio.h /usr/include/bits/types.h \
- /usr/include/bits/typesizes.h /usr/include/libio.h \
- /usr/include/_G_config.h /usr/include/wchar.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/cpu_defines.h \
+ /usr/include/c++/4.6/bits/functexcept.h \
+ /usr/include/c++/4.6/bits/exception_defines.h \
+ /usr/include/c++/4.6/bits/cpp_type_traits.h \
+ /usr/include/c++/4.6/ext/type_traits.h \
+ /usr/include/c++/4.6/ext/numeric_traits.h \
+ /usr/include/c++/4.6/bits/stl_pair.h /usr/include/c++/4.6/bits/move.h \
+ /usr/include/c++/4.6/bits/concept_check.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_types.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_funcs.h \
+ /usr/include/c++/4.6/bits/stl_iterator.h \
+ /usr/include/c++/4.6/debug/debug.h /usr/include/c++/4.6/bits/allocator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++allocator.h \
+ /usr/include/c++/4.6/ext/new_allocator.h /usr/include/c++/4.6/new \
+ /usr/include/c++/4.6/exception /usr/include/c++/4.6/bits/stl_function.h \
+ /usr/include/c++/4.6/backward/binders.h \
+ /usr/include/c++/4.6/bits/stl_map.h \
+ /usr/include/c++/4.6/initializer_list \
+ /usr/include/c++/4.6/bits/stl_multimap.h \
+ /usr/include/c++/4.6/bits/range_access.h /usr/include/c++/4.6/utility \
+ /usr/include/c++/4.6/bits/stl_relops.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stdarg.h /usr/include/stdio.h \
+ /usr/include/bits/types.h /usr/include/bits/typesizes.h \
+ /usr/include/libio.h /usr/include/_G_config.h /usr/include/wchar.h \
 /usr/include/bits/stdio_lim.h /usr/include/bits/sys_errlist.h \
 /usr/include/bits/stdio.h /usr/include/bits/stdio2.h \
- /usr/include/c++/4.4/iostream /usr/include/c++/4.4/ostream \
- /usr/include/c++/4.4/ios /usr/include/c++/4.4/iosfwd \
- /usr/include/c++/4.4/bits/stringfwd.h \
- /usr/include/c++/4.4/bits/postypes.h /usr/include/c++/4.4/cwchar \
+ /usr/include/c++/4.6/iostream /usr/include/c++/4.6/ostream \
+ /usr/include/c++/4.6/ios /usr/include/c++/4.6/iosfwd \
+ /usr/include/c++/4.6/bits/stringfwd.h \
+ /usr/include/c++/4.6/bits/postypes.h /usr/include/c++/4.6/cwchar \
 /usr/include/bits/wchar.h /usr/include/xlocale.h \
- /usr/include/bits/wchar2.h /usr/include/c++/4.4/bits/char_traits.h \
- /usr/include/c++/4.4/bits/localefwd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++locale.h \
- /usr/include/c++/4.4/clocale /usr/include/locale.h \
- /usr/include/bits/locale.h /usr/include/c++/4.4/cctype \
+ /usr/include/bits/wchar2.h /usr/include/c++/4.6/bits/char_traits.h \
+ /usr/include/c++/4.6/bits/localefwd.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++locale.h \
+ /usr/include/c++/4.6/clocale /usr/include/locale.h \
+ /usr/include/bits/locale.h /usr/include/c++/4.6/cctype \
 /usr/include/ctype.h /usr/include/endian.h /usr/include/bits/endian.h \
- /usr/include/bits/byteswap.h /usr/include/c++/4.4/bits/ios_base.h \
- /usr/include/c++/4.4/ext/atomicity.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr-default.h \
+ /usr/include/bits/byteswap.h /usr/include/c++/4.6/bits/ios_base.h \
+ /usr/include/c++/4.6/ext/atomicity.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr-default.h \
 /usr/include/pthread.h /usr/include/sched.h /usr/include/time.h \
 /usr/include/bits/sched.h /usr/include/bits/time.h /usr/include/signal.h \
 /usr/include/bits/sigset.h /usr/include/bits/pthreadtypes.h \
@ -56,24 +56,24 @@ one2many.o: one2many.cpp /usr/lib/openmpi/include/mpi.h \
 /usr/include/bits/posix_opt.h /usr/include/bits/environments.h \
 /usr/include/bits/confname.h /usr/include/getopt.h \
 /usr/include/bits/unistd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/atomic_word.h \
- /usr/include/c++/4.4/bits/locale_classes.h /usr/include/c++/4.4/string \
- /usr/include/c++/4.4/bits/ostream_insert.h \
- /usr/include/c++/4.4/cxxabi-forced.h \
- /usr/include/c++/4.4/bits/basic_string.h \
- /usr/include/c++/4.4/bits/basic_string.tcc \
- /usr/include/c++/4.4/bits/locale_classes.tcc \
- /usr/include/c++/4.4/streambuf /usr/include/c++/4.4/bits/streambuf.tcc \
- /usr/include/c++/4.4/bits/basic_ios.h \
- /usr/include/c++/4.4/bits/locale_facets.h /usr/include/c++/4.4/cwctype \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/atomic_word.h \
+ /usr/include/c++/4.6/bits/locale_classes.h /usr/include/c++/4.6/string \
+ /usr/include/c++/4.6/bits/ostream_insert.h \
+ /usr/include/c++/4.6/bits/cxxabi_forced.h \
+ /usr/include/c++/4.6/bits/basic_string.h \
+ /usr/include/c++/4.6/bits/basic_string.tcc \
+ /usr/include/c++/4.6/bits/locale_classes.tcc \
+ /usr/include/c++/4.6/streambuf /usr/include/c++/4.6/bits/streambuf.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.h \
+ /usr/include/c++/4.6/bits/locale_facets.h /usr/include/c++/4.6/cwctype \
 /usr/include/wctype.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_base.h \
- /usr/include/c++/4.4/bits/streambuf_iterator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_inline.h \
- /usr/include/c++/4.4/bits/locale_facets.tcc \
- /usr/include/c++/4.4/bits/basic_ios.tcc \
- /usr/include/c++/4.4/bits/ostream.tcc /usr/include/c++/4.4/istream \
- /usr/include/c++/4.4/bits/istream.tcc \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_base.h \
+ /usr/include/c++/4.6/bits/streambuf_iterator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_inline.h \
+ /usr/include/c++/4.6/bits/locale_facets.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.tcc \
+ /usr/include/c++/4.6/bits/ostream.tcc /usr/include/c++/4.6/istream \
+ /usr/include/c++/4.6/bits/istream.tcc \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/constants.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/functions.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/datatype.h \
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/one2many.o
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/one2many.o
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/send2one.d
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/send2one.d
@ -1,54 +1,54 @@
 send2one.o: send2one.cpp /usr/lib/openmpi/include/mpi.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stddef.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stddef.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/mpicxx.h \
- /usr/include/c++/4.4/map /usr/include/c++/4.4/bits/stl_tree.h \
- /usr/include/c++/4.4/bits/stl_algobase.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++config.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/os_defines.h \
+ /usr/include/c++/4.6/map /usr/include/c++/4.6/bits/stl_tree.h \
+ /usr/include/c++/4.6/bits/stl_algobase.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++config.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/os_defines.h \
 /usr/include/features.h /usr/include/bits/predefs.h \
 /usr/include/sys/cdefs.h /usr/include/bits/wordsize.h \
 /usr/include/gnu/stubs.h /usr/include/gnu/stubs-64.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/cpu_defines.h \
- /usr/include/c++/4.4/cstddef /usr/include/c++/4.4/bits/functexcept.h \
- /usr/include/c++/4.4/exception_defines.h \
- /usr/include/c++/4.4/bits/cpp_type_traits.h \
- /usr/include/c++/4.4/ext/type_traits.h \
- /usr/include/c++/4.4/ext/numeric_traits.h \
- /usr/include/c++/4.4/bits/stl_pair.h /usr/include/c++/4.4/bits/move.h \
- /usr/include/c++/4.4/bits/concept_check.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_types.h \
- /usr/include/c++/4.4/bits/stl_iterator_base_funcs.h \
- /usr/include/c++/4.4/bits/stl_iterator.h \
- /usr/include/c++/4.4/debug/debug.h /usr/include/c++/4.4/bits/allocator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++allocator.h \
- /usr/include/c++/4.4/ext/new_allocator.h /usr/include/c++/4.4/new \
- /usr/include/c++/4.4/exception /usr/include/c++/4.4/bits/stl_function.h \
- /usr/include/c++/4.4/backward/binders.h \
- /usr/include/c++/4.4/bits/stl_map.h \
- /usr/include/c++/4.4/initializer_list \
- /usr/include/c++/4.4/bits/stl_multimap.h /usr/include/c++/4.4/utility \
- /usr/include/c++/4.4/bits/stl_relops.h \
- /usr/lib/gcc/x86_64-linux-gnu/4.4.3/include/stdarg.h \
- /usr/include/stdio.h /usr/include/bits/types.h \
- /usr/include/bits/typesizes.h /usr/include/libio.h \
- /usr/include/_G_config.h /usr/include/wchar.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/cpu_defines.h \
+ /usr/include/c++/4.6/bits/functexcept.h \
+ /usr/include/c++/4.6/bits/exception_defines.h \
+ /usr/include/c++/4.6/bits/cpp_type_traits.h \
+ /usr/include/c++/4.6/ext/type_traits.h \
+ /usr/include/c++/4.6/ext/numeric_traits.h \
+ /usr/include/c++/4.6/bits/stl_pair.h /usr/include/c++/4.6/bits/move.h \
+ /usr/include/c++/4.6/bits/concept_check.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_types.h \
+ /usr/include/c++/4.6/bits/stl_iterator_base_funcs.h \
+ /usr/include/c++/4.6/bits/stl_iterator.h \
+ /usr/include/c++/4.6/debug/debug.h /usr/include/c++/4.6/bits/allocator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++allocator.h \
+ /usr/include/c++/4.6/ext/new_allocator.h /usr/include/c++/4.6/new \
+ /usr/include/c++/4.6/exception /usr/include/c++/4.6/bits/stl_function.h \
+ /usr/include/c++/4.6/backward/binders.h \
+ /usr/include/c++/4.6/bits/stl_map.h \
+ /usr/include/c++/4.6/initializer_list \
+ /usr/include/c++/4.6/bits/stl_multimap.h \
+ /usr/include/c++/4.6/bits/range_access.h /usr/include/c++/4.6/utility \
+ /usr/include/c++/4.6/bits/stl_relops.h \
+ /usr/lib/gcc/x86_64-linux-gnu/4.6/include/stdarg.h /usr/include/stdio.h \
+ /usr/include/bits/types.h /usr/include/bits/typesizes.h \
+ /usr/include/libio.h /usr/include/_G_config.h /usr/include/wchar.h \
 /usr/include/bits/stdio_lim.h /usr/include/bits/sys_errlist.h \
 /usr/include/bits/stdio.h /usr/include/bits/stdio2.h \
- /usr/include/c++/4.4/iostream /usr/include/c++/4.4/ostream \
- /usr/include/c++/4.4/ios /usr/include/c++/4.4/iosfwd \
- /usr/include/c++/4.4/bits/stringfwd.h \
- /usr/include/c++/4.4/bits/postypes.h /usr/include/c++/4.4/cwchar \
+ /usr/include/c++/4.6/iostream /usr/include/c++/4.6/ostream \
+ /usr/include/c++/4.6/ios /usr/include/c++/4.6/iosfwd \
+ /usr/include/c++/4.6/bits/stringfwd.h \
+ /usr/include/c++/4.6/bits/postypes.h /usr/include/c++/4.6/cwchar \
 /usr/include/bits/wchar.h /usr/include/xlocale.h \
- /usr/include/bits/wchar2.h /usr/include/c++/4.4/bits/char_traits.h \
- /usr/include/c++/4.4/bits/localefwd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/c++locale.h \
- /usr/include/c++/4.4/clocale /usr/include/locale.h \
- /usr/include/bits/locale.h /usr/include/c++/4.4/cctype \
+ /usr/include/bits/wchar2.h /usr/include/c++/4.6/bits/char_traits.h \
+ /usr/include/c++/4.6/bits/localefwd.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/c++locale.h \
+ /usr/include/c++/4.6/clocale /usr/include/locale.h \
+ /usr/include/bits/locale.h /usr/include/c++/4.6/cctype \
 /usr/include/ctype.h /usr/include/endian.h /usr/include/bits/endian.h \
- /usr/include/bits/byteswap.h /usr/include/c++/4.4/bits/ios_base.h \
- /usr/include/c++/4.4/ext/atomicity.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/gthr-default.h \
+ /usr/include/bits/byteswap.h /usr/include/c++/4.6/bits/ios_base.h \
+ /usr/include/c++/4.6/ext/atomicity.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/gthr-default.h \
 /usr/include/pthread.h /usr/include/sched.h /usr/include/time.h \
 /usr/include/bits/sched.h /usr/include/bits/time.h /usr/include/signal.h \
 /usr/include/bits/sigset.h /usr/include/bits/pthreadtypes.h \
@ -56,24 +56,24 @@ send2one.o: send2one.cpp /usr/lib/openmpi/include/mpi.h \
 /usr/include/bits/posix_opt.h /usr/include/bits/environments.h \
 /usr/include/bits/confname.h /usr/include/getopt.h \
 /usr/include/bits/unistd.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/atomic_word.h \
- /usr/include/c++/4.4/bits/locale_classes.h /usr/include/c++/4.4/string \
- /usr/include/c++/4.4/bits/ostream_insert.h \
- /usr/include/c++/4.4/cxxabi-forced.h \
- /usr/include/c++/4.4/bits/basic_string.h \
- /usr/include/c++/4.4/bits/basic_string.tcc \
- /usr/include/c++/4.4/bits/locale_classes.tcc \
- /usr/include/c++/4.4/streambuf /usr/include/c++/4.4/bits/streambuf.tcc \
- /usr/include/c++/4.4/bits/basic_ios.h \
- /usr/include/c++/4.4/bits/locale_facets.h /usr/include/c++/4.4/cwctype \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/atomic_word.h \
+ /usr/include/c++/4.6/bits/locale_classes.h /usr/include/c++/4.6/string \
+ /usr/include/c++/4.6/bits/ostream_insert.h \
+ /usr/include/c++/4.6/bits/cxxabi_forced.h \
+ /usr/include/c++/4.6/bits/basic_string.h \
+ /usr/include/c++/4.6/bits/basic_string.tcc \
+ /usr/include/c++/4.6/bits/locale_classes.tcc \
+ /usr/include/c++/4.6/streambuf /usr/include/c++/4.6/bits/streambuf.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.h \
+ /usr/include/c++/4.6/bits/locale_facets.h /usr/include/c++/4.6/cwctype \
 /usr/include/wctype.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_base.h \
- /usr/include/c++/4.4/bits/streambuf_iterator.h \
- /usr/include/c++/4.4/x86_64-linux-gnu/bits/ctype_inline.h \
- /usr/include/c++/4.4/bits/locale_facets.tcc \
- /usr/include/c++/4.4/bits/basic_ios.tcc \
- /usr/include/c++/4.4/bits/ostream.tcc /usr/include/c++/4.4/istream \
- /usr/include/c++/4.4/bits/istream.tcc \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_base.h \
+ /usr/include/c++/4.6/bits/streambuf_iterator.h \
+ /usr/include/c++/4.6/x86_64-linux-gnu/./bits/ctype_inline.h \
+ /usr/include/c++/4.6/bits/locale_facets.tcc \
+ /usr/include/c++/4.6/bits/basic_ios.tcc \
+ /usr/include/c++/4.6/bits/ostream.tcc /usr/include/c++/4.6/istream \
+ /usr/include/c++/4.6/bits/istream.tcc \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/constants.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/functions.h \
 /usr/lib/openmpi/include/openmpi/ompi/mpi/cxx/datatype.h \
--- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/send2one.o
+++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/twoWayM2M/library/send2one.o
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/Archimedes/Archimedes.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/Archimedes/Archimedes.C
@ -100,13 +100,7 @@ Archimedes::~Archimedes()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void Archimedes::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void Archimedes::setForce() const
 {
    vector force(0,0,0);

@ -129,9 +123,14 @@ void Archimedes::setForce
                }
            }

-            if(treatDEM_) for(int j=0;j<3;j++) DEMForces[index][j] += force[j];
-            else if(treatExplicit_) for(int j=0;j<3;j++) expForces[index][j] += force[j];
-            else  for(int j=0;j<3;j++) impForces[index][j] += force[j];
+            if(!treatDEM_)
+            {            
+                if(treatExplicit_)
+                    for(int j=0;j<3;j++) expForces()[index][j] += force[j];
+                else
+                    for(int j=0;j<3;j++) impForces()[index][j] += force[j];
+            }
+            for(int j=0;j<3;j++) DEMForces()[index][j] += force[j];
        //}
    }
 }
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/Archimedes/Archimedes.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/Archimedes/Archimedes.H
@ -97,13 +97,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/ArchimedesIB/ArchimedesIB.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/ArchimedesIB/ArchimedesIB.C
@ -81,14 +81,8 @@ ArchimedesIB::ArchimedesIB
    }

    if (propsDict_.found("treatExplicit")) treatExplicit_=true;
-    if (modelType_=="A"){
-        treatDEM_=true;
-        Info << "accounting for Archimedes only on DEM side!" << endl;
-    }
-    if (modelType_=="B"){
-        treatDEM_=false;
-        Info << "accounting for Archimedes on DEM and CFD side!" << endl;
-    }
+    treatDEM_=true;
+    Info << "accounting for Archimedes only on DEM side!" << endl;
    particleCloud_.checkCG(true);
 }

@ -101,13 +95,7 @@ ArchimedesIB::~ArchimedesIB()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void ArchimedesIB::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void ArchimedesIB::setForce() const
 {
    vector force;
    for(int index = 0;index <  particleCloud_.numberOfParticles(); ++index)
@ -126,9 +114,12 @@ void ArchimedesIB::setForce
            }
            // set force on particle
            if(twoDimensional_) Warning<<"ArchimedesIB model doesn't work for 2D right now!!\n"<< endl;
-            if(treatDEM_) for(int j=0;j<3;j++) DEMForces[index][j] += force[j];
-            else if(treatExplicit_) for(int j=0;j<3;j++) expForces[index][j] += force[j];
-            else  for(int j=0;j<3;j++) impForces[index][j] += force[j];
+            if(!treatDEM_)
+            {
+                if(treatExplicit_) for(int j=0;j<3;j++) expForces()[index][j] += force[j];
+                else for(int j=0;j<3;j++) impForces()[index][j] += force[j];
+            }
+            for(int j=0;j<3;j++) DEMForces()[index][j] += force[j];
        //}
    }
 }
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/ArchimedesIB/ArchimedesIB.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/ArchimedesIB/ArchimedesIB.H
@ -103,13 +103,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/DiFeliceDrag/DiFeliceDrag.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/DiFeliceDrag/DiFeliceDrag.C
@ -92,13 +92,7 @@ DiFeliceDrag::~DiFeliceDrag()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void DiFeliceDrag::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void DiFeliceDrag::setForce() const
 {
    // get viscosity field
    #ifdef comp
@ -188,8 +182,9 @@ void DiFeliceDrag::setForce
                }
            }
            // set force on particle
-            if(treatExplicit_) for(int j=0;j<3;j++) expForces[index][j] += drag[j];
-            else  for(int j=0;j<3;j++) impForces[index][j] += drag[j];
+            if(treatExplicit_) for(int j=0;j<3;j++) expForces()[index][j] += drag[j];
+            else  for(int j=0;j<3;j++) impForces()[index][j] += drag[j];
+            for(int j=0;j<3;j++) DEMForces()[index][j] += drag[j];
        }
    //}
 }
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/DiFeliceDrag/DiFeliceDrag.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/DiFeliceDrag/DiFeliceDrag.H
@ -97,13 +97,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/GidaspowDrag/GidaspowDrag.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/GidaspowDrag/GidaspowDrag.C
@ -83,13 +83,7 @@ GidaspowDrag::~GidaspowDrag()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void GidaspowDrag::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void GidaspowDrag::setForce() const
 {
    // get viscosity field
    #ifdef comp
@ -143,8 +137,9 @@ void GidaspowDrag::setForce
            }

            // set force on particle
-            if(treatExplicit_) for(int j=0;j<3;j++) expForces[index][j] += drag[j];
-            else  for(int j=0;j<3;j++) impForces[index][j] += drag[j];
+            if(treatExplicit_) for(int j=0;j<3;j++) expForces()[index][j] += drag[j];
+            else  for(int j=0;j<3;j++) impForces()[index][j] += drag[j];
+            for(int j=0;j<3;j++) DEMForces()[index][j] += drag[j];
        //}// end if mask
    }// end loop particles
 }
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/GidaspowDrag/GidaspowDrag.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/GidaspowDrag/GidaspowDrag.H
@ -90,13 +90,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/KochHillDrag/KochHillDrag.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/KochHillDrag/KochHillDrag.C
@ -35,7 +35,7 @@ Description
 #include "addToRunTimeSelectionTable.H"
 #include "dataExchangeModel.H"

-#include "mpi.h"
+//#include "mpi.h"
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

 namespace Foam
@ -76,6 +76,13 @@ KochHillDrag::KochHillDrag
    if (propsDict_.found("verbose")) verbose_=true;
    if (propsDict_.found("treatExplicit")) treatExplicit_=true;
    if (propsDict_.found("interpolation")) interpolation_=true;
+    if (propsDict_.found("implDEM"))
+    {
+        treatExplicit_=false;
+        implDEM_=true;
+        setImpDEMdrag();
+        Info << "Using implicit DEM drag formulation." << endl;
+    }
    particleCloud_.checkCG(true);
 }

@ -88,13 +95,7 @@ KochHillDrag::~KochHillDrag()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void KochHillDrag::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void KochHillDrag::setForce() const
 {
    // get viscosity field
    #ifdef comp
@ -118,6 +119,7 @@ void KochHillDrag::setForce
    scalar Rep(0);
 	scalar Vs(0);
 	scalar volumefraction(0);
+    scalar betaP(0);

    interpolationCellPoint<scalar> voidfractionInterpolator_(voidfraction_);
    interpolationCellPoint<vector> UInterpolator_(U_);
@ -128,6 +130,9 @@ void KochHillDrag::setForce
        //{
            cellI = particleCloud_.cellIDs()[index][0];
            drag = vector(0,0,0);
+            betaP = 0;
+            Vs = 0;
+            Ufluid =vector(0,0,0);

            if (cellI > -1) // particle Found
            {
@ -182,7 +187,7 @@ void KochHillDrag::setForce
                    scalar F = voidfraction * (F0 + 0.5*F3*Rep);

                    // calc drag model coefficient betaP
-                    scalar betaP = 18.*nuf*rho/(ds/cg()*ds/cg())*voidfraction*F;
+                    betaP = 18.*nuf*rho/(ds/cg()*ds/cg())*voidfraction*F;

                    // calc particle's drag
                    drag = Vs*betaP*Ur;
@ -206,8 +211,23 @@ void KochHillDrag::setForce
                }
            }
            // set force on particle
-            if(treatExplicit_) for(int j=0;j<3;j++) expForces[index][j] += drag[j];
-            else  for(int j=0;j<3;j++) impForces[index][j] += drag[j];
+            if(treatExplicit_) for(int j=0;j<3;j++) expForces()[index][j] += drag[j];
+            else  for(int j=0;j<3;j++) impForces()[index][j] += drag[j];
+
+            // set Cd
+            if(implDEM_)
+            {
+                for(int j=0;j<3;j++) fluidVel()[index][j]=Ufluid[j];
+
+                if (modelType_=="B")
+                    Cds()[index][0] = Vs*betaP/voidfraction;
+                else
+                    Cds()[index][0] = Vs*betaP;
+
+            }else{
+                for(int j=0;j<3;j++) DEMForces()[index][j] += drag[j];
+            }
+
        //}
    }
 }
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/KochHillDrag/KochHillDrag.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/KochHillDrag/KochHillDrag.H
@ -32,6 +32,7 @@ Description
    based on Koch Hill 2001,"Inertial effects in suspensions and porous-media 
    flows", Annual Review of fluid mechanics.
    including interpolation of the velocity to the exact position
+    including drag coefficient for implicit drag for DEM

 Class
    KochHillDrag
@ -100,13 +101,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/LaEuScalarTemp/LaEuScalarTemp.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/LaEuScalarTemp/LaEuScalarTemp.C
@ -108,13 +108,7 @@ void LaEuScalarTemp::allocateMyArrays() const
 }
 // * * * * * * * * * * * * * * * public Member Functions  * * * * * * * * * * * * * //

-void LaEuScalarTemp::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void LaEuScalarTemp::setForce() const
 {
    // do nothing
 }
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/LaEuScalarTemp/LaEuScalarTemp.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/LaEuScalarTemp/LaEuScalarTemp.H
@ -116,13 +116,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;

        void manipulateScalarField(volScalarField&) const;
 };
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/MeiLift/MeiLift.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/MeiLift/MeiLift.C
@ -86,13 +86,7 @@ MeiLift::~MeiLift()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void MeiLift::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void MeiLift::setForce() const
 {
    // get viscosity field
    #ifdef comp
@ -171,9 +165,11 @@ void MeiLift::setForce
                }
            }
            // set force on particle
-            if(treatDEM_) for(int j=0;j<3;j++) DEMForces[index][j] += lift[j];
-            else if(!treatExplicit_) for(int j=0;j<3;j++) impForces[index][j] += lift[j];
-            else  for(int j=0;j<3;j++) expForces[index][j] += lift[j];
+            if(!treatDEM_){
+                if(!treatExplicit_) for(int j=0;j<3;j++) impForces()[index][j] += lift[j];
+                else  for(int j=0;j<3;j++) expForces()[index][j] += lift[j];
+            }
+            for(int j=0;j<3;j++) DEMForces()[index][j] += lift[j];
        //}
    }

--- a/src/lagrangian/cfdemParticle/subModels/forceModel/MeiLift/MeiLift.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/MeiLift/MeiLift.H
@ -96,13 +96,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/SchillerNaumannDrag/SchillerNaumannDrag.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/SchillerNaumannDrag/SchillerNaumannDrag.C
@ -84,13 +84,7 @@ SchillerNaumannDrag::~SchillerNaumannDrag()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void SchillerNaumannDrag::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void SchillerNaumannDrag::setForce() const
 {
    // get viscosity field
    #ifdef comp
@ -149,8 +143,9 @@ void SchillerNaumannDrag::setForce
                }
            }
            // set force on particle
-            if(treatExplicit_) for(int j=0;j<3;j++) expForces[index][j] += drag[j];
-            else  for(int j=0;j<3;j++) impForces[index][j] += drag[j];
+            if(treatExplicit_) for(int j=0;j<3;j++) expForces()[index][j] += drag[j];
+            else  for(int j=0;j<3;j++) impForces()[index][j] += drag[j];
+            for(int j=0;j<3;j++) DEMForces()[index][j] += drag[j];
        //}
    }

--- a/src/lagrangian/cfdemParticle/subModels/forceModel/SchillerNaumannDrag/SchillerNaumannDrag.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/SchillerNaumannDrag/SchillerNaumannDrag.H
@ -91,13 +91,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/ShirgaonkarIB/ShirgaonkarIB.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/ShirgaonkarIB/ShirgaonkarIB.C
@ -94,13 +94,7 @@ ShirgaonkarIB::~ShirgaonkarIB()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void ShirgaonkarIB::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void ShirgaonkarIB::setForce() const
 {

    label cellI;
@ -136,10 +130,11 @@ void ShirgaonkarIB::setForce

            // set force on particle
            if(twoDimensional_) drag /= depth_;
-            if(treatExplicit_) for(int j=0;j<3;j++) expForces[index][j] += drag[j];
-            else  for(int j=0;j<3;j++) impForces[index][j] += drag[j];
+            if(treatExplicit_) for(int j=0;j<3;j++) expForces()[index][j] += drag[j];
+            else  for(int j=0;j<3;j++) impForces()[index][j] += drag[j];
+            for(int j=0;j<3;j++) DEMForces()[index][j] += drag[j];

-            if(verbose_) Info << "impForces = " << impForces[index][0]<<","<<impForces[index][1]<<","<<impForces[index][2] << endl;
+            if(verbose_) Info << "impForces = " << impForces()[index][0]<<","<<impForces()[index][1]<<","<<impForces()[index][2] << endl;
        //}
    }
 }
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/ShirgaonkarIB/ShirgaonkarIB.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/ShirgaonkarIB/ShirgaonkarIB.H
@ -99,13 +99,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/checkCouplingInterval/checkCouplingInterval.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/checkCouplingInterval/checkCouplingInterval.C
@ -78,13 +78,7 @@ checkCouplingInterval::~checkCouplingInterval()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void checkCouplingInterval::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void checkCouplingInterval::setForce() const
 {
    if(particleCloud_.mesh().time().write())
    {
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/checkCouplingInterval/checkCouplingInterval.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/checkCouplingInterval/checkCouplingInterval.H
@ -89,13 +89,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/forceModel/forceModel.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/forceModel/forceModel.C
@ -58,6 +58,7 @@ forceModel::forceModel
    particleCloud_(sm),
    treatExplicit_(false),
    treatDEM_(false),
+    implDEM_(false),
    impParticleForces_
    (   IOobject
        (
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/forceModel/forceModel.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/forceModel/forceModel.H
@ -64,6 +64,8 @@ protected:

        Switch treatDEM_;                           // marker to use the force only on DEM side

+        Switch implDEM_;                           // marker to use the implicit force on DEM side
+
        mutable volVectorField impParticleForces_;  // sum of implicit particle forces [N]

        mutable volVectorField expParticleForces_;  // sum of explicit particle forces [N]
@ -119,7 +121,7 @@ public:


    // Member Functions
-        virtual void setForce(double** const& mask,double**&,double**&,double**&) const = 0;
+        virtual void setForce() const = 0;

        //tmp<volScalarField> provideScalarField();

@ -132,6 +134,16 @@ public:

        inline volVectorField& expParticleForces() const { return expParticleForces_;};

+        inline double ** impForces() const { return particleCloud_.impForces_;};
+
+        inline double ** expForces() const { return particleCloud_.expForces_;};
+
+        inline double ** DEMForces() const { return particleCloud_.DEMForces_;};
+
+        inline double ** Cds() const { return particleCloud_.Cds_;};
+
+        inline double ** fluidVel() const { return particleCloud_.fluidVel_;};
+
        inline const bool& coupleForce() const { return coupleForce_;};

        inline const scalar cg() const { return cg_; };
@ -141,6 +153,10 @@ public:
            cg_ = cg;
            Info << "cg is set to: " << cg_ << endl;
        };
+
+        void const setImpDEMdrag() const {particleCloud_.impDEMdrag_=true;};
+
+
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/gradPForce/gradPForce.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/gradPForce/gradPForce.C
@ -105,13 +105,7 @@ gradPForce::~gradPForce()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void gradPForce::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void gradPForce::setForce() const
 {
    volVectorField gradPField = fvc::grad(p_);
    /*if (useU_)
@ -172,9 +166,11 @@ void gradPForce::setForce
            }

            // set force on particle
-            if(treatDEM_) for(int j=0;j<3;j++) DEMForces[index][j] += force[j];
-            else if(!treatExplicit_) for(int j=0;j<3;j++) impForces[index][j] += force[j];
-            else  for(int j=0;j<3;j++) expForces[index][j] += force[j];
+            if(!treatDEM_){
+                if(!treatExplicit_) for(int j=0;j<3;j++) impForces()[index][j] += force[j];
+                else  for(int j=0;j<3;j++) expForces()[index][j] += force[j];
+            }
+            for(int j=0;j<3;j++) DEMForces()[index][j] += force[j];

        //}
    }
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/gradPForce/gradPForce.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/gradPForce/gradPForce.H
@ -104,13 +104,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/interface/interface.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/interface/interface.C
@ -93,13 +93,7 @@ interface::~interface()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void interface::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void interface::setForce() const
 {
 Info << "interface::setForce" << endl;
    for(int index = 0;index <  particleCloud_.numberOfParticles(); ++index)
@ -185,8 +179,9 @@ Info << "interface::setForce" << endl;
                    Info << "interface force is limited to " << interfaceForce << endl;
                }*/

-               if(treatExplicit_) for(int j=0;j<3;j++) expForces[index][j] += interfaceForce[j];
-               else  for(int j=0;j<3;j++) impForces[index][j] += interfaceForce[j];
+               if(treatExplicit_) for(int j=0;j<3;j++) expForces()[index][j] += interfaceForce[j];
+               else  for(int j=0;j<3;j++) impForces()[index][j] += interfaceForce[j];
+               for(int j=0;j<3;j++) DEMForces()[index][j] += interfaceForce[j];
            } // end if particle found on proc domain
        //}// end if in mask
    }// end loop particles
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/interface/interface.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/interface/interface.H
@ -105,13 +105,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/noDrag/noDrag.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/noDrag/noDrag.C
@ -74,13 +74,7 @@ noDrag::~noDrag()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void noDrag::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void noDrag::setForce() const
 {
    // Do nothing
    Info << "noDrag::setForce" << endl;
@ -89,8 +83,8 @@ void noDrag::setForce
        //if(mask[index][0])
        //{
            // set force on particle
-            if(treatExplicit_) for(int j=0;j<3;j++) expForces[index][j] = 0.;
-            else  for(int j=0;j<3;j++) impForces[index][j] = 0.;
+            if(treatExplicit_) for(int j=0;j<3;j++) expForces()[index][j] = 0.;
+            else  for(int j=0;j<3;j++) impForces()[index][j] = 0.;
        //}
    }
 }
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/noDrag/noDrag.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/noDrag/noDrag.H
@ -76,13 +76,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModel/virtualMassForce/virtualMassForce.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/virtualMassForce/virtualMassForce.C
@ -92,13 +92,7 @@ virtualMassForce::~virtualMassForce()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void virtualMassForce::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void virtualMassForce::setForce() const
 {
    reAllocArrays();

@ -131,8 +125,9 @@ void virtualMassForce::setForce

            }
            // set force on particle
-            if(treatExplicit_) for(int j=0;j<3;j++) expForces[index][j] += virtualMassForce[j];
-            else  for(int j=0;j<3;j++) impForces[index][j] += virtualMassForce[j];
+            if(treatExplicit_) for(int j=0;j<3;j++) expForces()[index][j] += virtualMassForce[j];
+            else  for(int j=0;j<3;j++) impForces()[index][j] += virtualMassForce[j];
+            for(int j=0;j<3;j++) DEMForces()[index][j] += virtualMassForce[j];
        //}
    }

--- a/src/lagrangian/cfdemParticle/subModels/forceModel/virtualMassForce/virtualMassForce.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/virtualMassForce/virtualMassForce.H
@ -91,13 +91,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;

        void reAllocArrays() const;
 };
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/viscForce/viscForce.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/viscForce/viscForce.C
@ -96,13 +96,7 @@ viscForce::~viscForce()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void viscForce::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void viscForce::setForce() const
 {

    // get viscosity field
@ -167,9 +161,11 @@ void viscForce::setForce
            }

            // set force on particle
-            if(treatDEM_) for(int j=0;j<3;j++) DEMForces[index][j] += force[j];
-            else if(!treatExplicit_) for(int j=0;j<3;j++) impForces[index][j] += force[j];
-            else  for(int j=0;j<3;j++) expForces[index][j] += force[j];
+            if(!treatDEM_){
+                if(!treatExplicit_) for(int j=0;j<3;j++) impForces()[index][j] += force[j];
+                else  for(int j=0;j<3;j++) expForces()[index][j] += force[j];
+            }
+             for(int j=0;j<3;j++) DEMForces()[index][j] += force[j];

        //}
    }
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/viscForce/viscForce.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/viscForce/viscForce.H
@ -92,13 +92,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/src/lagrangian/cfdemParticle/subModels/forceModelMS/DiFeliceDragMS/DiFeliceDragMS.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModelMS/DiFeliceDragMS/DiFeliceDragMS.C
@ -93,13 +93,7 @@ DiFeliceDragMS::~DiFeliceDragMS()

 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

-void DiFeliceDragMS::setForce
-(
-    double** const& mask,
-    double**& impForces,
-    double**& expForces,
-    double**& DEMForces
-) const
+void DiFeliceDragMS::setForce() const
 {
    // get viscosity field
    #ifdef comp
@ -218,8 +212,8 @@ void DiFeliceDragMS::setForce
                Warning <<"A BUG occurred in DiFeliceDragMS::setForce!!! nrigidC = " << nrigidC <<", ind = " << ind <<", index=" << index <<"\n" << endl;
                nrigidC = 1000;
            }
-            if(treatExplicit_) for(int j=0;j<3;j++) expForces[index][j] += cloudRefMS().expForcesCM()[ind][j] / nrigidC;
-            else  for(int j=0;j<3;j++) impForces[index][j] += cloudRefMS().impForcesCM()[ind][j] / nrigidC;
+            if(treatExplicit_) for(int j=0;j<3;j++) expForces()[index][j] += cloudRefMS().expForcesCM()[ind][j] / nrigidC;
+            else  for(int j=0;j<3;j++) impForces()[index][j] += cloudRefMS().impForcesCM()[ind][j] / nrigidC;
        }
    }
 }
--- a/src/lagrangian/cfdemParticle/subModels/forceModelMS/DiFeliceDragMS/DiFeliceDragMS.H
+++ b/src/lagrangian/cfdemParticle/subModels/forceModelMS/DiFeliceDragMS/DiFeliceDragMS.H
@ -100,13 +100,7 @@ public:


    // Member Functions
-        void setForce
-        (
-            double** const& mask,
-            double**&,
-            double**&,
-            double**&
-        ) const;
+        void setForce() const;
 };


--- a/tutorials/cfdemSolverPiso/ErgunTestMPI/DEM/liggghts.restart
+++ b/tutorials/cfdemSolverPiso/ErgunTestMPI/DEM/liggghts.restart
--- a/tutorials/cfdemSolverPiso/ErgunTestMPI_restart/DEM/liggghts.restart
+++ b/tutorials/cfdemSolverPiso/ErgunTestMPI_restart/DEM/liggghts.restart
--- a/tutorials/cfdemSolverPiso/ErgunTestMPI_restart/DEM/post/dummy
+++ b/tutorials/cfdemSolverPiso/ErgunTestMPI_restart/DEM/post/dummy
@ -37,3 +37,4 @@ dummyfile
 dummyfile
 dummyfile
 dummyfile
+dummyfile
--- a/tutorials/cfdemSolverPiso/settlingTestMPI/CFD/constant/couplingProperties
+++ b/tutorials/cfdemSolverPiso/settlingTestMPI/CFD/constant/couplingProperties
@ -38,7 +38,7 @@ meshMotionModel noMeshMotion;

 regionModel allRegion;

-IOModel basicIO; //trackIO; // 
+IOModel trackIO; //basicIO; // 

 dataExchangeModel twoWayMPI;//twoWayFiles;//oneWayVTK;//

--- a/tutorials/cfdemSolverPiso/settlingTestMPI/CFD/constant/particleTrackProperties
+++ b/tutorials/cfdemSolverPiso/settlingTestMPI/CFD/constant/particleTrackProperties
@ -15,7 +15,7 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

-cloudName       cfdemCloud1;
+cloudName       particleCloud;

 sampleFrequency 1;