ede112.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"
#include "com04_c.inc"
#include "com08_c.inc"
#include "vect01_c.inc"
#include "tabsiz_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	ede112 (pm, v, x, ixq, ale_connect, wyz, dyz, dzy, eyy, ezz)

Function/Subroutine Documentation

ede112()

subroutine ede112	(	dimension(npropm,nummat), intent(in)	pm,
		dimension(3,numnod), intent(in)	v,
		dimension(3,sx/3), intent(in)	x,
		integer, dimension(nixq,sixq/nixq), intent(in)	ixq,
		type(t_ale_connectivity), intent(in)	ale_connect,
		dimension(*), intent(inout)	wyz,
		dimension(*), intent(inout)	dyz,
		dimension(*), intent(inout)	dzy,
		dimension(*), intent(inout)	eyy,
		dimension(*), intent(inout)	ezz )

Definition at line 30 of file ede112.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE ale_connectivity_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
#include      "com04_c.inc"
#include      "com08_c.inc"
#include      "vect01_c.inc"
#include      "tabsiz_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
! SPMD CASE : SIXQ >= NIXQ*NUMELQ    (SIXQ = NIXQ*NUMELQ_L+NIXQ*NQVOIS_L)
! IXQ(1:NIXQ, 1:NUMELQ) local elems
!    (1:NIXQ, NUMELQ+1:) additional elems (also on adjacent domains but connected to the boundary of the current domain)
      INTEGER,INTENT(IN) :: IXQ(NIXQ,SIXQ/NIXQ)
! SPMD CASE : SX >= 3*NUMNOD    (SX = 3*(NUMNOD_L+NRCVVOIS_L))
! X(1:3,1:NUMNOD) : local nodes
!  (1:3, NUMNOD+1:) additional nodes (also on adjacent domains but connected to the boundary of the current domain)
      my_real,INTENT(IN) :: x(3,sx/3)
!
      my_real,INTENT(IN) :: v(3,numnod), pm(npropm,nummat)
      my_real,INTENT(INOUT) :: eyy(*),ezz(*),dyz(*), dzy(*),wyz(*)
      TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, II,IADJ,ML,NC1,NC2,NC3,NC4,J,IAD2
      my_real vy13, vy24, vz13,vz24,py1,py2,pz1,pz2,aire
C-----------------------------------------------
C   S o u r c e   L i n e s 
C-----------------------------------------------
      DO i=lft,llt
       ii=i+nft
       iad2 = ale_connect%ee_connect%iad_connect(ii)
       DO j=1,4
         iadj=ale_connect%ee_connect%connected(iad2 + j - 1)
         IF(iadj /= 0) THEN
           ml=nint(pm(19,ixq(1,iadj)))
           IF(ml /= 11) THEN
             nc1=ixq(2,iadj)
             nc2=ixq(3,iadj)
             nc3=ixq(4,iadj)
             nc4=ixq(5,iadj)
             py1=half*(x(3,nc2)-x(3,nc4))
             py2=half*(x(3,nc3)-x(3,nc1))
             pz1=half*(x(2,nc4)-x(2,nc2))
             pz2=half*(x(2,nc1)-x(2,nc3))    
             aire=two*(py1*pz2-pz1*py2)
             vy13=v(2,nc1)-v(2,nc3)
             vy24=v(2,nc2)-v(2,nc4)
             vz13=v(3,nc1)-v(3,nc3)
             vz24=v(3,nc2)-v(3,nc4)
             eyy(i)=(py1*vy13+py2*vy24)/aire
             ezz(i)=(pz1*vz13+pz2*vz24)/aire
             dzy(i)=(py1*vz13+py2*vz24)/aire
             dyz(i)=(pz1*vy13+pz2*vy24)/aire
             ezz(i)=zero
             wyz(i)=half*dt1*(dzy(i)-dyz(i))
           ENDIF
         ENDIF
       enddo!next J
      enddo!next I
C-----------------------------------------------
      RETURN