applications/utilities/preProcessing/mdInitialise/mdInitialise.C

   1 /*---------------------------------------------------------------------------*\
   2  =========                   |
   3  \\      /   F ield          | OpenFOAM: The Open Source CFD Toolbox
   4   \\    /    O peration      |
   5    \\  /     A nd            | Copyright (C) 2011 OpenFOAM Foundation
   6     \\/      M anipulation   |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software: you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Description
  25     Initialises fields for a molecular dynamics (MD) simulation.
  26
  27 \*---------------------------------------------------------------------------*/
  28
  29 #include "md.H"
  30 #include "fvCFD.H"
  31
  32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  33
  34 int main(int argc, char *argv[])
  35 {
  36 #   include "setRootCase.H"
  37 #   include "createTime.H"
  38 #   include "createMesh.H"
  39
  40     IOdictionary mdInitialiseDict
  41     (
  42         IOobject
  43         (
  44             "mdInitialiseDict",
  45             runTime.system(),
  46             runTime,
  47             IOobject::MUST_READ_IF_MODIFIED,
  48             IOobject::NO_WRITE,
  49             false
  50         )
  51     );
  52
  53     IOdictionary idListDict
  54     (
  55         IOobject
  56         (
  57             "idList",
  58             mesh.time().constant(),
  59             mesh,
  60             IOobject::NO_READ,
  61             IOobject::AUTO_WRITE
  62         )
  63     );
  64
  65     potential pot(mesh, mdInitialiseDict, idListDict);
  66
  67     moleculeCloud molecules(mesh, pot, mdInitialiseDict);
  68
  69     label totalMolecules = molecules.size();
  70
  71     if (Pstream::parRun())
  72     {
  73         reduce(totalMolecules, sumOp<label>());
  74     }
  75
  76     Info<< nl << "Total number of molecules added: " << totalMolecules
  77         << nl << endl;
  78
  79     IOstream::defaultPrecision(15);
  80
  81     if (!mesh.write())
  82     {
  83         FatalErrorIn(args.executable())
  84             << "Failed writing moleculeCloud."
  85             << nl << exit(FatalError);
  86     }
  87
  88     Info<< nl << "ClockTime = " << runTime.elapsedClockTime() << " s"
  89         << nl << endl;
  90
  91     Info<< "\nEnd\n" << endl;
  92
  93     return 0;
  94 }
  95
  96
  97 // ************************************************************************* //