sdlensh14.F File Reference

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

Go to the source code of this file.

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

Function/Subroutine Documentation

sdlensh14()

subroutine sdlensh14	(	intent(in)	voln,
		intent(out)	llsh,
		intent(out)	area,
		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)	ics,
		integer, intent(in)	nel )

Definition at line 28 of file sdlensh14.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,ICS
      my_real, DIMENSION(MVSIZ) , INTENT(OUT)  :: area,llsh 
      my_real, DIMENSION(MVSIZ) , INTENT(IN)  ::  
     .   voln,
     .   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(mvsiz),ry(mvsiz),rz(mvsiz),sx(mvsiz),sy(mvsiz),sz(mvsiz),
     .   xn(mvsiz,4) , yn(mvsiz,4) , zn(mvsiz,4) 
      my_real
     .   al1,al2,e3x,e3y,e3z,det
C=======================================================================
        SELECT CASE(ics)
          CASE (1)
            DO i=1,nel
              xn(i,1) = half*(x1(i)+x4(i))
              yn(i,1) = half*(y1(i)+y4(i))
              zn(i,1) = half*(z1(i)+z4(i))
              xn(i,2) = half*(x2(i)+x3(i))
              yn(i,2) = half*(y2(i)+y3(i))
              zn(i,2) = half*(z2(i)+z3(i))
              xn(i,3) = half*(x6(i)+x7(i))
              yn(i,3) = half*(y6(i)+y7(i))
              zn(i,3) = half*(z6(i)+z7(i))
              xn(i,4) = half*(x5(i)+x8(i))
              yn(i,4) = half*(y5(i)+y8(i))
              zn(i,4) = half*(z5(i)+z8(i))
            ENDDO 
          CASE (10)
            DO i=1,nel
              xn(i,1) = half*(x1(i)+x5(i))
              yn(i,1) = half*(y1(i)+y5(i))
              zn(i,1) = half*(z1(i)+z5(i))
              xn(i,2) = half*(x2(i)+x6(i))
              yn(i,2) = half*(y2(i)+y6(i))
              zn(i,2) = half*(z2(i)+z6(i))
              xn(i,3) = half*(x3(i)+x7(i))
              yn(i,3) = half*(y3(i)+y7(i))
              zn(i,3) = half*(z3(i)+z7(i))
              xn(i,4) = half*(x4(i)+x8(i))
              yn(i,4) = half*(y4(i)+y8(i))
              zn(i,4) = half*(z4(i)+z8(i))
            ENDDO 
          CASE (100)
            DO i=1,nel
              xn(i,1) = half*(x1(i)+x2(i))
              yn(i,1) = half*(y1(i)+y2(i))
              zn(i,1) = half*(z1(i)+z2(i))
              xn(i,2) = half*(x5(i)+x6(i))
              yn(i,2) = half*(y5(i)+y6(i))
              zn(i,2) = half*(z5(i)+z6(i))
              xn(i,3) = half*(x8(i)+x7(i))
              yn(i,3) = half*(y8(i)+y7(i))
              zn(i,3) = half*(z8(i)+z7(i))
              xn(i,4) = half*(x4(i)+x3(i))
              yn(i,4) = half*(y4(i)+y3(i))
              zn(i,4) = half*(z4(i)+z3(i))
            ENDDO 
        END SELECT
C------g1,g2 :
        DO i=1,nel
          rx(i)=xn(i,2)+xn(i,3)-xn(i,1)-xn(i,4)
          ry(i)=yn(i,2)+yn(i,3)-yn(i,1)-yn(i,4)
          rz(i)=zn(i,2)+zn(i,3)-zn(i,1)-zn(i,4)
          sx(i)=xn(i,3)+xn(i,4)-xn(i,1)-xn(i,2)
          sy(i)=yn(i,3)+yn(i,4)-yn(i,1)-yn(i,2)
          sz(i)=zn(i,3)+zn(i,4)-zn(i,1)-zn(i,2)
        ENDDO 
        DO i=1,nel
          e3x = ry(i) * sz(i) - rz(i) * sy(i) 
          e3y = rz(i) * sx(i) - rx(i) * sz(i) 
          e3z = rx(i) * sy(i) - ry(i) * sx(i) 
          det = sqrt(e3x*e3x + e3y*e3y + e3z*e3z)
          area(i)=fourth*det
          al1 = rx(i)*rx(i) + ry(i)*ry(i)+ rz(i)*rz(i)
          al2 = sx(i)*sx(i) + sy(i)*sy(i)+ sz(i)*sz(i)
          llsh(i) = fourth*max(al1,al2)
        ENDDO 
C
      RETURN