dtsph.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "vect01_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	dtsph (ssp, pm, geo, pid, mat, rho0, vis, deltax, vol, dtx)

Function/Subroutine Documentation

dtsph()

subroutine dtsph	(		ssp,
			pm,
			geo,
		integer, dimension(*)	pid,
		integer, dimension(*)	mat,
			rho0,
			vis,
			deltax,
			vol,
			dtx )

Definition at line 41 of file dtsph.F.

C------------------------
C     CALCUL DT ELEMENT
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   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "vect01_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real
     .     ssp(*), pm(npropm,*),geo(npropg,*), rho0(*), vis(*), 
     .     deltax(*), vol(*), dtx(*)
      INTEGER PID(*),MAT(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .   cx(mvsiz), qx(mvsiz), al(mvsiz), 
     .    qa, qb, cns1, cns2
C=======================================================================
      DO i=lft,llt
        cx(i)=ssp(i)
      END DO
 
      DO i=lft,llt
      al(i)=vol(i)**third
      END DO
 
      DO i=lft,llt
        qa  =geo(14,pid(i))
        qb  =geo(15,pid(i))
        cns1=geo(16,pid(i))
        cns2=geo(17,pid(i))*ssp(i)*al(i)*rho0(i)
        qx(i)=(qb+cns1)*ssp(i)+(two*vis(i)+cns2)/(rho0(i)*deltax(i))
      ENDDO
 
      DO i=lft,llt
         dtx(i)=deltax(i)/
     .          max(em20,qx(i)+sqrt(qx(i)*qx(i)+cx(i)*cx(i)))
      END DO
c-----------
      RETURN