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

!-------------------------------------------------------------
! Store Interface tractions into Elres_t
!-------------------------------------------------------------
DO nd=1, NdofC
DO n=1, NbelC(nr)
ne=ListEC(nr,n)
DO m=1, Ndofe
IF (nd == LdesteR(ne,m)) THEN
Elres_t(ne,m)= Elres_t(ne,m) + tcxf(nd)
END IF
END DO
END DO
END DO
!---------------------------------------------------------
! Store Results of free nodes into Elres_u or Elres_t
!---------------------------------------------------------
DO nd=NdofC+1, NdofR
DO n=1, NbelF(nr)
ne=ListEF(nr,n)
DO m=1, Ndofe
IF (nd == LdesteR(ne,m)) THEN
IF (Ncode(nr,nd) == 0) THEN
Elres_u(ne,m)= Elres_u(ne,m) + tcxf(nd)
ELSE IF (Bcode(ne,m) == 1) THEN
Elres_t(ne,m)= Elres_t(ne,m) + tcxf(nd)
END IF
END IF
END DO
END DO
END DO
DEALLOCATE (tc,xf,tcxf,ucr)
END DO &
Regions_loop_3
!--------------------------
! Print out results
!--------------------------
CLOSE (UNIT=2)
OPEN (UNIT=2,FILE= 'BERESULTS',FORM='FORMATTED')
Elements_all:&
DO nel=1,Maxe
WRITE (2,*) ' Results, Element ',nel
WRITE (2,*) 'u=' , (Elres_u(nel,m), m=1,Ndofe)
WRITE (2,*) 't=' , (Elres_t(nel,m), m=1,Ndofe)
END DO &
Elements_all
END PROGRAM General_purpose_MRBEM