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

added surfaceScalarField functionality

Browse Source
This commit is contained in:
andy
2008-05-16 15:39:14 +01:00
parent 9acc7a8b21
commit 31cea98bda
1 changed files with 35 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files

44
applications/utilities/postProcessing/patch/patchIntegrate/patchIntegrate.C
View File

@ -77,21 +77,45 @@ int main(int argc, char *argv[])
{ {
mesh.readUpdate(); mesh.readUpdate();
Info<< " Reading field " << fieldName << endl; label patchi = mesh.boundaryMesh().findPatchID(patchName);
if (patchi < 0)
{
FatalError
<< "Unable to find patch " << patchName << nl
<< exit(FatalError);
}
if (fieldHeader.headerClassName() == "volScalarField")
{
Info<< " Reading volScalarField " << fieldName << endl;
volScalarField field(fieldHeader, mesh); volScalarField field(fieldHeader, mesh);
label patchi = mesh.boundaryMesh().findPatchID(patchName); vector sumField = sum
if (patchi >= 0)
{
Info<< " Integral of " << fieldName << " over patch "
<< patchName << '[' << patchi << ']' << " = "
<< sum
( (
mesh.Sf().boundaryField()[patchi] mesh.Sf().boundaryField()[patchi]
*field.boundaryField()[patchi] *field.boundaryField()[patchi]
) );
<< endl;
Info<< " Integral of " << fieldName << " over patch "
<< patchName << '[' << patchi << ']' << " = "
<< sumField << nl;
}
else if (fieldHeader.headerClassName() == "surfaceScalarField")
{
Info<< " Reading surfaceScalarField " << fieldName << endl;
surfaceScalarField field(fieldHeader, mesh);
scalar sumField = sum(field.boundaryField()[patchi]);
Info<< " Integral of " << fieldName << " over patch "
<< patchName << '[' << patchi << ']' << " = "
<< sumField << nl;
}
else
{
FatalError
<< "Only possible to integrate volScalarFields "
<< "and surfaceScalarFields" << nl << exit(FatalError);
} }
} }
else else