meint.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	meint (output, off, snew, qold, eint, vnew, espe, sold1, sold2, sold3, sold4, sold5, sold6, d1, d2, d3, d4, d5, d6, psh, dvol, df, qnew, pnew, vis, tmu, einc, mln, vol_avg, nel, jtur, jlag, jpor)

Function/Subroutine Documentation

meint()

subroutine meint	(	type(output_), intent(inout)	output,
			off,
			snew,
			qold,
			eint,
			vnew,
			espe,
			sold1,
			sold2,
			sold3,
			sold4,
			sold5,
			sold6,
			d1,
			d2,
			d3,
			d4,
			d5,
			d6,
			psh,
			dvol,
			df,
			qnew,
			pnew,
			vis,
			tmu,
			einc,
		integer	mln,
			vol_avg,
		integer	nel,
		integer, intent(in)	jtur,
		integer, intent(in)	jlag,
		integer, intent(in)	jpor )

Definition at line 31 of file meint.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE ale_mod
      USE output_mod , ONLY : output_ 
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      "com06_c.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-----------------------------------------------
      TYPE(OUTPUT_), INTENT(INOUT) :: OUTPUT
      INTEGER, INTENT(IN) :: JTUR
      INTEGER, INTENT(IN) :: JLAG
      INTEGER, INTENT(IN) :: JPOR
      INTEGER MLN,NEL
      my_real off(*), snew(nel,6), qold(*), eint(*)
      my_real vnew(*), espe(*), qnew(*), pnew(*), sold1(*), sold2(*),
     .        sold3(*), sold4(*), sold5(*), sold6(*), 
     .        dvol(*), df(*), d1(*), d2(*), d3(*), d4(*), d5(*),
     .        d6(*), psh(*), einc(*),
     .        vis(*), tmu(*), vol_avg(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com08_c.inc"
#include      "scr06_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real pold(mvsiz), e1, e2, e3, e4, e5, e6, dta, wfextt
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      dta = half*dt1
 
      IF(ale%GLOBAL%INCOMP == 1 .AND. jlag == 0)THEN
        DO i=1,nel
          pold(i)=(sold1(i)+sold2(i)+sold3(i))*third
          eint(i)=eint(i)+(half*dvol(i)*(pold(i)-psh(i)))*off(i)
          espe(i)=eint(i)*df(i) / max(em15,vnew(i))
          qold(i)=qnew(i)
        ENDDO
        RETURN
      ENDIF
 
      DO i=1,nel
        pold(i)=(sold1(i)+sold2(i)+sold3(i)) * third
        sold1(i)=sold1(i)-pold(i)
        sold2(i)=sold2(i)-pold(i)
        sold3(i)=sold3(i)-pold(i)
        e1=d1(i)*(sold1(i)+snew(i,1))
        e2=d2(i)*(sold2(i)+snew(i,2))
        e3=d3(i)*(sold3(i)+snew(i,3))
        e4=d4(i)*(sold4(i)+snew(i,4))
        e5=d5(i)*(sold5(i)+snew(i,5))
        e6=d6(i)*(sold6(i)+snew(i,6))
        einc(i) = vol_avg(i) * (e1+e2+e3+e4+e5+e6) * dta
      ENDDO
 
      wfextt = zero
      DO i=1,nel
        wfextt = wfextt - dvol(i)*psh(i)
      ENDDO
!$OMP ATOMIC
      output%TH%WFEXT = output%TH%WFEXT + wfextt
 
C
      IF(jtur == 0 .OR. jpor == 2)THEN
       DO i=1,nel
         eint(i)=eint(i)+(einc(i)+half*dvol(i)*(pold(i)-psh(i)-qold(i)-qnew(i)))*off(i)
       ENDDO
      ELSE
       DO i=1,nel
         eint(i)=eint(i) + (einc(i)*(vis(i)-tmu(i))/vis(i) + half*dvol(i)*(pold(i)-psh(i)-qold(i)-qnew(i)))*off(i)
       ENDDO
      ENDIF
C
      DO i=1,nel
        qold(i)=qnew(i)
        espe(i)=eint(i)*df(i) / max(em15,vnew(i)) !ESPE = 'rho0.e'
      ENDDO
 
      !NO NEED TO UPDATE PRESSURE CALLING EOSMAIN SUBROUTINE FOR THESE TWO MATERIAL LAWS SINCE PRESSURE IS NOT ENERGY DEPENDANT.
      ! embedded pressure function
      IF(mln==21)THEN
        DO i=1,nel
          eint(i) = eint(i) - half*dvol(i)*(pnew(i)+psh(i))
        ENDDO
      ENDIF
C
      RETURN