r5def3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"
#include "com08_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	r5def3 (xl, vx2l, ry1l, rz1l, rx2l, ry2l, rz2l, fr_wave, fr_w_e, eint, fx, xmom, ymom, zmom, fy, fz, partsav, ipartr, exx, eyx, ezx, exy, eyy, ezy, exz, eyz, ezz, rx1, ry1, rz1, rx2, ry2, rz2, vx1, vx2, vy1, vy2, vz1, vz2, nc1, nc2, nel)

Function/Subroutine Documentation

r5def3()

subroutine r5def3	(		xl,
			vx2l,
			ry1l,
			rz1l,
			rx2l,
			ry2l,
			rz2l,
			fr_wave,
			fr_w_e,
			eint,
			fx,
			xmom,
			ymom,
			zmom,
			fy,
			fz,
			partsav,
		integer, dimension(*)	ipartr,
			exx,
			eyx,
			ezx,
			exy,
			eyy,
			ezy,
			exz,
			eyz,
			ezz,
			rx1,
			ry1,
			rz1,
			rx2,
			ry2,
			rz2,
			vx1,
			vx2,
			vy1,
			vy2,
			vz1,
			vz2,
		integer, dimension(*)	nc1,
		integer, dimension(*)	nc2,
		integer, intent(in)	nel )

Definition at line 28 of file r5def3.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-----------------------------------------------
#include      "com01_c.inc"
#include      "com08_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER IPARTR(*),NC1(*),NC2(*)
C     REAL
      my_real
     .   xl(*), vx2l(*),ry1l(*), rz1l(*), rx2l(*), ry2l(*), rz2l(*),
     .   fr_wave(*) ,fr_w_e(*),eint(*) ,  
     .   fx(*), fy(*), fz(*), xmom(*), ymom(*),zmom(*),partsav(npsav,*),
     .   exx(mvsiz), eyx(mvsiz), ezx(mvsiz),
     .   exy(mvsiz), eyy(mvsiz), ezy(mvsiz),
     .   exz(mvsiz), eyz(mvsiz), ezz(mvsiz), rx1(mvsiz), rx2(mvsiz),
     .   ry1(mvsiz), ry2(mvsiz), rz1(mvsiz), rz2(mvsiz), vx1(mvsiz),
     .   vx2(mvsiz), vy1(mvsiz), vy2(mvsiz), vz1(mvsiz), vz2(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,MX
C     REAL
      my_real
     .     vx1l, vy1l, vy2l, vz1l, vz2l,theta , xldemi,xsign
C-----------------------------------------------
C--------------------------------------------
C     VITESSES REPERE CONVECTEE
C--------------------------------------------
      DO i=1,nel
        rx2l(i) = exx(i)*(rx2(i)-rx1(i))
     .          + eyx(i)*(ry2(i)-ry1(i))
     .          + ezx(i)*(rz2(i)-rz1(i))
        ry1l(i) = exy(i)*rx1(i)+eyy(i)*ry1(i)+ezy(i)*rz1(i)
        ry2l(i) = exy(i)*rx2(i)+eyy(i)*ry2(i)+ezy(i)*rz2(i)
        rz1l(i) = exz(i)*rx1(i)+eyz(i)*ry1(i)+ezz(i)*rz1(i)
        rz2l(i) = exz(i)*rx2(i)+eyz(i)*ry2(i)+ezz(i)*rz2(i)
        vx2l(i) = exx(i)*(vx2(i)-vx1(i))
     .          + eyx(i)*(vy2(i)-vy1(i))
     .          + ezx(i)*(vz2(i)-vz1(i))
        vy2l   = exy(i)*(vx2(i)-vx1(i))
     .         + eyy(i)*(vy2(i)-vy1(i))
     .         + ezy(i)*(vz2(i)-vz1(i))
        xsign  = sign(one, xl(i) - half*vx2l(i)*dt1)
        xldemi = xsign/max(em15,abs(xl(i) - half*vx2l(i)*dt1))
        theta  = vy2l * xldemi
        rz1l(i) = rz1l(i) - theta
        rz2l(i) = rz2l(i) - theta
        vz2l   = exz(i)*(vx2(i)-vx1(i))
     .         + eyz(i)*(vy2(i)-vy1(i))
     .         + ezz(i)*(vz2(i)-vz1(i))
        theta  = vz2l * xldemi
        ry1l(i) = ry1l(i) + theta
        ry2l(i) = ry2l(i) + theta
        vx2l(i) = vx2l(i)
     .    - half*dt1*xldemi*(vy2l*vy2l+vz2l*vz2l)
      ENDDO
C--------------------------------------------
C     Energy
C--------------------------------------------
      DO i=1,nel
        eint(i) = eint(i)
     .+ half*dt1 * (vx2l(i) * fx(i) + rx2l(i) * xmom(i) 
     .          + (ry2l(i) - ry1l(i)) * ymom(i) 
     .          + (rz2l(i) - rz1l(i)) * zmom(i)
     .          + half * (ry2l(i) + ry1l(i)) * fz(i) * xl(i)
     .          - half * (rz2l(i) + rz1l(i)) * fy(i) * xl(i) )
      ENDDO
C
      IF (npsav >= 21) THEN
        DO i=1,nel
          mx = ipartr(i)
          partsav(23,mx)=partsav(23,mx)
     .          + half*dt1 * (rx2l(i) * xmom(i)
     .          + (ry2l(i) - ry1l(i)) * ymom(i) 
     .          + (rz2l(i) - rz1l(i)) * zmom(i)
     .          + half * (ry2l(i) + ry1l(i)) * fz(i) * xl(i)
     .          - half * (rz2l(i) + rz1l(i)) * fy(i) * xl(i) )
        ENDDO
      ENDIF
C--------------------------------------------
C     Front wave
C--------------------------------------------
      IF (ifrwv /= 0) THEN
        DO i=1,nel
          fr_w_e(i)=max(fr_wave(nc1(i)),fr_wave(nc2(i)),zero)
        ENDDO
      ENDIF
C---
      RETURN