read_dfs_wave_shaper.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "units_c.inc"
#include "param_c.inc"
#include "tabsiz_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	read_dfs_wave_shaper (detonators, igrnod, ipm, itabm1, unitab, lsubmodel, itab)

Function/Subroutine Documentation

read_dfs_wave_shaper()

subroutine read_dfs_wave_shaper	(	type(detonators_struct_), intent(inout), target	detonators,
		type (group_), dimension(ngrnod), intent(in)	igrnod,
		integer, dimension(npropmi,nummat), intent(in)	ipm,
		integer, dimension(sitabm1), intent(in)	itabm1,
		type (unit_type_), intent(in)	unitab,
		type(submodel_data), dimension(nsubmod), intent(in)	lsubmodel,
		integer, dimension(numnod), intent(in)	itab )

Definition at line 46 of file read_dfs_wave_shaper.F.

C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE unitab_mod
      USE message_mod
      USE detonators_mod
      USE groupdef_mod
      USE hm_option_read_mod
      USE submodel_mod
      USE detonators_mod
      USE names_and_titles_mod , ONLY : nchartitle, ncharkey
      USE format_mod , ONLY : fmw_10i
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      "com04_c.inc"
#include      "units_c.inc"
#include      "param_c.inc"
#include      "tabsiz_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE (UNIT_TYPE_),INTENT(IN)                  :: UNITAB
      INTEGER,INTENT(IN)                            :: ITABM1(SITABM1),ITAB(NUMNOD)
      INTEGER,INTENT(IN)                            :: IPM(NPROPMI,NUMMAT)
      TYPE(SUBMODEL_DATA),INTENT(IN)                :: LSUBMODEL(NSUBMOD)
      TYPE(DETONATORS_STRUCT_),INTENT(INOUT),TARGET :: DETONATORS
      TYPE (GROUP_),DIMENSION(NGRNOD),INTENT(IN)    :: IGRNOD
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER              :: I, MAT, J, K,IGU,IGS,JJ,MDET,DET_ID,IDET
      INTEGER              :: IBID,  NODE_ID1, NODE_ID2,uID1,uID2, IOPT, IUNIT, UID
      INTEGER              :: FLAG_FMT,IMAT,IFLAGUNIT,UNUSED
      INTEGER              :: STAT,NPE
      my_real              :: xc, yc, zc, alt, xc1, yc1, zc1, xc2, yc2, zc2, nx, ny, nz, bid, vcj, vdet
      CHARACTER*40         :: MESS
      CHARACTER*64         :: chain1,chain2
      CHARACTER(LEN=NCHARKEY) :: KEY, KEY2
      CHARACTER(LEN=NCHARTITLE) :: TITR
      LOGICAL              :: IS_ENCRYPTED, IS_AVAILABLE
      INTEGER,POINTER,DIMENSION(:) :: IECRAN
      my_real,POINTER,DIMENSION(:) :: decran
      INTEGER              :: NDETPS,NDETSG,NECRAN,NDETPL,NDETCORD
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER,EXTERNAL  :: NODGRNR5, USR2SYS, NINTRI
      INTEGER           :: UNUSED_MAT_DETONATOR
      DATA mess/'DETONATORS DEFINITION                   '/      
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
 
      CALL hm_option_start('/DFS/WAV_SHA')
 
      DO idet=1,detonators%N_DET_WAVE_SHAPER
 
        CALL hm_option_read_key(lsubmodel,option_id=det_id, unit_id=uid,keyword2=key,keyword3=key2)
        IF (len_trim(key) > 0) key = key(1:7)
        IF (len_trim(key2) > 0) key2 = key2(1:4)
 
        iecran => detonators%WAVE_SHAPER(idet)%NODES(1:)
        decran => detonators%WAVE_SHAPER(idet)%TIME(1:)
 
        is_encrypted= .false.
        is_available = .false.
        CALL hm_option_is_encrypted(is_encrypted)
        !---------------------------------!
        !             READING             !
        !---------------------------------!
        CALL hm_get_floatv('rad_det_locationA_X', vdet, is_available, lsubmodel, unitab)
        CALL hm_get_floatv('rad_det_locationa_y', YC1, IS_AVAILABLE, LSUBMODEL, UNITAB)
        CALL HM_GET_FLOATV('rad_det_locationa_z', ZC1, IS_AVAILABLE, LSUBMODEL, UNITAB)
        CALL HM_GET_FLOATV('rad_det_time', ALT, IS_AVAILABLE, LSUBMODEL,UNITAB)
        CALL HM_GET_INTV('rad_det_materialid', MAT, IS_AVAILABLE, LSUBMODEL)
        CALL HM_GET_INTV('entityid', IGU, IS_AVAILABLE, LSUBMODEL)
        !---------------------------------!
        !      CHECKING USER FLAGS        !
        !      +INTERNAL ID               !
        !---------------------------------!
        MDET=MAT !bak
        IF (ALT > INFINITY) ALT=INFINITY
        IF(MAT > 0)UNUSED=UNUSED_MAT_DETONATOR(MAT,NUMMAT,IPM)
        UNUSED=0
        IF (MAT < 0) THEN
          CALL ANCMSG(MSGID=102,MSGTYPE=MSGERROR,ANMODE=ANINFO,
     .                I1=DET_ID,
     .                C1='detonator is referring to a negative material id',
     .                C2='/dfs/wav_sha',
     .                I2=MDET)
        ELSEIF (UNUSED==1) THEN
          CALL ANCMSG(MSGID=102,MSGTYPE=MSGERROR,ANMODE=ANINFO,
     .                I1=DET_ID,
     .                C1='detonator is referring to an unknown material id',
     .                C2='/dfs/wav_sha',
     .                I2=MDET)
        ELSEIF (UNUSED==2) THEN
          CALL ANCMSG(MSGID=102,MSGTYPE=MSGERROR,ANMODE=ANINFO,
     .                I1=DET_ID,
     .                C1='detonator must refer to a jwl material law(laws 5, 51, 97, 151)',
     .                C2='/dfs/wav_sha',
     .                I2=MDET)
        ELSE
          !---------------------------------!
          !        LISTING PRINTOUT         !
          !---------------------------------!
          NPE = NODGRNR5(IGU ,IGS,IECRAN(1),IGRNOD ,ITABM1 ,MESS)
          !Nodes in group are ordered from 1 to NPE. last point is the nearest from the detonation origin. Points are composing the screen lines (guard lines). Screen lines is the boundary of the obstacle (shadow area
          IF(IS_ENCRYPTED) WRITE(IOUT,1001)
.NOT.          IF(IS_ENCRYPTED)WRITE(IOUT,1550) DET_ID,VDET,YC1,ZC1,ALT,MDET,IGU,NPE
.NOT.          IF(IS_ENCRYPTED)WRITE(IOUT,FMT=FMW_10I) (ITAB(IECRAN(I)),I=1,NPE)
          DO I=1,NPE
            CALL IFRONTPLUS(IECRAN(I),1)
          END DO
          DETONATORS%WAVE_SHAPER(IDET)%TDET = ALT
          DETONATORS%WAVE_SHAPER(IDET)%MAT = MAT
          DETONATORS%WAVE_SHAPER(IDET)%VDET = VDET
          DETONATORS%WAVE_SHAPER(IDET)%XDET = ZERO
          DETONATORS%WAVE_SHAPER(IDET)%YDET = YC1
          DETONATORS%WAVE_SHAPER(IDET)%ZDET = ZC1
          DETONATORS%WAVE_SHAPER(IDET)%NUMNOD = NPE
        END IF
 
      ENDDO !next IDET
C-----------------------------------------------
C   O u t p u t   F o r m a t
C-----------------------------------------------
 1001 FORMAT(///5X,
     &          'shadow line detonation    ',I10,/5X,
     &          '----------------------    ',/5X,
     &          'confidential data')    
 1550 FORMAT(///5X,
     &          'shadow line detonation    =',I10,/5X,
     &          '----------------------     ',/5X,
     &          'OPTIONAL velocity         =',1PG20.13,/5X,
     &          'y-coordinate              =',1PG20.13,/5X,
     &          'z-coordinate              =',1PG20.13,/5X,
     &          'lighting time             =',1PG20.13,/5X,
     &          'explosive material number =',I10,/5X,
     &          'shadow line node group id =',I10,/5X,
     &          'number of points(shadow)  =',I10,/5X,
     &          'shadow line definition :   ')
 
C-----------------------------------------------