zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: ensight output failing with dimensioned field (fixes #244)

Browse Source 
- The new field needs initialization with a dimensioned<Type> not just
  the dimensionSet.

- The new field was also incorrectly being registered, which could
  cause issues later.
This commit is contained in:
Mark Olesen
2016-09-23 13:13:12 +02:00
parent 3c25f26aab
commit 06699f17f5
1 changed files with 7 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
applications/utilities/postProcessing/dataConversion/foamToEnsight/ensightField.C
View File

@ -81,7 +81,9 @@ volField
// Construct volField (with zeroGradient) from dimensioned field // Construct volField (with zeroGradient) from dimensioned field
IOobject io(df); IOobject io(df);
io.readOpt() = IOobject::NO_READ; io.readOpt() = IOobject::NO_READ;
io.writeOpt() = IOobject::NO_WRITE;
io.registerObject() = false;
tmp<GeometricField<Type, fvPatchField, volMesh>> tvf tmp<GeometricField<Type, fvPatchField, volMesh>> tvf
( (
@ -89,16 +91,17 @@ volField
( (
io, io,
df.mesh(), df.mesh(),
df.dimensions(), dimensioned<Type>("0", df.dimensions(), Zero),
zeroGradientFvPatchField<scalar>::typeName zeroGradientFvPatchField<Type>::typeName
) )
); );
tvf.ref().internalField() = df; tvf.ref().internalField() = df;
tvf.ref().correctBoundaryConditions(); tvf.ref().correctBoundaryConditions();
const GeometricField<Type, fvPatchField, volMesh>& vf = tvf();
if (meshSubsetter.hasSubMesh()) if (meshSubsetter.hasSubMesh())
{ {
const GeometricField<Type, fvPatchField, volMesh>& vf = tvf();
tmp<GeometricField<Type, fvPatchField, volMesh>> tfld tmp<GeometricField<Type, fvPatchField, volMesh>> tfld
( (
meshSubsetter.interpolate(vf) meshSubsetter.interpolate(vf)