00001 /*---------------------------------------------------------------------------*\ 00002 ========= | 00003 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox 00004 \\ / O peration | 00005 \\ / A nd | Copyright (C) 1991-2010 OpenCFD Ltd. 00006 \\/ M anipulation | 00007 ------------------------------------------------------------------------------- 00008 License 00009 This file is part of OpenFOAM. 00010 00011 OpenFOAM is free software: you can redistribute it and/or modify it 00012 under the terms of the GNU General Public License as published by 00013 the Free Software Foundation, either version 3 of the License, or 00014 (at your option) any later version. 00015 00016 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT 00017 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 00018 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 00019 for more details. 00020 00021 You should have received a copy of the GNU General Public License 00022 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>. 00023 00024 \*---------------------------------------------------------------------------*/ 00025 00026 #include "basicPsiThermo.H" 00027 00028 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00029 00030 Foam::autoPtr<Foam::basicPsiThermo> Foam::basicPsiThermo::New 00031 ( 00032 const fvMesh& mesh 00033 ) 00034 { 00035 word thermoTypeName; 00036 00037 // Enclose the creation of the thermophysicalProperties to ensure it is 00038 // deleted before the turbulenceModel is created otherwise the dictionary 00039 // is entered in the database twice 00040 { 00041 IOdictionary thermoDict 00042 ( 00043 IOobject 00044 ( 00045 "thermophysicalProperties", 00046 mesh.time().constant(), 00047 mesh, 00048 IOobject::MUST_READ, 00049 IOobject::NO_WRITE 00050 ) 00051 ); 00052 00053 thermoDict.lookup("thermoType") >> thermoTypeName; 00054 } 00055 00056 Info<< "Selecting thermodynamics package " << thermoTypeName << endl; 00057 00058 fvMeshConstructorTable::iterator cstrIter = 00059 fvMeshConstructorTablePtr_->find(thermoTypeName); 00060 00061 if (cstrIter == fvMeshConstructorTablePtr_->end()) 00062 { 00063 FatalErrorIn("basicPsiThermo::New(const fvMesh&)") 00064 << "Unknown basicPsiThermo type " << thermoTypeName << nl << nl 00065 << "Valid basicPsiThermo types are:" << nl 00066 << fvMeshConstructorTablePtr_->sortedToc() << nl 00067 << exit(FatalError); 00068 } 00069 00070 return autoPtr<basicPsiThermo>(cstrIter()(mesh)); 00071 } 00072 00073 00074 // ************************ vim: set sw=4 sts=4 et: ************************ //