fe_close.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	fe_close (hh, rho, voln, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nel)

Function/Subroutine Documentation

fe_close()

subroutine fe_close	(		hh,
			rho,
			voln,
			vx1,
			vx2,
			vx3,
			vx4,
			vx5,
			vx6,
			vx7,
			vx8,
			vy1,
			vy2,
			vy3,
			vy4,
			vy5,
			vy6,
			vy7,
			vy8,
			vz1,
			vz2,
			vz3,
			vz4,
			vz5,
			vz6,
			vz7,
			vz8,
			f11,
			f21,
			f31,
			f12,
			f22,
			f32,
			f13,
			f23,
			f33,
			f14,
			f24,
			f34,
			f15,
			f25,
			f35,
			f16,
			f26,
			f36,
			f17,
			f27,
			f37,
			f18,
			f28,
			f38,
		integer, intent(in)	nel )

Definition at line 29 of file fe_close.F.

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      "mvsiz_p.inc"
#include      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      my_real
     .   f11(mvsiz),f21(mvsiz),f31(mvsiz),
     .   f12(mvsiz),f22(mvsiz),f32(mvsiz),
     .   f13(mvsiz),f23(mvsiz),f33(mvsiz),
     .   f14(mvsiz),f24(mvsiz),f34(mvsiz),
     .   f15(mvsiz),f25(mvsiz),f35(mvsiz),
     .   f16(mvsiz),f26(mvsiz),f36(mvsiz),
     .   f17(mvsiz),f27(mvsiz),f37(mvsiz),
     .   f18(mvsiz),f28(mvsiz),f38(mvsiz),
     .   hh(*),rho(*),voln(*),
     .   vx1(*),vx2(*),vx3(*),vx4(*),vx5(*),vx6(*),vx7(*),vx8(*),
     .   vy1(*),vy2(*),vy3(*),vy4(*),vy5(*),vy6(*),vy7(*),vy8(*),
     .   vz1(*),vz2(*),vz3(*),vz4(*),vz5(*),vz6(*),vz7(*),vz8(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     . fac  
C      
      IF(dt1==zero)RETURN
      DO i=1,nel
        fac=-one_over_8*hh(i)*rho(i)*voln(i)/dt1
          f11(i) = f11(i) + fac*vx1(i)
          f12(i) = f12(i) + fac*vx2(i)
          f13(i) = f13(i) + fac*vx3(i)
          f14(i) = f14(i) + fac*vx4(i)
          f15(i) = f15(i) + fac*vx5(i)
          f16(i) = f16(i) + fac*vx6(i)
          f17(i) = f17(i) + fac*vx7(i)
          f18(i) = f18(i) + fac*vx8(i)
 
          f21(i) = f21(i) + fac*vy1(i)
          f22(i) = f22(i) + fac*vy2(i)
          f23(i) = f23(i) + fac*vy3(i)
          f24(i) = f24(i) + fac*vy4(i)
          f25(i) = f25(i) + fac*vy5(i)
          f26(i) = f26(i) + fac*vy6(i)
          f27(i) = f27(i) + fac*vy7(i)
          f28(i) = f28(i) + fac*vy8(i)
 
          f31(i) = f31(i) + fac*vz1(i)
          f32(i) = f32(i) + fac*vz2(i)
          f33(i) = f33(i) + fac*vz3(i)
          f34(i) = f34(i) + fac*vz4(i)
          f35(i) = f35(i) + fac*vz5(i)
          f36(i) = f36(i) + fac*vz6(i)
          f37(i) = f37(i) + fac*vz7(i)
          f38(i) = f38(i) + fac*vz8(i)
      ENDDO
      RETURN