int main(int argc, char *argv[])
{
argList::noParallel();
timeSelector::addOptions();
# include "setRootCase.H"
# include "createTime.H"
// Get times list
instantList timeDirs = timeSelector::select0(runTime, args);
# include "createMesh.H"
# include "readTransportProperties.H"
const word& gFormat = runTime.graphFormat();
// Setup channel indexing for averaging over channel down to a line
IOdictionary channelDict
(
IOobject
(
"postSedimentDict",
mesh.time().constant(),
mesh,
IOobject::MUST_READ_IF_MODIFIED,
IOobject::NO_WRITE
)
);
channelIndex channelIndexing(mesh, channelDict);
// For each time step read all fields
forAll(timeDirs, timeI)
{
runTime.setTime(timeDirs[timeI], timeI);
Info<< "Collapsing fields for time " << runTime.timeName() << endl;
# include "readFields.H"
# include "calculateFields.H"
// Average fields over channel down to a line
# include "collapse.H"
}
int main(int argc, char *argv[])
{
argList::noParallel();
timeSelector::addOptions();
# include "setRootCase.H"
# include "createTime.H"
// Get times list
instantList timeDirs = timeSelector::select0(runTime, args);
# include "createMesh.H"
# include "readTransportProperties.H"
const word& gFormat = runTime.graphFormat();
// Setup channel indexing for averaging over channel down to a line
IOdictionary channelDict
(
IOobject
(
"postChannelExtDict",
mesh.time().constant(),
mesh,
IOobject::MUST_READ_IF_MODIFIED,
IOobject::NO_WRITE
)
);
channelIndex channelIndexing(mesh, channelDict);
List<Triple<word> > fields(channelDict.lookup("fieldsAndIdentifiers"));
bool backwardsCompatibility(channelDict.lookupOrDefault("backwardsCompatibility", true));
// For each time step read all fields
forAll(timeDirs, timeI)
{
runTime.setTime(timeDirs[timeI], timeI);
Info << "Processing fields for time " << runTime.timeName() << endl;
Info << endl;
fileName path(runTime.rootPath()/runTime.caseName()/"graphs"/runTime.timeName());
mkDir(path);
forAll(fields, I)
{
const word& fieldName = fields[I].first();
const word& identifierName = fields[I].second();
const word& moment = fields[I].third();
IOobject fieldHeader
(
fieldName,
runTime.timeName(),
mesh,
IOobject::MUST_READ
);
if (!fieldHeader.headerOk())
{
Info << endl;
Info<< "No " << fieldName <<" field" << endl;
Info << endl;
continue;
}
Info << endl;
Info << fieldName << endl;
if(fieldHeader.headerClassName() == volScalarField::typeName)
{
Info << endl;
Info << " Reading field" << endl;
volScalarField field(fieldHeader, mesh);
Info << " Collapsing field" << endl;
Info << endl;
processField(field, channelIndexing, identifierName, path, gFormat, moment);
}
if(fieldHeader.headerClassName() == volVectorField::typeName)
{
Info << endl;
Info << " Reading field" << endl;
volVectorField field(fieldHeader, mesh);
Info << " Collapsing field" << endl;
Info << endl;
if((fieldHeader.name() == "U" || fieldHeader.name() == "UMean") && backwardsCompatibility)
{
processField(field.component(vector::X)(), channelIndexing, identifierName, path, gFormat, "mean");
}
else
{
processField(field.component(vector::X)(), channelIndexing, identifierName+"_x", path, gFormat, moment);
processField(field.component(vector::Y)(), channelIndexing, identifierName+"_y", path, gFormat, moment);
processField(field.component(vector::Z)(), channelIndexing, identifierName+"_z", path, gFormat, moment);
}
}
if(fieldHeader.headerClassName() == volSymmTensorField::typeName)
{
Info << endl;
Info << " Reading field" << endl;
volSymmTensorField field(fieldHeader, mesh);
Info << " Collapsing field" << endl;
Info << endl;
if((fieldHeader.name() == "UPrime2Mean") && backwardsCompatibility)
{
processField(field.component(symmTensor::XX)(), channelIndexing, "u", path, gFormat, "rms");
processField(field.component(symmTensor::YY)(), channelIndexing, "v", path, gFormat, "rms");
processField(field.component(symmTensor::ZZ)(), channelIndexing, "w", path, gFormat, "rms");
processField(field.component(symmTensor::XY)(), channelIndexing, "uv", path, gFormat, "mean", true);
}
else
{
processField(field.component(symmTensor::XX)(), channelIndexing, identifierName+"_xx", path, gFormat, moment);
processField(field.component(symmTensor::XY)(), channelIndexing, identifierName+"_xy", path, gFormat, moment, true);
processField(field.component(symmTensor::XZ)(), channelIndexing, identifierName+"_xz", path, gFormat, moment, true);
processField(field.component(symmTensor::YY)(), channelIndexing, identifierName+"_yy", path, gFormat, moment);
processField(field.component(symmTensor::YZ)(), channelIndexing, identifierName+"_yz", path, gFormat, moment, true);
processField(field.component(symmTensor::ZZ)(), channelIndexing, identifierName+"_zz", path, gFormat, moment);
}
}
if(fieldHeader.headerClassName() == volTensorField::typeName)
{
Info << endl;
Info << " Reading field" << endl;
volTensorField field(fieldHeader, mesh);
Info << " Collapsing field" << endl;
Info << endl;
processField(field.component(tensor::XX)(), channelIndexing, identifierName+"_xx", path, gFormat, moment);
processField(field.component(tensor::XY)(), channelIndexing, identifierName+"_xy", path, gFormat, moment, true);
processField(field.component(tensor::XZ)(), channelIndexing, identifierName+"_xz", path, gFormat, moment, true);
processField(field.component(tensor::YX)(), channelIndexing, identifierName+"_yx", path, gFormat, moment, true);
processField(field.component(tensor::YY)(), channelIndexing, identifierName+"_yy", path, gFormat, moment);
processField(field.component(tensor::YZ)(), channelIndexing, identifierName+"_yz", path, gFormat, moment, true);
processField(field.component(tensor::XZ)(), channelIndexing, identifierName+"_xz", path, gFormat, moment, true);
processField(field.component(tensor::YZ)(), channelIndexing, identifierName+"_yz", path, gFormat, moment, true);
processField(field.component(tensor::ZZ)(), channelIndexing, identifierName+"_zz", path, gFormat, moment);
}
}
}
