1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | foam-extend: Open Source CFD
4 \\ / O peration | Version: 3.2
5 \\ / A nd | Web: http://www.foam-extend.org
6 \\/ M anipulation | For copyright notice see file Copyright
7 -------------------------------------------------------------------------------
9 This file is part of foam-extend.
11 foam-extend is free software: you can redistribute it and/or modify it
12 under the terms of the GNU General Public License as published by the
13 Free Software Foundation, either version 3 of the License, or (at your
14 option) any later version.
16 foam-extend is distributed in the hope that it will be useful, but
17 WITHOUT ANY WARRANTY; without even the implied warranty of
18 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
19 General Public License for more details.
21 You should have received a copy of the GNU General Public License
22 along with foam-extend. If not, see <http://www.gnu.org/licenses/>.
24 \*---------------------------------------------------------------------------*/
27 void Foam::calcTypes::addSubtract::writeAddSubtractField
29 const IOobject& baseHeader,
30 const IOobject& addHeader,
35 typedef GeometricField<Type, fvPatchField, volMesh> fieldType;
39 baseHeader.headerClassName() == fieldType::typeName
40 && baseHeader.headerClassName() == addHeader.headerClassName()
43 if (resultName_ == "")
47 resultName_ = baseHeader.name() + "_add_" + addHeader.name();
51 resultName_ = baseHeader.name() + "_subtract_"
56 Info<< " Reading " << baseHeader.name() << endl;
57 fieldType baseField(baseHeader, mesh);
59 Info<< " Reading " << addHeader.name() << endl;
60 fieldType addField(addHeader, mesh);
62 if (baseField.dimensions() == addField.dimensions())
64 Info<< " Calculating " << resultName_ << endl;
71 mesh.time().timeName(),
75 calcMode_ == ADD ? baseField + addField : baseField - addField
81 Info<< " Cannot calculate " << resultName_ << nl
82 << " - inconsistent dimensions: "
83 << baseField.dimensions() << " - " << addField.dimensions()
92 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //