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

BCode = 0 ! Default BC= Neumann Condition
READ (1,*)NE_u
IF( NE_u > 0) THEN
Elem_presc_displ: &
DO n=1,NE_u
READ (1,*) Nel,(Elres_u(Nel,m),m=1,Ndofe)
BCode(Nel,:)=1
WRITE (2,*)'Element ',Nel,' Prescribed values: '
Na= 1
Nodes: &
DO M= 1,Nodel
WRITE(2,*) Elres_u(Nel,na:na+ndof-1)
Na= na+Ndof
END DO &
Nodes
END DO &
Elem_presc_displ
END IF
WRITE(2,*)''
WRITE (2,*)'Elements with Neuman BC´s: '
WRITE (2,*)''
Elres_t(:,:)=0 ! Default prescribed values = 0.0
READ (1,*)NE_t
Elem_presc_trac: &
DO n=1,NE_t
READ (1,*) Nel,(Elres_t(Nel,m),m=1,Ndofe)
WRITE (2,*)'Element ',Nel,' Prescribed values: '
Na= 1
Nodes1: &
DO M= 1,Nodel
WRITE(2,*) Elres_t(Nel,na:na+ndof-1)
Na= na+Ndof
END DO &
Nodes1
END DO &
Elem_presc_trac
RETURN
END SUBROUTINE BCInput
7.4.1 User's manual
The input data which have to be supplied in file INPUT are described below. Free field
input is used, that is, numbers are separated by blanks. However, all numbers, including
zero entries must be specified.
The input is divided into two parts. First, general information about the problem is
read in, then the mesh geometry is specified. The problem may consist of linear and
quadratic elements, as shown in Figure 7.8. The sequence in which node numbers have
to be entered when specifying incidences is also shown. Note that this order determines
the direction of the outward normal, which has to point away from the material. For 3-D