Postprocessing - The Boundary Element Method with Programming - page 250

END DO &
Gauss_points_xsi
Eta1= Eta2
END DO Subdivisions_eta
Xsi1= Xsi2
END DO Subdivisions_xsi
END DO Symmetry_loop1
END DO Element_loop1
WRITE (2,'(A,3F10.3)') ' u: ',uPnt
!------------------------------------------------------------
! Computation of Fluxes/Stresses at Points inside a region
!------------------------------------------------------------
SPnt= 0.0
Element_loop2: &
DO NEL= 1,MAXE
Symmetry_loop2: &
DO nsy=1,Nsym
Inci= Incie(nel,:)
Elcor= xp(:,Inci(:))
IF (ldim == 2) THEN
ELengx= Dist((Elcor(:,3)+Elcor(:,2))/2.&
,(Elcor(:,4)+Elcor(:,1))/2.,Cdim) ! Lxsi
ELenge= Dist((Elcor(:,2)+Elcor(:,1))/2.&
,(Elcor(:,3)+Elcor(:,4))/2.,Cdim) ! Leta
ELSE
Call Elength(Elengx,Elcor,nodel,ldim)
END IF
Ldest= 1
Fac= 1.0
El_ue(:,:)=El_u(Nel,:,:)
El_te(:,:)=El_t(Nel,:,:)
IF (Isym > 0) THEN
DO Nod=1,Nodel
dofa= (nod-1)*Ndof+1
dofe= dofa+Ndof-1
El_trac(dofa:dofe)= El_te(Nod,:)
EL_disp(dofa:dofe)= El_ue(Nod,:)
END DO
CALL Mirror(Isym,nsy,Nodes,Elcor,Fac,Inci&
,Ldest,El_trac,EL_disp &
,nodel,ndof,Cdim)
DO Nod=1,Nodel
dofa= (nod-1)*Ndof+1
dofe= dofa+Ndof-1
El_te(Nod,:)= El_trac(dofa:dofe)
El_ue(Nod,:)= El_disp(dofa:dofe)
Fac_nod(Nod,:)= Fac(dofa:dofe)
END DO
End IF
Rmin= Min_dist(Elcor,xPnt,Nodel,ldim,Inci)
Mi= Ngaus(Rmin/Elengx,Cdim,Rlim)