unused_mat_detonator.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
integer function	unused_mat_detonator (mdet, nummat, listmat)

Function/Subroutine Documentation

unused_mat_detonator()

integer function unused_mat_detonator	(	integer, intent(inout)	mdet,
		integer, intent(in)	nummat,
		integer, dimension(npropmi,*), intent(in)	listmat )

Definition at line 32 of file unused_mat_detonator.F.

C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C This subroutine is checking that detonator material id
C does exist in input file
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      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER               :: UNUSED_MAT_DETONATOR
        INTEGER,intent(in)    :: LISTMAT(NPROPMI,*)     !material used in the model : {LISTMAT(1,I)},I=1..NUMMAT
        INTEGER,intent(inout) :: MDET                   !detonator flag
        INTEGER,intent(in)    :: NUMMAT                 !Number of Material cards defined in user file
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER JJ, I, MID, MLN
C-----------------------------------------------
          unused_mat_detonator=0
          jj=0
          IF (mdet /= 0) THEN
            DO i=1,nummat
              mid= listmat(1,i)
              mln= listmat(2,i)
              IF (mdet == mid) THEN !MDET refer to an used material
                IF(mln /= 5 .AND. mln /= 51 .AND. mln /= 97 .AND. mln /= 151)THEN
                  unused_mat_detonator=2
                  RETURN
                ENDIF
                mdet=i !Internal ID
                EXIT
              ELSEIF(mid == 0)THEN  !end of list reached
                jj=jj+1
                EXIT
              ELSE                  !MDET does not correspond to this current material
                jj=jj+1
              END IF
            END DO
            !!JJ=I => All material were compared to Mdet and do not correspond to it.
            IF (jj==i) THEN
              unused_mat_detonator=1
            END IF
          END IF