I want to divide my domain into 5 layers and calculate the average velocity at each layer. Can anyone suggest me how to do it ? I am thinking about using computeAverage(computeVelocityNorm) , but I am not sure how to divide into 5 layers.
is your domain 2d or 3d? what exactly do you mean by layer?
my domain is 3D , my layer means region. I divided my domain into 5 rectangular regions and I would like to know the average velocity in that region.
I guess you can use
Array<T,nDim> computeAverage(MultiTensorField3D<T,nDim>& tensorField, Box3D domain)
dataProcessors/dataAnalysisWrapper3D.hh , where domain is each one of your 5 domains.
It would look like:
Array<T,3> avgVel = computeAverage(*computeVelocity(lattice), domain);