Provide a dedicated subroutine to evaluate the integrated value of a field

[mathrack]

Several cases estimate the integrated value of a field on

• a face, for instance the inlet or the outlet
• a volume, for instance to estimate the total flow rate in the domain

In case of stretching, the integration is often approximate :

Incompact3d/src/Case-PTBL.f90

Line 1350 in 5565268

ut1=ut1+(yp(j+1)-yp(j))*(ux2(1,j+1,k)-half*(ux2(1,j+1,k)-ux2(1,j,k)))

It would be best to use the stretching (ppy) to integrate exactly fields :

Incompact3d/src/Case-Channel.f90

Line 240 in 5565268

ub = ub + ux(i,jloc,k) / ppy(j)

It might be best to provide a dedicated module with subroutines to integrate fields on faces and volumes to avoid code duplication.