[TomWaits]

I have the following code in my fvOptions to write "USource" to a file. This works, however it writes "USource" for every cell in my domain, which is very slow. I am only interested in writing "USource" for the CellSet I define in topoSetDict. Could someone please help with a way to do this?

    
    if (mesh_.time().outputTime())
    {   
        volVectorField momentumSourceField
        (
            IOobject
            (
                "momentumSourceField",
                mesh_.time().timeName(),
                mesh_,
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            mesh_,
            vector::zero
        );        
        forAll(cells, i)
        {
            momentumSourceField[cells[i]] = Usource[cells[i]];
            momentumSourceField.write();
        }

Many thanks in advance,

Tom Waits

[zhangyan]

Are you looking for this?
https://github.com/StachuraMichal/Op...sampleCellZone

[TomWaits]

Originally Posted by zhangyan

Are you looking for this?
https://github.com/StachuraMichal/Op...sampleCellZone

Thank you for your response. I do not think so (although I am new to OpenFoam so I am unsure). "USource" is not a velocity or pressure that I am sampling in post, it is a variable that is calculated during fvOptions (see actuationDiskSourceTemplates.C in the source code).

I was hoping for a minor modification to the code I posted to loop over the cells defined in toposet as opposed to the entire grid.

As an aside, when I try to compile the code you sent, I get the following error:

cellZoneSet/cellZoneSet.C:139:20: error: ‘const class Foam::dictionary’ has no member named ‘get’; did you mean ‘set’?
     zoneName_(dict.get<word>("cellZone"))
                    ^~~
                    set
cellZoneSet/cellZoneSet.C:139:28: error: expected primary-expression before ‘>’ token
     zoneName_(dict.get<word>("cellZone"))

Any ideas as to why?

Regards,

Tom Waits