multifluid_global_tdet.F File Reference

#include "my_real.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "param_c.inc"
#include "vect01_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	multifluid_global_tdet (iparg, elbuf_tab, multi_fvm, ipm)

Function/Subroutine Documentation

multifluid_global_tdet()

subroutine multifluid_global_tdet	(	integer, dimension(nparg,ngroup), intent(in)	iparg,
		type(elbuf_struct_), dimension(ngroup), intent(inout), target	elbuf_tab,
		type(multi_fvm_struct), intent(inout)	multi_fvm,
		integer, dimension(npropmi,nummat), intent(in)	ipm )

Definition at line 30 of file multifluid_global_tdet.F.

C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C CASE OF SUBMATERIAL BASED ON JWL MATERIAL LAW (HIGH EXPLISIVE WITH TIME CONTROL DETONATION)
C
C This subroutine is computing detonation times in global buffer from detonation times of layers (submaterials
C Global detonation time is relevant when there are several explosive submaterial laws.
C
C GBUF%TB : global burning time (may not be allocated in case of JWL is not used, check GBUF%G_TB > 0)
C LBUF%TB : local burning time of current layer
C If initial volume fraction is 0.0 then detonation time of current layer is not taken into account to calculate global detonation time.
C
C This subroutine must be called after INIVIOL treatment (INITIAL VOLUME FRACTION)
C Precondition : MULTI_FVM%IS_USED=.TRUE.
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod 
      USE multi_fvm_mod
      USE message_mod
      USE multi_fvm_mod
      USE constant_mod , ONLY : zero, ep21
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
          implicit none
#include      "my_real.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "param_c.inc"
#include      "vect01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER,INTENT(IN) :: IPARG(NPARG,NGROUP),IPM(NPROPMI,NUMMAT)
      TYPE(ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP),INTENT(INOUT) :: ELBUF_TAB
      TYPE(MULTI_FVM_STRUCT), INTENT(INOUT) :: MULTI_FVM     
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      LOGICAL IS_JWL
      INTEGER :: NEL,IG,OFFSET,NG,SUBMATLAW,ILAY,NLAY,II
      my_real :: vfrac
      TYPE(L_BUFEL_) ,POINTER :: LBUF     
      TYPE(G_BUFEL_) ,POINTER :: GBUF  
      TYPE(BUF_MAT_) ,POINTER :: MBUF
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
      IF(multi_fvm%IS_USED)THEN
        DO ng=1,ngroup    
              mtn     = iparg(1,ng)
              nel     = iparg(2,ng)
              nft     = iparg(3,ng)
              ity     = iparg(5,ng)
              gbuf => elbuf_tab(ng)%GBUF
              !---skip if Burn Fraction is not allocated
              IF(gbuf%G_TB > 0)THEN                
                !---skip if material law is not #151
                IF (mtn == 151) THEN
                  nlay = elbuf_tab(ng)%NLAY
                  DO ig=1,nel,nvsiz                                                                               
                    offset   = ig - 1                                                                                
                    lft      = 1                                             
                    llt      = min(nvsiz,nel-offset)                                                                 
                    nft      = iparg(3,ng) + offset  
                    is_jwl = .false.
                    !Number of layers ( = number of material in law 151)
                    IF (nlay > 1) THEN
                       gbuf%TB(lft:llt) = -ep21
                       DO ilay = 1, nlay
                          !SUBMATLAW = IPM(2, IPM(20 + ILAY, MTN))
                          submatlaw = elbuf_tab(ng)%BUFLY(ilay)%ILAW
                          IF (submatlaw == 5) THEN
                             lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(1,1,1)
                             DO ii = lft, llt
                                vfrac = lbuf%VOL(ii)/gbuf%VOL(ii)
                                IF(vfrac >= zero)THEN
                                  gbuf%TB(ii) = max(gbuf%TB(ii), lbuf%TB(ii))
                                ENDIF
                             ENDDO
                             is_jwl = .true.
                          ENDIF
                       ENDDO
                       IF(.NOT.is_jwl)THEN
                         DO ii = lft, llt
                           gbuf%TB(ii) = zero
                         ENDDO   
                       ELSE
                         DO ii = lft, llt
                           IF(gbuf%TB(ii) <= -ep21)THEN
                             gbuf%TB(ii) = zero
                           ENDIF
                         ENDDO                                        
                       ENDIF
                     ENDIF  ! IF(NLAY > 1)
                  ENDDO ! next IG                          
                ENDIF ! IF(MTN == 151)                                      
              ENDIF ! IF(GBUF%G_TB > 0)
         ENDDO ! next NG                                 
 
      ENDIF ! IF(MULTI_FVM%IS_USED)THEN