sdef43.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	sdef43 (nel, npg, hh, dzx, dyz, dzz, vxloc, vyloc, vzloc)

Function/Subroutine Documentation

sdef43()

subroutine sdef43	(	integer	nel,
		integer	npg,
			hh,
			dzx,
			dyz,
			dzz,
			vxloc,
			vyloc,
			vzloc )

Definition at line 28 of file sdef43.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
CC-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include "mvsiz_p.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL,NPG
C     REAL
      my_real 
     .   vxloc(mvsiz,8),vyloc(mvsiz,8),vzloc(mvsiz,8),hh(4,npg),
     .   dzx(mvsiz,npg),dyz(mvsiz,npg),dzz(mvsiz,npg)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IEL,IPG
C     REAL
      my_real
     .   pg,p1,p2,p3,dxinf,dyinf,dzinf,dxsup,dysup,dzsup
      DATA   pg/0.5773502691896257645091488/
C=======================================================================
      p1 = (one + pg)*(one - pg) * fourth
      p2 = (one + pg)*(one + pg) * fourth
      p3 = (one - pg)*(one - pg) * fourth
C
      hh(1,1) = p2
      hh(2,1) = p1
      hh(3,1) = p3
      hh(4,1) = p1
C
      hh(1,2) = p1
      hh(2,2) = p2
      hh(3,2) = p1
      hh(4,2) = p3
C
      hh(1,3) = p3
      hh(2,3) = p1
      hh(3,3) = p2
      hh(4,3) = p1
C
      hh(1,4) = p1
      hh(2,4) = p3
      hh(3,4) = p1
      hh(4,4) = p2
C
      DO ipg = 1,npg
        DO iel = 1,nel
          dxsup   = hh(1,ipg)*vxloc(iel,5)+hh(2,ipg)*vxloc(iel,6)
     .            + hh(3,ipg)*vxloc(iel,7)+hh(4,ipg)*vxloc(iel,8)
          dxinf   = hh(1,ipg)*vxloc(iel,1)+hh(2,ipg)*vxloc(iel,2)
     .            + hh(3,ipg)*vxloc(iel,3)+hh(4,ipg)*vxloc(iel,4)
          dysup   = hh(1,ipg)*vyloc(iel,5)+hh(2,ipg)*vyloc(iel,6)
     .            + hh(3,ipg)*vyloc(iel,7)+hh(4,ipg)*vyloc(iel,8)
          dyinf   = hh(1,ipg)*vyloc(iel,1)+hh(2,ipg)*vyloc(iel,2)
     .            + hh(3,ipg)*vyloc(iel,3)+hh(4,ipg)*vyloc(iel,4)
          dzsup   = hh(1,ipg)*vzloc(iel,5)+hh(2,ipg)*vzloc(iel,6)
     .            + hh(3,ipg)*vzloc(iel,7)+hh(4,ipg)*vzloc(iel,8)
          dzinf   = hh(1,ipg)*vzloc(iel,1)+hh(2,ipg)*vzloc(iel,2)
     .            + hh(3,ipg)*vzloc(iel,3)+hh(4,ipg)*vzloc(iel,4)
C
          dzz(iel,ipg) = dzsup - dzinf
          dyz(iel,ipg) = dysup - dyinf
          dzx(iel,ipg) = dxsup - dxinf
        ENDDO
      ENDDO
C-----------
      RETURN