applications/test/readCHEMKINIII/readCHEMKINIII.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2010 OpenCFD Ltd.
   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 \*---------------------------------------------------------------------------*/
  25
  26 #include "chemkinReader.H"
  27 #include "argList.H"
  28 #include "IFstream.H"
  29 #include "OFstream.H"
  30
  31 using namespace Foam;
  32
  33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  34 // Main program:
  35
  36 int main(int argc, char *argv[])
  37 {
  38     argList::validArgs.clear();
  39     argList::validArgs.append("CHEMKINIIIFile");
  40     argList::validOptions.insert("thermo", "fileName");
  41     argList args(argc, argv);
  42
  43     fileName thermoFileName = fileName::null;
  44     if (args.options().found("thermo"))
  45     {
  46         thermoFileName = args.options()["thermo"];
  47     }
  48
  49     fileName CHEMKINFileName(args.additionalArgs()[0]);
  50
  51     chemkinReader ck(CHEMKINFileName, thermoFileName);
  52
  53     //Info<< ck.isotopeAtomicWts() << endl;
  54     //Info<< ck.specieNames() << endl;
  55     //Info<< ck.speciePhase() << endl;
  56     //Info<< ck.specieThermo() << endl;
  57     //Info<< ck.reactions() << endl;
  58
  59     const SLPtrList<gasReaction>& reactions = ck.reactions();
  60
  61     {
  62         OFstream reactionStream("reactions");
  63         reactionStream<< reactions << endl;
  64     }
  65
  66     {
  67         IFstream reactionStream("reactions");
  68
  69         label nReactions(readLabel(reactionStream));
  70         reactionStream.readBeginList(args.executable().c_str());
  71
  72         PtrList<gasReaction> testReactions(nReactions);
  73
  74         forAll (testReactions, i)
  75         {
  76             testReactions.set
  77             (
  78                 i,
  79                 gasReaction::New
  80                 (
  81                     ck.species(),
  82                     ck.speciesThermo(),
  83                     reactionStream
  84                 )
  85             );
  86         }
  87
  88         reactionStream.readEndList(args.executable().c_str());
  89
  90         Info<< testReactions << endl;
  91     }
  92
  93     Info << "End\n" << endl;
  94
  95     return 0;
  96 }
  97
  98
  99 // ************************************************************************* //