[heksel8i]

Hey!

I would like to create a postprocessing script for lagrangian data and I have been trying to look for examples (like particleTracks or lagrangianMap) but they use the PassiveParticle approach which won't reveal example diameter, age etc. information. So this far I have implemented the following and the question is :

How can I access the parcels of the SparyCloud, which I have?

Code this far:

 Info<< "Read lagrangian data"<< endl;


        fileNameList cloudDirs
        (
            readDir
            (
                runTime.timePath()/cloud::prefix,
                fileName::DIRECTORY
            )
        );

        forAll(cloudDirs, cloudI)
        {
            // Search for list of lagrangian objects for this time
            IOobjectList objects
            (
                mesh,
                runTime.timeName(),
                cloud::prefix/cloudDirs[cloudI]
            );

            //This is the way to access the IOobject
            IOobject*  diametersOwn = objects.lookup("d");

     
/*
// SHOULD IT BE SOMETHING LIKE THIS???

            Cloud<basicSprayCloud> Parcels
            (
                mesh,
                cloudDirs[cloudI],
                false
            );
            Info << Parcels.d() << endl;
            
*/

Thanks in advance.

[heksel8i]

To reply myself and others having this problem in the future: One can proceed in the following way:

        Info<< "Read lagrangian data"<< endl;

        fileNameList cloudDirs
        (
            readDir
            (
                runTime.timePath()/cloud::prefix,
                fileName::DIRECTORY
            )
        );


        forAll(cloudDirs, cloudI)
        {
            // Search for list of lagrangian objects for this time
            IOobjectList objects
            (
                mesh,
                runTime.timeName(),
                cloud::prefix/cloudDirs[cloudI]
            );

            //accessing to the particle data, example diameter
            const scalarField& d = readParticleField<scalar>("d", objects);