Multiple regions - The Boundary Element Method with Programming - page 309

DO nr=1,Nregs
node= 0
Elements_of_region: &
DO nb=1,Nbel(nr)
ne= ListR(nr,nb)
Interface_elements: &
IF (Bcode(ne,1) == 2) THEN
NbelC(nr)= NbelC(nr) + 1
ListEC(nr,NbelC(nr))= ne
Nodes_of_Elem: &
DO n=1,Nodel
! check if node has allready been entered
is=0
DO no=1,node
IF(InciR(nr,no) == Incie(ne,n)) THEN
is= 1
EXIT
END IF
END DO
IF (is == 0) THEN
node=node+1
InciR(nr,node)= Incie(ne,n)
IncieR(ne,n)= node
ELSE
IncieR(ne,n)= no
END IF
END DO &
Nodes_of_Elem
END IF &
Interface_elements
END DO &
Elements_of_region
NodeC(nr)= Node ! No of interface nodes of Region nr
NdofC= NodeC(nr)*Ndof ! D.o.F. at interface of Region nr
Elements_of_region1: &
DO nb=1,Nbel(nr)
ne= ListR(nr,nb)
Free_elements: &
IF (Bcode(ne,1) /= 2) THEN
NbelF(nr)= NbelF(nr) + 1
ListEF(nr,NbelF(nr))= ne
Nodes_of_Elem1: &
DO n=1,Nodel
is=0
DO no=1,node
IF (InciR(nr,no) == Incie(ne,n)) THEN
is= 1
EXIT
END IF
END DO
IF (is == 0) THEN