in_out_side.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	in_out_side (inod, inod2surf, knod2surf, nnod2surf, x, xn, dist, nseg, surf_eltyp, nod_normal, surf_nodes, swiftsurf, idsurf, segtosurf)

Function/Subroutine Documentation

in_out_side()

subroutine in_out_side	(	integer	inod,
		integer, dimension(nnod2surf,*)	inod2surf,
		integer, dimension(*)	knod2surf,
		integer	nnod2surf,
			x,
			xn,
			dist,
		integer	nseg,
		integer, dimension(nseg)	surf_eltyp,
			nod_normal,
		integer, dimension(nseg,4)	surf_nodes,
		integer, dimension(nsurf)	swiftsurf,
		integer	idsurf,
		integer, dimension(*)	segtosurf )

Definition at line 31 of file in_out_side.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER INOD,KNOD2SURF(*),NNOD2SURF,
     .  INOD2SURF(NNOD2SURF,*),NSEG,SURF_ELTYP(NSEG),
     .  SURF_NODES(NSEG,4),SWIFTSURF(NSURF),IDSURF,SEGTOSURF(*)
      my_real
     . x(3,*),xn(3),dist,
     . nod_normal(3,*),xfas(3,4)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,K,OK,FAC,FAC1,FAC2,NULL_DIST,
     . IPL,IXPL(4),
     . p,r,p1,p2,dd1(4),dd2(4),ddd1(3),ddd2(3),ITYP,SH_ELEM,ISURF
      my_real
     . dist_pl(3),vx_nod_inod,vy_nod_inod,vz_nod_inod,xsign(3),
     . v1x,v1y,v1z,v2x,v2y,v2z,v3x,v3y,v3z,v12x,v12y,v12z,dot,
     . nx,ny,nz
      DATA dd1/4,1,2,3/,dd2/2,3,4,1/,ddd1/3,1,2/,ddd2/2,3,1/
C-----------------------------------------------
        ok = 0
        null_dist = 0
        DO 60 k=1,knod2surf(inod)
          ipl = inod2surf(k,inod)
          isurf = segtosurf(ipl)
          IF (isurf /= idsurf) GOTO 60
          ipl = ipl - swiftsurf(isurf)
          IF (ipl <= 0 .OR. ipl > nseg .OR. inod > numnod) GOTO 60
          ityp = surf_eltyp(ipl)
!
          IF(ityp==3)THEN
            ixpl(1) = surf_nodes(ipl,1)
            ixpl(2) = surf_nodes(ipl,2)
            ixpl(3) = surf_nodes(ipl,3)
            ixpl(4) = surf_nodes(ipl,4)
          ELSEIF(ityp==7)THEN
            ixpl(1) = surf_nodes(ipl,1)
            ixpl(2) = surf_nodes(ipl,2)
            ixpl(3) = surf_nodes(ipl,3)
            ixpl(4) = surf_nodes(ipl,3)
          ENDIF
C---
          p = 0
C---
          IF(ityp==3)THEN
           DO r=1,4
             IF(inod == ixpl(r))THEN
               p = r
               EXIT
             ENDIF
           ENDDO
          ELSEIF(ityp==7)THEN
           DO r=1,3
             IF(inod == ixpl(r))THEN
               p = r
               EXIT
             ENDIF
           ENDDO
          ENDIF
C
          IF(p == 0)GOTO 60
C
          IF(ityp==3)THEN
            p1 = dd1(p)
            p2 = dd2(p)
          ELSEIF(ityp==7)THEN
            p1 = ddd1(p)
            p2 = ddd2(p)
          ENDIF
C
          xfas(1,1) = x(1,ixpl(1))
          xfas(2,1) = x(2,ixpl(1))
          xfas(3,1) = x(3,ixpl(1))
          xfas(1,2) = x(1,ixpl(2))
          xfas(2,2) = x(2,ixpl(2))
          xfas(3,2) = x(3,ixpl(2))
          xfas(1,3) = x(1,ixpl(3))
          xfas(2,3) = x(2,ixpl(3))
          xfas(3,3) = x(3,ixpl(3))
          xfas(1,4) = x(1,ixpl(4))
          xfas(2,4) = x(2,ixpl(4))
          xfas(3,4) = x(3,ixpl(4))
C
          CALL nfacette(xfas,nx,ny,nz)
C
          dist_pl(1:3) = zero
C vect NOD --> INOD
          vx_nod_inod = xn(1) - x(1,inod)
          vy_nod_inod = xn(2) - x(2,inod)
          vz_nod_inod = xn(3) - x(3,inod)
C PLANE 1
C V1xV2
          v12x = nx
          v12y = ny
          v12z = nz
C dist of N to plane 1
          dot = vx_nod_inod*v12x+vy_nod_inod*v12y+vz_nod_inod*v12z
          dist_pl(1) = dot
C PLANE 2
C V_NOR x V1
          v1x = nod_normal(1,inod)
          v1y = nod_normal(2,inod)
          v1z = nod_normal(3,inod)
          v2x = xfas(1,p)-xfas(1,p1)
          v2y = xfas(2,p)-xfas(2,p1)
          v2z = xfas(3,p)-xfas(3,p1)
          v12x = v1y*v2z-v1z*v2y
          v12y = v1z*v2x-v1x*v2z
          v12z = v1x*v2y-v1y*v2x
C dist of N to plane 2
          dot = vx_nod_inod*v12x+vy_nod_inod*v12y+vz_nod_inod*v12z
          dist_pl(2) = dot
C PLANE 3
C V_NOR x V2 
          v1x = nod_normal(1,inod)
          v1y = nod_normal(2,inod)
          v1z = nod_normal(3,inod)
          v2x = xfas(1,p2)-xfas(1,p)
          v2y = xfas(2,p2)-xfas(2,p)
          v2z = xfas(3,p2)-xfas(3,p)
          v12x = v1y*v2z-v1z*v2y
          v12y = v1z*v2x-v1x*v2z
          v12z = v1x*v2y-v1y*v2x
C dist of N to plane 3
          dot = vx_nod_inod*v12x+vy_nod_inod*v12y+vz_nod_inod*v12z
          dist_pl(3) = dot
C---
          DO r=1,3
            xsign(r) = sign(one,dist_pl(r))
            IF(dist_pl(r) == zero)xsign(r) = one !ZERO
          ENDDO
          fac1 = 0
          fac2 = 0
C
          IF(null_dist == 1)EXIT
C
          DO r=1,3
            IF(xsign(r) >= zero)THEN
              fac1 = fac1 + 1
            ENDIF
          ENDDO
          IF(fac1 == 3)THEN
            dist = one
            ok = ok + 1
          ENDIF
          IF(ok > 0)EXIT
C---
 60     CONTINUE
C---
c        IF(NULL_DIST == 1) DIST = ZERO
        IF(ok == 0) dist = -one  ! external node
        IF(null_dist == 1) dist = zero
C-----------------------------------------------
      RETURN