Postprocessing - The Boundary Element Method with Programming - page 252

END IF
Node_loop2:&
DO Node=1,Nodel
Disp= El_ue(Node,:)* Fac_nod(Node,:)
Trac= El_te(Node,:)* Fac_nod(Node,:)
SPnt= SPnt + (MATMUL(TS,Trac)- &
MATMUL(US,Disp))* Ni(Node)* Fact
END DO Node_loop2
END DO Gauss_points_eta2
END DO Gauss_points_xsi2
Eta1= Eta2
END DO Subdivisions_eta1
Xsi1= Xsi2
END DO Subdivisions_xsi1
END DO Symmetry_loop2
END DO Element_loop2
IF (Ndof == 1) THEN
WRITE (2,'(A,6F10.3)') ' Flux: ',SPnt
ELSE
IF (CDIM == 2) THEN
IF (IPS==1) THEN
SPnt(4)=SPnt(3)
SPnt(3)= ny*(SPnt(1)+SPnt(2))
ELSE
SPnt(4)=SPnt(3)
SPnt(3)= 0
END IF
END IF
WRITE(2,'(A,6F10.3)') ' Stress: ',SPnt
END IF
END DO Internal_points
END PROGRAM Post_processor

9.4.1 Input specification

INPUT DATA SPECIFICATION FOR Postprocessor

1.0 Boundary results

Nan, Nen First element and last element on

which boundary results are to be computed

2.0 Internal point specification loop

x, y, (z) Coordinates of internal points

Specify as many as required.