dtel.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	dtel (ssp, pm, geo, pid, mat, rho0, vis, deltax, aire, vol, dtx)
subroutine	dt70el (ssp, pm, geo, pid, mat, rho0, vis, deltax, aire, vol, dtx, igeo)

Function/Subroutine Documentation

dt70el()

subroutine dt70el	(		ssp,
			pm,
			geo,
		integer, dimension(*)	pid,
		integer, dimension(*)	mat,
			rho0,
			vis,
			deltax,
			aire,
			vol,
			dtx,
		integer, dimension(npropgi,*)	igeo )

Definition at line 107 of file dtel.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      "com01_c.inc"
#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(*), aire(*), vol(*), dtx(*)
      INTEGER PID(*),MAT(*),IGEO(NPROPGI,*)
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) 
      my_real
     .   qa, qb, cns1, cns2,facq
C-----------------------------------------------
C=======================================================================
       facq = one
!! qa=qb=zero by default for law70     
      IF(igeo(31,pid(1)) == 1) facq = zero
      DO i=lft,llt
        cx(i)=ssp(i)
        al(i)=vol(i)** third
        IF(n2d>0) al(i)=sqrt(aire(i))
      ENDDO
C
      DO i=lft,llt                                               
        qa =facq*geo(14,pid(i))                                       
        qb =facq*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)*max(em20,deltax(i))) 
      ENDDO                                                      
C
      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

dtel()

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

Definition at line 43 of file dtel.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      "com01_c.inc"
#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(*), aire(*), 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) 
      my_real
     .   qa, qb, cns1, cns2
C-----------------------------------------------
C=======================================================================
      DO i=lft,llt
        cx(i)=ssp(i)
        al(i)=vol(i)** third
        IF(n2d>0) al(i)=sqrt(aire(i))
      ENDDO
C
      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)*max(em20,deltax(i))) 
      ENDDO                                                      
C
      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