Hi I am using the subroutine DISP.
I defined a variable P used to determine when to fix the boundary.
if (KSTEP .EQ. 1) then
IF((COORDS(3) .EQ. 0.0000) .AND. (JDOF .EQ. 1) .AND. (P .GE. 0.10)) then
U(1)=0.0000000000
!write(*,*) 'P_conver',P_conver
END IF
IF((COORDS(3) .EQ. 0.0000)...