a4momtn3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	a4momtn3 (rho, vol, vis, r, s, t, deltax, vmx, vmy, vmz, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, px, py, pz, vdx, vdy, vdz, fx, fy, fz, nel)

Function/Subroutine Documentation

a4momtn3()

subroutine a4momtn3	(		rho,
			vol,
			vis,
			r,
			s,
			t,
			deltax,
			vmx,
			vmy,
			vmz,
			dxx,
			dxy,
			dxz,
			dyx,
			dyy,
			dyz,
			dzx,
			dzy,
			dzz,
			px,
			py,
			pz,
			vdx,
			vdy,
			vdz,
			fx,
			fy,
			fz,
		integer, intent(in)	nel )

Definition at line 30 of file a4momtn3.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-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      my_real
     .   rho(*),vol(*),vis(*),vdx(*),vdy(*),vdz(*),vmx(*),vmy(*),vmz(*),
     .   r(mvsiz,3),s(mvsiz,3),t(mvsiz,3),deltax(*),px(*),py(*),pz(*),
     .   fx(*),fy(*),fz(*),
     .   dxx(*),dxy(*),dxz(*),
     .   dyx(*),dyy(*),dyz(*),
     .   dzx(*),dzy(*),dzz(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I
      my_real :: a(mvsiz),gam(mvsiz) 
C-----------------------------------------------
C   S o u r c e   C o d e
C-----------------------------------------------
      CALL upwind4(
     1   rho,     vis,     vdx,     vdy,
     2   vdz,     r,       s,       t,
     3   gam,     nel)
      DO i=1,nel
        a(i)=fourth*rho(i)*vol(i)*(one+four*gam(i)*( px(i)*vdx(i)+py(i)*vdy(i)+pz(i)*vdz(i)))
        fx(i)=fx(i)-a(i)*(vdx(i)*dxx(i)+vdy(i)*dxy(i)+vdz(i)*dxz(i))
        fy(i)=fy(i)-a(i)*(vdx(i)*dyx(i)+vdy(i)*dyy(i)+vdz(i)*dyz(i))
        fz(i)=fz(i)-a(i)*(vdx(i)*dzx(i)+vdy(i)*dzy(i)+vdz(i)*dzz(i))
      ENDDO
      RETURN