sdlenmax.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	sdlenmax (lm_max, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, nel)

Function/Subroutine Documentation

sdlenmax()

subroutine sdlenmax	(	intent(out)	lm_max,
		intent(in)	x1,
		intent(in)	x2,
		intent(in)	x3,
		intent(in)	x4,
		intent(in)	x5,
		intent(in)	x6,
		intent(in)	x7,
		intent(in)	x8,
		intent(in)	y1,
		intent(in)	y2,
		intent(in)	y3,
		intent(in)	y4,
		intent(in)	y5,
		intent(in)	y6,
		intent(in)	y7,
		intent(in)	y8,
		intent(in)	z1,
		intent(in)	z2,
		intent(in)	z3,
		intent(in)	z4,
		intent(in)	z5,
		intent(in)	z6,
		intent(in)	z7,
		intent(in)	z8,
		integer, intent(in)	nel )

Definition at line 32 of file sdlenmax.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   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      my_real, DIMENSION(MVSIZ) , INTENT(OUT)  :: lm_max 
      my_real, DIMENSION(MVSIZ) , INTENT(IN)  ::  
     .   x1, x2, x3, x4, x5, x6, x7, x8,
     .   y1, y2, y3, y4, y5, y6, y7, y8,  
     .   z1, z2, z3, z4, z5, z6, z7, z8
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, N
      my_real
     .   rx, ry, rz, sx, sy, sz, tx, ty, tz,norm
C=======================================================================
      DO i=1,nel
        rx = x1(i)+x2(i)+x5(i)+x6(i)-x3(i)-x4(i)-x7(i)-x8(i)
        ry = y1(i)+y2(i)+y5(i)+y6(i)-y3(i)-y4(i)-y7(i)-y8(i)
        rz = z1(i)+z2(i)+z5(i)+z6(i)-z3(i)-z4(i)-z7(i)-z8(i)
        sx = x5(i)+x6(i)+x7(i)+x8(i)-x1(i)-x2(i)-x3(i)-x4(i)
        sy = y5(i)+y6(i)+y7(i)+y8(i)-y1(i)-y2(i)-y3(i)-y4(i)
        sz = z5(i)+z6(i)+z7(i)+z8(i)-z1(i)-z2(i)-z3(i)-z4(i)
        tx = x3(i)+x2(i)+x7(i)+x6(i)-x1(i)-x4(i)-x5(i)-x8(i)
        ty = y3(i)+y2(i)+y7(i)+y6(i)-y1(i)-y4(i)-y5(i)-y8(i)
        tz = z3(i)+z2(i)+z7(i)+z6(i)-z1(i)-z4(i)-z5(i)-z8(i)
        norm = rx*rx + ry*ry + rz*rz
        lm_max(i) = norm
        norm = sx*sx + sy*sy + sz*sz
        lm_max(i) = max(lm_max(i),norm)
        norm = tx*tx + ty*ty + tz*tz
        lm_max(i) = max(lm_max(i),norm)
      ENDDO
      DO i=1,nel
        lm_max(i) = fourth*sqrt(lm_max(i))
      ENDDO 
C
      RETURN