edefo2.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com04_c.inc"
#include "com01_c.inc"
#include "com08_c.inc"
#include "vect01_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	edefo2 (vol, v, veul, y1, y2, y3, y4, z1, z2, z3, z4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, py1, py2, pz1, pz2, wyz, dyz, dzy, eyy, ezz, ett, eyz, eyt, ezt, voln, aire, deltax, vdy, vdz, vd2, nc1, nc2, nc3, nc4)

Function/Subroutine Documentation

edefo2()

subroutine edefo2	(		vol,
			v,
			veul,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4,
			vy1,
			vy2,
			vy3,
			vy4,
			vz1,
			vz2,
			vz3,
			vz4,
			py1,
			py2,
			pz1,
			pz2,
			wyz,
			dyz,
			dzy,
			eyy,
			ezz,
			ett,
			eyz,
			eyt,
			ezt,
			voln,
			aire,
			deltax,
			vdy,
			vdz,
			vd2,
		integer, dimension(*)	nc1,
		integer, dimension(*)	nc2,
		integer, dimension(*)	nc3,
		integer, dimension(*)	nc4 )

Definition at line 29 of file edefo2.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"
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real
     .   vol(*), v(3,numnod), veul(10,*),
     .   y1(*),y2(*),y3(*),y4(*),z1(*),z2(*),z3(*),z4(*),
     .   vy1(*),vy2(*),vy3(*),vy4(*),vz1(*),vz2(*),vz3(*),vz4(*),
     .   py1(*), py2(*), pz1(*), pz2(*),
     .   wyz(*), dyz(*), dzy(*),
     .   eyy(*),ezz(*), ett(*), eyz(*), eyt(*), ezt(*),
     .   voln(*),aire(*),deltax(*),vdy(*),vdz(*),vd2(*)
      INTEGER NC1(*), NC2(*), NC3(*), NC4(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com08_c.inc"
#include      "vect01_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, II
      my_real vy13(mvsiz), vy24(mvsiz), vz13(mvsiz),vz24(mvsiz),yavg(mvsiz)
C-----------------------------------------------
      DO i=lft,llt
        ii=i+nft
        py1(i)=half*(z2(i)-z4(i))
        py2(i)=half*(z3(i)-z1(i))
        pz1(i)=half*(y4(i)-y2(i))
        pz2(i)=half*(y1(i)-y3(i))      
        aire(i)=veul(1,ii)
        deltax(i)=veul(2,ii)
      ENDDO !next I
 
      IF(n2d == 1) THEN
       DO i=lft,llt
         voln(i)=vol(i)
       ENDDO
      ELSE
       DO i=lft,llt
         voln(i)=aire(i)
       ENDDO
      ENDIF
 
      DO i=lft,llt
        vy1(i)=v(2,nc1(i))
        vz1(i)=v(3,nc1(i))
        vy2(i)=v(2,nc2(i))
        vz2(i)=v(3,nc2(i))
        vy3(i)=v(2,nc3(i))
        vz3(i)=v(3,nc3(i))
        vy4(i)=v(2,nc4(i))
        vz4(i)=v(3,nc4(i))
      ENDDO  
 
      IF(mtn /= 11)THEN
       DO i=lft,llt
         yavg(i) =y1(i)+y2(i)+y3(i)+y4(i)
         vy13(i)=vy1(i)-vy3(i)
         vy24(i)=vy2(i)-vy4(i)
         vz13(i)=vz1(i)-vz3(i)
         vz24(i)=vz2(i)-vz4(i)
         vdy(i)=fourth*(vy1(i)+vy2(i)+vy3(i)+vy4(i))
         vdz(i)=fourth*(vz1(i)+vz2(i)+vz3(i)+vz4(i))
         vd2(i)=vdy(i)**2+vdz(i)**2
       ENDDO 
 
       DO i=lft,llt
         eyy(i)=(py1(i)*vy13(i)+py2(i)*vy24(i))/aire(i)
         ezz(i)=(pz1(i)*vz13(i)+pz2(i)*vz24(i))/aire(i)
         ett(i)=zero
         dzy(i)=(py1(i)*vz13(i)+py2(i)*vz24(i))/aire(i)
         dyz(i)=(pz1(i)*vy13(i)+pz2(i)*vy24(i))/aire(i)
         eyt(i)=zero
         ezt(i)=zero
       ENDDO 
 
       IF(n2d == 1) THEN
        DO i=lft,llt
         ett(i)=(vy1(i)+vy2(i)+vy3(i)+vy4(i))/yavg(i)
        ENDDO
       ENDIF
 
       DO i=lft,llt
         eyz(i)= dzy(i)+dyz(i)
         wyz(i)=half*dt1*(dzy(i)-dyz(i))
       ENDDO
      ENDIF
C-----------------------------------------------
      RETURN