Assembly and Solution - The Boundary Element Method with Programming - page 189

SUBROUTINE Geomin(Nodes,Maxe,xp,Inci,Nodel,Cdim)
!------------------------------------
! Inputs mesh geometry
!-------------------------------------
INTEGER, INTENT(IN) :: Nodes ! Number of nodes
INTEGER, INTENT(IN) :: Maxe ! Number of elements
INTEGER, INTENT(IN) :: Nodel ! Number of Nodes of elements
INTEGER, INTENT(IN) :: Cdim ! Cartesian Dimension
REAL, INTENT(OUT) :: xP(:,:) ! Node co-ordinates
REAL :: xmax(Cdim),xmin(Cdim),delta_x(Cdim)
INTEGER, INTENT(OUT):: Inci(:,:) ! Element incidences
INTEGER :: Node,Nel,M,n
!-------------------------------------------------------
! Read Node Co-ordinates from Inputfile
!-------------------------------------------------------
DO Node=1,Nodes
READ( 1,*) (xP(M,Node),M=1,Cdim)
WRITE (2,'(A5,I5,A8,3F8.2)') 'Node ',Node,&
' Coor ',(xP(M,Node),M=1,Cdim)
END DO
!-------------------------------------------------------
! Read Incidences from Inputfile
!-------------------------------------------------------
WRITE (2,*)''
WRITE (2,*)'Incidences: '
WRITE (2,*)''
Elements_1:&
DO Nel=1,Maxe
READ(1,*) (Inci(Nel,n),n=1,Nodel)
WRITE(2,'(A3,I5,A8,24I5)')'EL ',Nel,' Inci ',Inci(Nel,:)
END DO &
Elements_1
RETURN
END SUBROUTINE Geomin
SUBROUTINE BCInput(Elres_u,Elres_t,Bcode,nodel,ndofe,ndof)
!------------------------------------------
! Reads boundary conditions
!-------------------------------------------
REAL,INTENT(OUT) :: Elres_u(:,:) ! Element results , u
REAL,INTENT(OUT) :: Elres_t(:,:) ! Element results , t
INTEGER,INTENT(OUT) :: BCode(:,:) ! Element BC´s
INTEGER,INTENT(IN) :: nodel ! Nodes per element
INTEGER,INTENT(IN) :: ndofe ! D.o.F. per Element
INTEGER,INTENT(IN) :: ndof ! D.o.F per Node
INTEGER :: NE_u,NE_t
WRITE (2,*)''
WRITE (2,*)'Elements with Dirichlet BC´s: '
WRITE (2,*)''
Elres_u(:,:)=0 ! Default prescribed values for u = 0.0