gray10.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	gray10 (off, snew, qold, eint, voln, dvol, sold1, sold2, sold3, sold4, sold5, sold6, d1, d2, d3, d4, d5, d6, nel, pold, einc)

Function/Subroutine Documentation

gray10()

subroutine gray10	(		off,
			snew,
			qold,
			eint,
			voln,
			dvol,
			sold1,
			sold2,
			sold3,
			sold4,
			sold5,
			sold6,
			d1,
			d2,
			d3,
			d4,
			d5,
			d6,
		integer	nel,
			pold,
			einc )

Definition at line 28 of file gray10.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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL
      my_real
     .   off(*),    snew(nel,6), qold(*),  eint(*),  voln(*),
     .   dvol(*),   sold1(*),    sold2(*), sold3(*),
     .   sold4(*),  sold5(*),    sold6(*), d1(*),    d2(*),  d3(*),
     .   d4(*),     d5(*),       d6(*),    pold(*),  einc(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com08_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   vavg(mvsiz), e1(mvsiz), e2(mvsiz), e3(mvsiz), e4(mvsiz),
     .   e5(mvsiz),   e6(mvsiz), dta
C-----------------------------------------------
     
      dta =fourth*dt1
 
      DO i=1,nel
        e1(i)=d1(i)*(sold1(i)+snew(i,1)+pold(i))
        e2(i)=d2(i)*(sold2(i)+snew(i,2)+pold(i))
        e3(i)=d3(i)*(sold3(i)+snew(i,3)+pold(i))
        e4(i)=d4(i)*(sold4(i)+snew(i,4))
        e5(i)=d5(i)*(sold5(i)+snew(i,5))
        e6(i)=d6(i)*(sold6(i)+snew(i,6))
      ENDDO
 
      DO i=1,nel
        vavg(i)=(two*voln(i)-dvol(i))*off(i)
        dvol(i)=dvol(i)*off(i)
        einc(i)=vavg(i)*(e1(i)+e2(i)+e3(i)+e4(i)+e5(i)+e6(i))*dta
      ENDDO
 
      DO i=1,nel
        einc(i)=einc(i)-dvol(i)*(pold(i)+qold(i))*half
        eint(i)=eint(i)+einc(i)
      ENDDO
 
      RETURN