[ljjouc]

Dear Foamers,

I really encounter an urgent problem in creating my own solver. Here is the problem:

I have obtained a VolScalarField rou(x,y,z,t), now I want to integrate it from water surface z=h to a water depth z so as to get a hydrostatic pressure head r(x,y,z,t), and the corresponding formula is stored in the attached figure. But, I just do not know how to code program to implement this integration?

Shall anyone kindly give me some advice or tell me how I should program by means of the OF programming language? The solving of this problem involves my Phd project and it is very important. Thanks in advance!