[lzhou]

Hi,

I am using Finite area method to simulate the scour on river bed.

When calculating the divergence of an "areaVectorField" Qsa (fac::div(Qsa)), I found out that the program suppose the Qsa value at the boundary of the finite area mesh is (0 0 0), even though I set the boundary condition of Qsa as zeroGradient in file 0/Qsa. So the divergence calculated is not correct at the boundary of finite area mesh.

I would like to ask how to solve this problem ? Or is there a way to get the label of the face cell at the boundary of the finite area mesh and correct it ? Thanks !

Lu ZHOU