[saralen]

Hi all,


I'm changing the code in such a way I can extract surface average of some variable at every deltat (no time step).

I don't want to use any kind of postProcess function but directly change the main code.



I'd like do something like that:

forAll(mesh.C().boundaryField(), patchi){

const scalarField& Ap = mesh.magSf().boundaryField()[patchi];
...
}

but referring in general to a surface at a fixed z, not only to boundaryField(). I want to replace this cycle with a cycle over height.


What I can use instead of .boundaryField() ?
What field of 'mesh' I have to access to?
Where I can find more detail about OpenFOAM classes/structure?


Hope to have been clear,
Sara