I was debuging the alelav.f90 subroutine and I found that may be there is a bug in the following line, because to get the type of boundary condition at the boundary face ifac is used itypfb(ifac).
Code: Select all
if (ialtyb(ifac) .eq. ifresf) then
Code: Select all
if (itypfb(ifac).eq.ifresf) then
I attach the modified alelav.f90, paste it in the SRC directory and check if its works. May be the developers can check if it is correct.
Regards,
Luciano