CFD Online Discussion Forums - runtime cuttingplanes in simpleFoam

Hello,

I'm trying to use sampling with cuttingplanes during run time.

For this, I modified the simpleFoam.C file by including some code taken from sample utility in the main while loop :

Code:

IOsampledSurfaces sSurfs

    (

        sampledSurfaces::typeName,

        mesh,

        "sampleDict",

        IOobject::MUST_READ,

        true

    );



        // Handle geometry/topology changes

        //polyMesh::readUpdateState state = mesh.readUpdate();



        //sSurfs.readUpdate(state);



        sSurfs.write();



        Info<< "write cuttingplanes at Time = " << runTime.timeName() << nl << endl;

All this works well, fields on cutting planes are written for each time step with the information given in sampleDict.

But now, I would like to use the fields on cuttingplane to compute statistics (mean of field, standard deviation on the cutting plane,...).

Can somenone help me to replace the sSurfs.write() statement with some code that could allow me to retrieve the fields value and to operate on them ?

Thanks