applications/test/readCHEMKINIII/Test-readCHEMKINIII.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 \*---------------------------------------------------------------------------*/
  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.append("CHEMKINIIIFile");
  39     argList::addOption("thermo", "fileName");
  40     argList args(argc, argv);
  41
  42     fileName thermoFileName = fileName::null;
  43     args.optionReadIfPresent("thermo", thermoFileName);
  44
  45     chemkinReader ck(args[1], thermoFileName);
  46
  47     //Info<< ck.isotopeAtomicWts() << nl
  48     //    << ck.specieNames() << nl
  49     //    << ck.speciePhase() << nl
  50     //    << ck.specieThermo() << nl
  51     //    << ck.reactions() << endl;
  52
  53     const SLPtrList<gasReaction>& reactions = ck.reactions();
  54
  55     {
  56         OFstream reactionStream("reactions");
  57         reactionStream<< reactions << endl;
  58     }
  59
  60     {
  61         IFstream reactionStream("reactions");
  62
  63         label nReactions(readLabel(reactionStream));
  64         reactionStream.readBeginList(args.executable().c_str());
  65
  66         PtrList<gasReaction> testReactions(nReactions);
  67
  68         forAll(testReactions, i)
  69         {
  70             testReactions.set
  71             (
  72                 i,
  73                 gasReaction::New
  74                 (
  75                     ck.species(),
  76                     ck.speciesThermo(),
  77                     reactionStream
  78                 )
  79             );
  80         }
  81
  82         reactionStream.readEndList(args.executable().c_str());
  83
  84         Info<< testReactions << endl;
  85     }
  86
  87     Info<< "End\n" << endl;
  88
  89     return 0;
  90 }
  91
  92
  93 // ************************************************************************* //