From fa851dde1e20b8b5d22a86628f69c039910b75cb Mon Sep 17 00:00:00 2001
From: graham <g.macpherson@opencfd.co.uk>
Date: Wed, 14 Apr 2010 16:03:11 +0100
Subject: [PATCH] BUG: incompressibleUncoupledKinematicParcelFoam. Using
headerOk to exclude H and HdotGradH fields from being read if not present.
---
.../createFields.H | 48 ++++++++++---------
1 file changed, 26 insertions(+), 22 deletions(-)
diff --git a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/createFields.H b/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/createFields.H
index e26fddee0b2..5ac8e24ef49 100644
--- a/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/createFields.H
+++ b/applications/solvers/lagrangian/incompressibleUncoupledKinematicParcelFoam/createFields.H
@@ -84,30 +84,34 @@
g
);
- Info<< "\nReading field H\n" << endl;
- volVectorField H
+ IOobject Hheader
(
- IOobject
- (
- "H",
- runTime.timeName(),
- mesh,
- IOobject::READ_IF_PRESENT,
- IOobject::AUTO_WRITE
- ),
- mesh
+ "H",
+ runTime.timeName(),
+ mesh,
+ IOobject::MUST_READ,
+ IOobject::AUTO_WRITE
);
- Info<< "\nReading field HdotGradH\n" << endl;
- volVectorField HdotGradH
+ if (Hheader.headerOk())
+ {
+
+ Info<< "\nReading field H\n" << endl;
+ volVectorField H(Hheader, mesh);
+ }
+
+ IOobject HdotGradHheader
(
- IOobject
- (
- "HdotGradH",
- runTime.timeName(),
- mesh,
- IOobject::READ_IF_PRESENT,
- IOobject::AUTO_WRITE
- ),
- mesh
+ "HdotGradH",
+ runTime.timeName(),
+ mesh,
+ IOobject::MUST_READ,
+ IOobject::AUTO_WRITE
);
+
+ if (HdotGradHheader.headerOk())
+ {
+
+ Info<< "\nReading field HdotGradH\n" << endl;
+ volVectorField HdotGradH(HdotGradHheader, mesh);
+ }
--
GitLab