fixed writing out entries in advective bc
[OpenFOAM-1.6-ext.git] / src / thermophysicalModels / basic / rhoThermo / basicRhoThermo / newBasicRhoThermo.C
blob2bd8c686a3a7a4bd654ba83d4e9ae118dadbb19b
1 /*---------------------------------------------------------------------------*\
2   =========                 |
3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
4    \\    /   O peration     |
5     \\  /    A nd           | Copyright held by original author
6      \\/     M anipulation  |
7 -------------------------------------------------------------------------------
8 License
9     This file is part of OpenFOAM.
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 the
13     Free Software Foundation; either version 2 of the License, or (at your
14     option) any later version.
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.
21     You should have received a copy of the GNU General Public License
22     along with OpenFOAM; if not, write to the Free Software Foundation,
23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
25 Description
26     Selection function for 'basic' density-based thermodynamics
28 \*---------------------------------------------------------------------------*/
30 #include "basicRhoThermo.H"
32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
34 Foam::autoPtr<Foam::basicRhoThermo> Foam::basicRhoThermo::New
36     const fvMesh& mesh
39     word thermoTypeName;
41     // Enclose the creation of the thermophysicalProperties to ensure it is
42     // deleted before the turbulenceModel is created otherwise the dictionary
43     // is entered in the database twice
44     {
45         IOdictionary thermoDict
46         (
47             IOobject
48             (
49                 "thermophysicalProperties",
50                 mesh.time().constant(),
51                 mesh,
52                 IOobject::MUST_READ,
53                 IOobject::NO_WRITE
54             )
55         );
57         thermoDict.lookup("thermoType") >> thermoTypeName;
58     }
60     Info<< "Selecting thermodynamics package " << thermoTypeName << endl;
62     fvMeshConstructorTable::iterator cstrIter =
63         fvMeshConstructorTablePtr_->find(thermoTypeName);
65     if (cstrIter == fvMeshConstructorTablePtr_->end())
66     {
67         FatalErrorIn("basicRhoThermo::New(const fvMesh&)")
68             << "Unknown basicRhoThermo type " << thermoTypeName << nl << nl
69             << "Valid basicRhoThermo types are:" << nl
70             << fvMeshConstructorTablePtr_->toc() << nl
71             << exit(FatalError);
72     }
74     return autoPtr<basicRhoThermo>(cstrIter()(mesh));
78 // ************************************************************************* //