i7pen3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	i7pen3 (marge, gapv, n1, n2, n3, pene, nx1, ny1, nz1, nx2, ny2, nz2, nx3, ny3, nz3, nx4, ny4, nz4, p1, p2, p3, p4, last)

Function/Subroutine Documentation

i7pen3()

subroutine i7pen3	(		marge,
			gapv,
		intent(inout)	n1,
		intent(inout)	n2,
		intent(inout)	n3,
		intent(inout)	pene,
		intent(inout)	nx1,
		intent(inout)	ny1,
		intent(inout)	nz1,
		intent(inout)	nx2,
		intent(inout)	ny2,
		intent(inout)	nz2,
		intent(inout)	nx3,
		intent(inout)	ny3,
		intent(inout)	nz3,
		intent(inout)	nx4,
		intent(inout)	ny4,
		intent(inout)	nz4,
		intent(inout)	p1,
		intent(inout)	p2,
		intent(inout)	p3,
		intent(inout)	p4,
		integer, intent(in)	last )

Parameters

[in]

last

last candidat

Definition at line 38 of file i7pen3.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) :: LAST !< last candidat
C     REAL
      my_real
     .        marge, gapv(*)
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: n1,n2,n3,pene
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: nx1,ny1,nz1
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: nx2,ny2,nz2
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: nx3,ny3,nz3
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: nx4,ny4,nz4
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: p1,p2,p3,p4
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .     d1,d2,d3,d4,s2
C-----------------------------------------------
       DO i=1,last
        d1 = sqrt(p1(i))
        p1(i) = max(zero, gapv(i) + marge - d1)
C
        d2 = sqrt(p2(i))
        p2(i) = max(zero, gapv(i) + marge - d2)
C
        d3 = sqrt(p3(i))
        p3(i) = max(zero, gapv(i) + marge - d3)
C
        d4 = sqrt(p4(i))
        p4(i) = max(zero, gapv(i) + marge - d4)
       ENDDO
C
       DO i=1,last
        pene(i) = max(p1(i),p2(i),p3(i),p4(i))
        if(p1(i)==pene(i))then
         n1(i) = nx1(i)
         n2(i) = ny1(i) 
         n3(i) = nz1(i) 
        elseif(p2(i)==pene(i))then
         n1(i) = nx2(i)
         n2(i) = ny2(i) 
         n3(i) = nz2(i) 
        elseif(p3(i)==pene(i))then
         n1(i) = nx3(i)
         n2(i) = ny3(i) 
         n3(i) = nz3(i) 
        elseif(p4(i)==pene(i))then
         n1(i) = nx4(i)
         n2(i) = ny4(i) 
         n3(i) = nz4(i) 
        endif
       ENDDO
       DO i=1,last
         s2 = one/max(em30,sqrt(n1(i)**2 + n2(i)**2 + n3(i)**2))
         n1(i) = n1(i)*s2
         n2(i) = n2(i)*s2
         n3(i) = n3(i)*s2
       ENDDO
C
      RETURN