sigeps37.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine

sigeps37 (nel, nuparam, nuvar, nfunc, ifunc, npf, tf, time, timestep, uparam, rho0, rho, volume, eint, epspxx, epspyy, epspzz, epspxy, epspyz, epspzx, depsxx, depsyy, depszz, depsxy, depsyz, depszx, epsxx, epsyy, epszz, epsxy, epsyz, epszx, sigoxx, sigoyy, sigozz, sigoxy, sigoyz, sigozx, signxx, signyy, signzz, signxy, signyz, signzx, sigvxx, sigvyy, sigvzz, sigvxy, sigvyz, sigvzx, soundsp, viscmax, uvar, off)

Function/Subroutine Documentation

sigeps37()

subroutine sigeps37	(	integer	nel,
		integer	nuparam,
		integer	nuvar,
		integer	nfunc,
		integer, dimension(nfunc)	ifunc,
		integer, dimension(*)	npf,
			tf,
			time,
			timestep,
			uparam,
			rho0,
			rho,
			volume,
			eint,
			epspxx,
			epspyy,
			epspzz,
			epspxy,
			epspyz,
			epspzx,
			depsxx,
			depsyy,
			depszz,
			depsxy,
			depsyz,
			depszx,
			epsxx,
			epsyy,
			epszz,
			epsxy,
			epsyz,
			epszx,
			sigoxx,
			sigoyy,
			sigozz,
			sigoxy,
			sigoyz,
			sigozx,
			signxx,
			signyy,
			signzz,
			signxy,
			signyz,
			signzx,
			sigvxx,
			sigvyy,
			sigvzz,
			sigvxy,
			sigvyz,
			sigvzx,
			soundsp,
			viscmax,
			uvar,
			off )

Definition at line 30 of file sigeps37.F.

C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C This subroutines aim is to compute initial sound speed
C-----------------------------------------------
C   I M P L I C I T   T Y P E S
C-----------------------------------------------
#include      "implicit_f.inc"
C----------------------------------------------------------------
C  I N P U T   A R G U M E N T S
C----------------------------------------------------------------
      INTEGER       NEL,     NUPARAM, NUVAR
 
      my_real
     .      time       , timestep   , uparam(nuparam),
     .      rho(nel), rho0(nel), volume(nel), eint(nel),
     .      epspxx(nel), epspyy(nel), epspzz(nel),
     .      epspxy(nel), epspyz(nel), epspzx(nel),
     .      depsxx(nel), depsyy(nel), depszz(nel),
     .      depsxy(nel), depsyz(nel), depszx(nel),
     .      epsxx(nel), epsyy(nel), epszz(nel),
     .      epsxy(nel), epsyz(nel), epszx(nel),
     .      sigoxx(nel), sigoyy(nel), sigozz(nel),
     .      sigoxy(nel), sigoyz(nel), sigozx(nel)
C----------------------------------------------------------------
C  O U T P U T   A R G U M E N T S
C----------------------------------------------------------------
      my_real
     .      signxx(nel), signyy(nel), signzz(nel),
     .      signxy(nel), signyz(nel), signzx(nel),
     .      sigvxx(nel), sigvyy(nel), sigvzz(nel),
     .      sigvxy(nel), sigvyz(nel), sigvzx(nel),
     .      soundsp(nel), viscmax(nel)
C----------------------------------------------------------------
C  I N P U T  O U T P U T   A R G U M E N T S
C----------------------------------------------------------------
      my_real uvar(nel,nuvar), off(nel) 
C----------------------------------------------------------------
C  VARIABLES FOR FUNCTION INTERPOLATION 
C----------------------------------------------------------------
      INTEGER NPF(*), NFUNC, IFUNC(NFUNC)
      my_real
     . finter,tf(*)
      EXTERNAL finter
C----------------------------------------------------------------
C  L O C A L  V A R I B L E S
C----------------------------------------------------------------
      INTEGER :: I
      my_real :: c1,gam,pmin,p0,rho10,rho20,a1,p,mu1p1,mu2p1,rho1,rho2,a
 
 
C----------------------------------------------------------------
      DO i=1,nel
        c1           =  uparam(4)
        gam          =  uparam(5)
        pmin         =  uparam(8)
        p0           =  uparam(9)
        rho10        =  uparam(11)
        rho20        =  uparam(12)
        a1           =  uparam(10)         
        p            = -p0  
        rho(i)       =  rho10 * a1 + (one-a1)*rho20     
        IF(gam*c1>=em30)THEN !if Liquid and gas correctly defined
          mu1p1      = one  !(P/C1+ONE)
          mu2p1      = one  !( ONE+P/P0)**(ONE/GAM)
          rho1       = rho10*mu1p1
          rho2       = rho20*mu2p1
          a          = (rho(i)-rho2)/(rho1-rho2)
          soundsp(i) = a*sqrt(c1/rho1) + (one-a) * sqrt(gam*p0/rho2) 
        ENDIF 
         
      ENDDO
 
 
      RETURN