sdlen3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	sdlen3 (voln, deltax, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, n1x, n2x, n3x, n4x, n5x, n6x, n1y, n2y, n3y, n4y, n5y, n6y, n1z, n2z, n3z, n4z, n5z, n6z, nel, mtn, jale, jeul)

Function/Subroutine Documentation

sdlen3()

subroutine sdlen3	(		voln,
			deltax,
			x1,
			x2,
			x3,
			x4,
			x5,
			x6,
			x7,
			x8,
			y1,
			y2,
			y3,
			y4,
			y5,
			y6,
			y7,
			y8,
			z1,
			z2,
			z3,
			z4,
			z5,
			z6,
			z7,
			z8,
			n1x,
			n2x,
			n3x,
			n4x,
			n5x,
			n6x,
			n1y,
			n2y,
			n3y,
			n4y,
			n5y,
			n6y,
			n1z,
			n2z,
			n3z,
			n4z,
			n5z,
			n6z,
		integer, intent(in)	nel,
		integer, intent(in)	mtn,
		integer, intent(in)	jale,
		integer, intent(in)	jeul )

Definition at line 37 of file sdlen3.F.

C-----------------------------------------------
C   M o d u l e  s
C-----------------------------------------------
      USE ale_mod
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, INTENT(IN) :: MTN
      INTEGER, INTENT(IN) :: JALE
      INTEGER, INTENT(IN) :: JEUL
      my_real
     .   voln(*),deltax(*),
     .   x1(*), x2(*), x3(*), x4(*), x5(*), x6(*), x7(*), x8(*),
     .   y1(*), y2(*), y3(*), y4(*), y5(*), y6(*), y7(*), y8(*),  
     .   z1(*), z2(*), z3(*), z4(*), z5(*), z6(*), z7(*), z8(*),
     .   n1x(*), n2x(*), n3x(*), n4x(*), n5x(*), n6x(*),
     .   n1y(*), n2y(*), n3y(*), n4y(*), n5y(*), n6y(*),
     .   n1z(*), n2z(*), n3z(*), n4z(*), n5z(*), n6z(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr17_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, NFAC
      my_real xioff(mvsiz), aream(mvsiz), atest(mvsiz), area(mvsiz,6)
C-----------------------------------------------
      DO i=1,nel
        xioff(i) = one
        aream(i) = em20
      ENDDO
 
      IF((ale%GLOBAL%ICAA==1 .OR. ale%GLOBAL%ISFINT>=2) .AND. jale+jeul/=0)THEN
        CALL slena(
     1   n1x,     n2x,     n3x,     n4x,
     2   n5x,     n6x,     n1y,     n2y,
     3   n3y,     n4y,     n5y,     n6y,
     4   n1z,     n2z,     n3z,     n4z,
     5   n5z,     n6z,     x1,      x2,
     6   x3,      x4,      x5,      x6,
     7   x7,      x8,      y1,      y2,
     8   y3,      y4,      y5,      y6,
     9   y7,      y8,      z1,      z2,
     a   z3,      z4,      z5,      z6,
     b   z7,      z8,      aream,   nel)
        DO i=1,nel
          deltax(i)=two*voln(i)*xioff(i)/sqrt(aream(i))
        ENDDO
      ELSE
        CALL slen(
     1   x1,      x2,      x3,      x4,
     2   y1,      y2,      y3,      y4,
     3   z1,      z2,      z3,      z4,
     4   1,       area,    aream,   nel)
        CALL slen(
     1   x5,      x6,      x7,      x8,
     2   y5,      y6,      y7,      y8,
     3   z5,      z6,      z7,      z8,
     4   2,       area,    aream,   nel)
        CALL slen(
     1   x1,      x2,      x6,      x5,
     2   y1,      y2,      y6,      y5,
     3   z1,      z2,      z6,      z5,
     4   3,       area,    aream,   nel)
        CALL slen(
     1   x2,      x3,      x7,      x6,
     2   y2,      y3,      y7,      y6,
     3   z2,      z3,      z7,      z6,
     4   4,       area,    aream,   nel)
        CALL slen(
     1   x3,      x4,      x8,      x7,
     2   y3,      y4,      y8,      y7,
     3   z3,      z4,      z8,      z7,
     4   5,       area,    aream,   nel)
        CALL slen(
     1   x4,      x1,      x5,      x8,
     2   y4,      y1,      y5,      y8,
     3   z4,      z1,      z5,      z8,
     4   6,       area,    aream,   nel)
        IF (idts6>0) THEN
          CALL sldege(
     1   x1,      x2,      x3,      x4,
     2   x5,      x6,      x7,      x8,
     3   y1,      y2,      y3,      y4,
     4   y5,      y6,      y7,      y8,
     5   z1,      z2,      z3,      z4,
     6   z5,      z6,      z7,      z8,
     7   area,    aream,   voln,    nel)
        ENDIF
        IF(mtn==5.OR.mtn==41)THEN
          DO i=1,nel
            atest(i)=em4*aream(i)
          ENDDO       
          DO i=1,nel                            
            nfac=0                                
            IF(area(i,1)<atest(i)) nfac=nfac+1 
            IF(area(i,2)<atest(i)) nfac=nfac+1 
            IF(area(i,3)<atest(i)) nfac=nfac+1 
            IF(area(i,4)<atest(i)) nfac=nfac+1 
            IF(area(i,5)<atest(i)) nfac=nfac+1 
            IF(area(i,6)<atest(i)) nfac=nfac+1 
            IF(nfac>=3) xioff(i)=ep03           
          ENDDO                                          
       END IF
       DO  i=1,nel
         deltax(i)=four*voln(i)*xioff(i)/sqrt(aream(i))
       ENDDO      
      ENDIF
C-----------------------------------------------      
      RETURN