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

Function/Subroutine Documentation

sesa10()

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

Definition at line 28 of file sesa10.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(*), dvol(*),
     .   sold1(*), sold2(*), sold3(*), sold4(*), sold5(*),
     .   sold6(*), d1(*), d2(*), d3(*), d4(*), d5(*), d6(*), pold(*), df(*), 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,voln(mvsiz)
C-----------------------------------------------
      dta =fourth*dt1
C
      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
C
      DO i=1,nel
        vavg(i)=(two*voln(i)-dvol(i))*off(i)
        dvol(i)=half*dvol(i)*off(i)
        einc(i)=vavg(i)*(e1(i)+e2(i)+e3(i)+e4(i)+e5(i)+e6(i))*dta
      ENDDO
C
      DO i=1,nel
        einc(i)=einc(i)-dvol(i)*(pold(i)+qold(i))
        eint(i)=eint(i)+einc(i)
      ENDDO
C
      RETURN