ssort_n3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	ssort_n3 (xi, yi, zi, x1, x2, x3, y1, y2, y3, z1, z2, z3, marge, stif, ifc1, nel)

Function/Subroutine Documentation

ssort_n3()

subroutine ssort_n3	(	intent(in)	xi,
		intent(in)	yi,
		intent(in)	zi,
		double precision, dimension(mvsiz), intent(in)	x1,
		double precision, dimension(mvsiz), intent(in)	x2,
		double precision, dimension(mvsiz), intent(in)	x3,
		double precision, dimension(mvsiz), intent(in)	y1,
		double precision, dimension(mvsiz), intent(in)	y2,
		double precision, dimension(mvsiz), intent(in)	y3,
		double precision, dimension(mvsiz), intent(in)	z1,
		double precision, dimension(mvsiz), intent(in)	z2,
		double precision, dimension(mvsiz), intent(in)	z3,
		intent(in)	marge,
		intent(in)	stif,
		integer, dimension(mvsiz), intent(inout)	ifc1,
		integer, intent(in)	nel )

Definition at line 28 of file ssort_n3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT (IN)  :: NEL
      INTEGER, DIMENSION(MVSIZ),INTENT (INOUT) :: IFC1
      my_real, DIMENSION(MVSIZ), INTENT (IN) :: marge,
     .                        xi,     yi,     zi,stif 
      DOUBLE PRECISION, DIMENSION(MVSIZ), INTENT (IN) ::
     .                        X1,     X2,      X3,
     .                        Y1,     Y2,      Y3,
     .                        Z1,     Z2,      Z3
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,K
C                                                                     12
      my_real
     .   rx, ry, rz, sx, sy, sz,nx,ny,nz,bbb,
     .   ux,uy,uz,tol,pene(mvsiz),norm
C----------------------------
         DO i=1,nel
           IF (ifc1(i)>0) cycle
           rx =x2(i)-x1(i)
           ry =y2(i)-y1(i)
           rz =z2(i)-z1(i)
           sx =x3(i)-x1(i)
           sy =y3(i)-y1(i)
           sz =z3(i)-z1(i)
           nx =ry*sz - rz*sy
           ny =rz*sx - rx*sz
           nz =rx*sy - ry*sx
           norm=one/max(em20,sqrt(nx*nx+ny*ny+nz*nz))
           bbb = ((x3(i)-xi(i))*nx + (y3(i)-yi(i))*ny +
     .            (z3(i)-zi(i))*nz)*norm 
           pene(i) = abs(bbb)  
           IF (pene(i)<marge(i).AND.stif(i)>zero) ifc1(i)=2
         ENDDO
C       
      RETURN