Merge branch 'master' of /home/noisy2/OpenFOAM/OpenFOAM-dev

This commit is contained in:
mattijs
2008-06-23 16:08:06 +01:00
5 changed files with 178 additions and 73 deletions

View File

@ -405,64 +405,71 @@ bool doCommand
backup(mesh, setName, currentSet, setName + "_old"); backup(mesh, setName, currentSet, setName + "_old");
} }
if (action == topoSetSource::CLEAR) switch (action)
{ {
// Already handled above by not reading case topoSetSource::CLEAR:
}
else if (action == topoSetSource::INVERT)
{
currentSet.invert(currentSet.maxSize(mesh));
}
else if (action == topoSetSource::LIST)
{
currentSet.writeDebug(Pout, mesh, 100);
Pout<< endl;
}
else if (action == topoSetSource::SUBSET)
{
if (is >> sourceType)
{ {
autoPtr<topoSetSource> setSource // Already handled above by not reading
( break;
topoSetSource::New
(
sourceType,
mesh,
is
)
);
// Backup current set.
topoSet oldSet
(
mesh,
currentSet.name() + "_old2",
currentSet
);
currentSet.clear();
currentSet.resize(oldSet.size());
setSource().applyToSet(topoSetSource::NEW, currentSet);
// Combine new value of currentSet with old one.
currentSet.subset(oldSet);
} }
} case topoSetSource::INVERT:
else
{
if (is >> sourceType)
{ {
autoPtr<topoSetSource> setSource currentSet.invert(currentSet.maxSize(mesh));
( break;
topoSetSource::New }
case topoSetSource::LIST:
{
currentSet.writeDebug(Pout, mesh, 100);
Pout<< endl;
break;
}
case topoSetSource::SUBSET:
{
if (is >> sourceType)
{
autoPtr<topoSetSource> setSource
( (
sourceType, topoSetSource::New
mesh, (
is sourceType,
) mesh,
); is
)
);
setSource().applyToSet(action, currentSet); // Backup current set.
topoSet oldSet
(
mesh,
currentSet.name() + "_old2",
currentSet
);
currentSet.clear();
currentSet.resize(oldSet.size());
setSource().applyToSet(topoSetSource::NEW, currentSet);
// Combine new value of currentSet with old one.
currentSet.subset(oldSet);
}
break;
}
default:
{
if (is >> sourceType)
{
autoPtr<topoSetSource> setSource
(
topoSetSource::New
(
sourceType,
mesh,
is
)
);
setSource().applyToSet(action, currentSet);
}
} }
} }
@ -532,7 +539,7 @@ bool doCommand
} }
} }
return error; return !error;
} }
@ -612,31 +619,36 @@ commandStatus parseType
switch(stat) switch(stat)
{ {
case polyMesh::UNCHANGED: case polyMesh::UNCHANGED:
{
Pout<< " mesh not changed." << endl; Pout<< " mesh not changed." << endl;
break; break;
}
case polyMesh::POINTS_MOVED: case polyMesh::POINTS_MOVED:
{
Pout<< " points moved; topology unchanged." << endl; Pout<< " points moved; topology unchanged." << endl;
break; break;
}
case polyMesh::TOPO_CHANGE: case polyMesh::TOPO_CHANGE:
{
Pout<< " topology changed; patches unchanged." << nl Pout<< " topology changed; patches unchanged." << nl
<< " "; << " ";
printMesh(runTime, mesh); printMesh(runTime, mesh);
break;
break; }
case polyMesh::TOPO_PATCH_CHANGE: case polyMesh::TOPO_PATCH_CHANGE:
{
Pout<< " topology changed and patches changed." << nl Pout<< " topology changed and patches changed." << nl
<< " "; << " ";
printMesh(runTime, mesh); printMesh(runTime, mesh);
break; break;
}
default: default:
{
FatalErrorIn("parseType") << "Illegal mesh update state " FatalErrorIn("parseType") << "Illegal mesh update state "
<< stat << abort(FatalError); << stat << abort(FatalError);
break; break;
}
} }
return INVALID; return INVALID;

85
bin/upgradeTurbulenceProperties Executable file
View File

@ -0,0 +1,85 @@
#!/bin/sh
#------------------------------------------------------------------------------
# ========= |
# \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
# \\ / O peration |
# \\ / A nd | Copyright (C) 1991-2007 OpenCFD Ltd.
# \\/ M anipulation |
#------------------------------------------------------------------------------
# License
# This file is part of OpenFOAM.
#
# OpenFOAM is free software; you can redistribute it and/or modify it
# under the terms of the GNU General Public License as published by the
# Free Software Foundation; either version 2 of the License, or (at your
# option) any later version.
#
# OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
# FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
# for more details.
#
# You should have received a copy of the GNU General Public License
# along with OpenFOAM; if not, write to the Free Software Foundation,
# Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
#
# Script
# upgradeTurbulenceProperties
#
# Description
# Upgrade the turbulenceProperties dictionary to the new format employed
# in OpenFOAM version 1.5
# - RAS turbulence models now defined by the RASProperties dictionary,
# and RASModel keyword, and
# - LES turbulence models now defined by the LESProperties dictionary,
# and LESModel keyword.
#
#------------------------------------------------------------------------------
printUsage()
{
echo "Usage: `basename $0` <turbulenceProperties>"
echo " Where <turbulenceProperties> is the full path to the"
echo " turbulenceProperties dictionary"
}
convertDict()
{
echo " Identified $1 turbulence model"
sed -e "s/turbulenceProperties/$1Properties/" \
-e "s/$2/$1Model/" \
-e "s/[a-zA-Z0-9]* [ ]*\[[0-9 ]*\]//" \
$3 > "$outputPath/$1Properties"
echo " written $1Properties to $outputPath/"
}
outputPath=`dirname $1`
if [ $# -ne 1 ]; then
printUsage
exit 1
elif [ ! -e $1 ]; then
echo " Error: file $1 does not exist"
echo ""
printUsage
exit 1
fi
# Identify type of turbulence model
RAS=`grep turbulenceModel $1`
LES=`grep LESmodel $1`
if [ -n "$RAS" ]; then
convertDict "RAS" "turbulenceModel" $1
elif [ -n "$LES" ]; then
convertDict "LES" "LESmodel" $1
else
echo "Unable to determine turbulence model type - nothing changed"
exit 1
fi
echo "done."
exit 0

View File

@ -111,6 +111,8 @@ void Foam::calcTypes::components::calc
mesh, mesh,
processed processed
); );
writeComponentFields<symmTensor>(fieldHeader, mesh, processed);
writeComponentFields<tensor>(fieldHeader, mesh, processed);
if (!processed) if (!processed)
{ {

View File

@ -25,13 +25,14 @@ FoamFile
interpolationScheme cellPoint; interpolationScheme cellPoint;
writeFormat raw; setFormat raw;
sampleSets sets
( (
uniform data
{ {
name data; type uniform;
axis x; axis x;
start (-4.995 0 0); start (-4.995 0 0);
end (4.995 0 0); end (4.995 0 0);
@ -39,6 +40,9 @@ sampleSets
} }
); );
surfaces
();
fields fields
( (
T T

View File

@ -25,13 +25,13 @@ FoamFile
interpolationScheme cellPoint; interpolationScheme cellPoint;
writeFormat raw; setFormat raw;
sampleSets sets
( (
uniform leftPatch
{ {
name leftPatch; type uniform;
axis y; axis y;
start (0 0.5 0.25); start (0 0.5 0.25);
end (0 2 0.25); end (0 2 0.25);
@ -39,10 +39,12 @@ sampleSets
} }
); );
fields surfaces
();
fields
( (
sigmaxx sigmaxx
); );
// ************************************************************************* // // ************************************************************************* //