findCloudFields.H

Go to the documentation of this file.

// check all time directories for the following:
 
// The fields for each cloud:
HashTable<HashTable<word>> cloudFields;
 
// Identify if lagrangian data exist at any time step.
if (timeDirs.size() && doLagrangian)
{
    const fileName& baseDir = mesh.time().path();
    const fileName cloudPrefix(regionPrefix/cloud::prefix);
 
    Info<< "Searching for lagrangian ... " << flush;
 
    for (const instant& inst : timeDirs)
    {
        const word& timeName = inst.name();
 
        // DO NOT USE -->>  runTime.setTime(timeDirs[timeI], timeI);  <<--
        // It incurs a large overhead when done so frequently.
 
        fileNameList cloudDirs
        (
            readDir
            (
                baseDir/timeName/cloudPrefix,
                fileName::DIRECTORY
            )
        );
 
        for (fileName& cloudDir : cloudDirs)
        {
            const word cloudName(std::move(cloudDir));
 
            IOobjectList cloudObjs
            (
                mesh,
                timeName,
                cloudPrefix/cloudName
            );
 
            // Clouds require "coordinates".
            // The "positions" are for v1706 and lower.
            // - detect and remove since these are treated specially
 
            bool isCloud = false;
            if (cloudObjs.erase("coordinates"))
            {
                isCloud = true;
            }
            if (cloudObjs.erase("positions"))
            {
                isCloud = true;
            }
 
            if (isCloud)
            {
                // Save the cloud fields on a per cloud basis
                auto& fieldsPerCloud = cloudFields(cloudName);
 
                forAllConstIters(cloudObjs, fieldIter)
                {
                    const IOobject* io = *fieldIter;
 
                    // Field name/type
                    fieldsPerCloud.insert(io->name(), io->headerClassName());
                }
            }
        }
    }
 
    if (Pstream::parRun())
    {
        Pstream::mapCombineGather(cloudFields, HashTableOps::plusEqOp<word>());
        Pstream::mapCombineScatter(cloudFields);
    }
 
    if (cloudFields.empty())
    {
        Info<< "none detected." << endl;
    }
}
 
 
// Sorted list of cloud names
const wordList cloudNames(cloudFields.sortedToc());
 
if (cloudNames.size())
{
    // Complete the echo information - as flatOutput
    cloudNames.writeList(Info) << endl;
}
 
 
// ************************************************************************* //