diff --git a/applications/utilities/parallelProcessing/reconstructPar/reconstructPar.C b/applications/utilities/parallelProcessing/reconstructPar/reconstructPar.C index 6b6e3bb535842350bf4f59e72f94c535c7308e0a..d545d64233106d57415cbb6edd17f7dcacb2e86d 100644 --- a/applications/utilities/parallelProcessing/reconstructPar/reconstructPar.C +++ b/applications/utilities/parallelProcessing/reconstructPar/reconstructPar.C @@ -3,7 +3,7 @@ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | \\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation - \\/ M anipulation | Copyright (C) 2015 OpenCFD Ltd. + \\/ M anipulation | Copyright (C) 2015-2018 OpenCFD Ltd. ------------------------------------------------------------------------------- License This file is part of OpenFOAM. @@ -728,7 +728,7 @@ int main(int argc, char *argv[]) ) { Info << "Reconstructing FA fields" << nl << endl; - + faMesh aMesh(mesh); processorFaMeshes procFaMeshes(procMeshes.meshes()); @@ -748,7 +748,7 @@ int main(int argc, char *argv[]) .reconstructFaAreaFields<sphericalTensor>(objects); faReconstructor.reconstructFaAreaFields<symmTensor>(objects); faReconstructor.reconstructFaAreaFields<tensor>(objects); - + faReconstructor.reconstructFaEdgeFields<scalar>(objects); } else diff --git a/src/parallel/reconstruct/reconstruct/reconstructLagrangianPositions.C b/src/parallel/reconstruct/reconstruct/reconstructLagrangianPositions.C index c8e8be4c5e0be65801dac7f2393ba1d8c080c8bf..30622183ec0c1b462716955ab04e5229088780cf 100644 --- a/src/parallel/reconstruct/reconstruct/reconstructLagrangianPositions.C +++ b/src/parallel/reconstruct/reconstruct/reconstructLagrangianPositions.C @@ -3,7 +3,7 @@ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | \\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation - \\/ M anipulation | + \\/ M anipulation | Copyright (C) 2018 OpenCFD Ltd. ------------------------------------------------------------------------------- License This file is part of OpenFOAM. @@ -60,7 +60,7 @@ void Foam::reconstructLagrangianPositions // Inverting sign if necessary and subtracting 1 from // faceProcAddressing - label mappedTetFace = mag(faceMap[ppi.tetFace()]) - 1; + const label mappedTetFace = mag(faceMap[ppi.tetFace()]) - 1; lagrangianPositions.append ( @@ -77,6 +77,16 @@ void Foam::reconstructLagrangianPositions } IOPosition<Cloud<passiveParticle>>(lagrangianPositions).write(); + + // Force writing of "positions" too, if specified via the InfoSwitch + if (particle::writeLagrangianPositions) + { + IOPosition<Cloud<passiveParticle>> + ( + lagrangianPositions, + IOPosition<Cloud<passiveParticle>>::geometryType::POSITIONS + ).write(); + } }