Postprocessing - The Boundary Element Method with Programming - page 245

points requiring the same number of integration points. In this case the number of
computations of the Jacobian can be reduced significantly. Using table 6.1 and element
subdivision it will be found later that the internal points may be placed quite close to the
boundary.
PROGRAM Post_processor
!------------------------------------------------
! General purpose Postprocessor
! for computing results at boundary and interior points
!------------------------------------------------------
USE Utility_lib; USE Elast_lib; USE Laplace_lib
USE Integration_lib
USE Postproc_lib
IMPLICIT NONE
INTEGER, ALLOCATABLE :: Inci(:) ! Incidences (one elem.)
INTEGER, ALLOCATABLE :: Incie(:,:) ! Incidences (all elem.)
INTEGER, ALLOCATABLE :: Ldest(:) ! Destinations (one elem.)
REAL, ALLOCATABLE :: Elcor(:,:) ! Element coordinates
REAL, ALLOCATABLE :: El_u(:,:,:)!
REAL, ALLOCATABLE :: El_t(:,:,:)! Results of System
REAL, ALLOCATABLE :: El_ue(:,:) ! Diplacements of Element
REAL, ALLOCATABLE :: El_te(:,:) ! Traction of Element
REAL, ALLOCATABLE :: Disp(:) ! Diplacement results Node
REAL, ALLOCATABLE :: Trac(:) ! Traction results of Node
REAL, ALLOCATABLE :: El_trac(:) ! Traction results Element
REAL, ALLOCATABLE :: El_disp(:) ! Displacement of Element
REAL, ALLOCATABLE :: xP(:,:) ! Node co-ordinates of BE
REAL, ALLOCATABLE :: xPnt(:) ! Co-ordinates of int. point
REAL, ALLOCATABLE :: Ni(:),GCcor(:),dxr(:),Vnorm(:)
CHARACTER (LEN=80) :: Title
REAL :: Elengx,Elenge,Rmin,Glcorx(8),Wix(8),Glcore(8),Wie(8)
REAL :: Jac
REAL :: Xsi1,Xsi2,Eta1,Eta2,RJacB,RonL
REAL, ALLOCATABLE :: Flow(:),Stress(:)! Results for bound.Point
REAL, ALLOCATABLE :: uPnt(:),SPnt(:) ! Results for int Point
REAL, ALLOCATABLE :: TU(:,:),UU(:,:) ! Kernels for u
REAL, ALLOCATABLE :: TS(:,:),US(:,:) ! Kernels for q,s
REAL, ALLOCATABLE :: Fac(:),Fac_nod(:,:) ! Fact. for symmetry
INTEGER :: Cdim,Node,M,N,Istat,Nodel,Nel,Ndof,Cod,Nreg
INTEGER :: Ltyp,Nodes,Maxe,Ndofe,Ndofs,Ncol,ndg,ldim
INTEGER :: nod,nd,Nstres,Nsym,Isym,nsy,IPS,Nan,Nen,Ios,dofa,dofe
INTEGER :: Mi,Ki,K,I,NDIVX,NDIVSX,NDIVE,NDIVSE,MAXDIVS
REAL : : Con,E,ny,Fact,G,C2,C3,C5,C6,C7
REAL :: xsi,eta,Weit,R,Rlim(2)
OPEN (UNIT=1,FILE='INPUT',FORM='FORMATTED')
OPEN (UNIT=2,FILE='OUTPUT',FORM='FORMATTED')
Call Jobin(Title,Cdim,Ndof,IPS,Nreg,Ltyp,Con,E,ny,&
Isym,nodel,nodes,maxe)
Ndofe= nodel*ndof
ldim= Cdim-1