From 51eb1ac5b643fc9f8885973e99951bb09395e658 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 7 Nov 2011 14:25:51 +0000
Subject: [PATCH] BUG: chemFoam - density output on start-up was misleading due
to presence of call to thermo.correct()
- corresponds to mantis bug report 333
---
.../combustion/chemFoam/readInitialConditions.H | 17 ++++++++---------
1 file changed, 8 insertions(+), 9 deletions(-)
diff --git a/applications/solvers/combustion/chemFoam/readInitialConditions.H b/applications/solvers/combustion/chemFoam/readInitialConditions.H
index 817ed264f9a..3fbb13a0fce 100644
--- a/applications/solvers/combustion/chemFoam/readInitialConditions.H
+++ b/applications/solvers/combustion/chemFoam/readInitialConditions.H
@@ -1,13 +1,5 @@
word constProp(initialConditions.lookup("constantProperty"));
- if (constProp == "pressure" || constProp == "volume")
- {
- Info << constProp << " will be held constant." << nl
- << " p = " << p[0] << " [Pa]" << nl
- << " T = " << thermo.T()[0] << " [K] " << nl
- << " rho = " << rho[0] << " [kg/m3]" << nl
- << endl;
- }
- else
+ if ((constProp != "pressure") && (constProp != "volume"))
{
FatalError << "in initialConditions, unknown constantProperty type "
<< constProp << nl << " Valid types are: pressure volume."
@@ -110,3 +102,10 @@
scalar integratedHeat = 0.0;
+ Info << constProp << " will be held constant." << nl
+ << " p = " << p[0] << " [Pa]" << nl
+ << " T = " << thermo.T()[0] << " [K] " << nl
+ << " rho = " << rho[0] << " [kg/m3]" << nl
+ << endl;
+
+
--
GitLab