read_dfs_detline.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_detline (detonators, x, ipm, itabm1, unitab, lsubmodel)

Function/Subroutine Documentation

read_dfs_detline()

subroutine read_dfs_detline	(	type(detonators_struct_), intent(inout)	detonators,
		dimension(3,numnod), intent(in)	x,
		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 )

Definition at line 44 of file read_dfs_detline.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 names_and_titles_mod , ONLY : nchartitle, ncharkey
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)
      INTEGER,INTENT(IN)                     :: IPM(NPROPMI,NUMMAT)
      my_real,INTENT(IN)                     :: x(3,numnod)
      TYPE(SUBMODEL_DATA),INTENT(IN)         :: LSUBMODEL(NSUBMOD)
      TYPE(DETONATORS_STRUCT_),INTENT(INOUT) :: DETONATORS
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER              :: I, MAT, J, NPEM,NPCM,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
      my_real              :: xc, yc, zc, alt, xc1, yc1, zc1, xc2, yc2, zc2, nx, ny, nz, bid, vcj
      CHARACTER*40         :: MESS
      CHARACTER*64         :: chain1,chain2
      CHARACTER(LEN=NCHARKEY) :: KEY, KEY2
      CHARACTER(LEN=NCHARTITLE) :: TITR
      LOGICAL              :: IS_ENCRYPTED, IS_AVAILABLE
      LOGICAL              :: IS_NODE_DEFINED
      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/DETLINE')
 
      DO idet=1,detonators%N_DET_LINE
 
        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)
 
        is_encrypted= .false.
        is_available = .false.
        is_node_defined = .false.
        IF(key2(1:4) == 'NODE')is_node_defined = .true.
        CALL hm_option_is_encrypted(is_encrypted)
        !---------------------------------!
        !             READING             !
        !---------------------------------!
        IF(is_node_defined)THEN
          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('rad_det_node1', uid1, is_available, lsubmodel)
          CALL hm_get_intv('rad_det_node2', uid2, is_available, lsubmodel)
          xc=zero
          yc=zero
          zc=zero
        ELSE
          CALL hm_get_floatv('rad_det_locationA_X', xc1, 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_locationB_X', xc2, is_available, lsubmodel, unitab)
          CALL hm_get_floatv('rad_det_locationB_Y', yc2, is_available, lsubmodel, unitab)
          CALL hm_get_floatv('rad_det_locationB_Z', zc2, 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)
        ENDIF
 
        ! COORDINATES
        IF(is_node_defined)THEN
          node_id1=usr2sys(uid1,itabm1,mess,det_id)
          IF(node_id1 > 0)THEN
            xc1 = x(1,node_id1)
            yc1 = x(2,node_id1)
            zc1 = x(3,node_id1)
          ENDIF
          node_id2=usr2sys(uid2,itabm1,mess,det_id)
          IF(node_id2 > 0)THEN
            xc2 = x(1,node_id2)
            yc2 = x(2,node_id2)
            zc2 = x(3,node_id2)
          ENDIF
          IF(node_id1 == 0 .AND. node_id2 == 0)THEN
            CALL ancmsg(msgid = 104,msgtype = msgerror,anmode  = aninfo,
     .                  c1='/DFS/DETLINE/NODE',
     .                  i1=det_id,
     .                  c2='INVALID NODE_ID')
          ENDIF
        ENDIF
        !---------------------------------!
        !      CHECKING USER FLAGS        !
        !      +INTERNAL ID               !
        !---------------------------------!
        mdet=mat !bak
        IF (alt >  infinity)alt= infinity
        IF (alt < -infinity)alt=-infinity
        unused=0
        IF(mat > 0)unused=unused_mat_detonator(mat,nummat,ipm)
        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/DETLINE',
     .                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/DETLINE',
     .                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/DETLINE',
     .                i2=mdet)
        ELSE
          detonators%LINE(idet)%IS_MAT_VALID = .true.
          !---------------------------------!
          !        LISTING PRINTOUT         !
          !---------------------------------!
          IF(is_node_defined)THEN
            IF(.NOT.is_encrypted)WRITE(iout,1501) det_id,node_id1,xc1,yc1,zc1,node_id2,xc2,yc2,zc2,alt,mdet
          ELSE
            IF(.NOT.is_encrypted)WRITE(iout,1500) det_id,xc1,yc1,zc1,xc2,yc2,zc2,alt,mdet
          ENDIF
          IF(is_encrypted)WRITE(iout,1001)
 
          detonators%LINE(idet)%TDET = alt
          detonators%LINE(idet)%MAT = mat
          detonators%LINE(idet)%XDET_1 = xc1
          detonators%LINE(idet)%YDET_1 = yc1
          detonators%LINE(idet)%ZDET_1 = zc1
          detonators%LINE(idet)%XDET_2 = xc2
          detonators%LINE(idet)%YDET_2 = yc2
          detonators%LINE(idet)%ZDET_2 = zc2
        END IF
 
      ENDDO
 
C-----------------------------------------------
C   O u t p u t   F o r m a t
C-----------------------------------------------      
 1001 FORMAT(///5x,
     &          'DETONATION LINE           ',i10,/5x,
     &          '---------------           ',/5x,
     &          'CONFIDENTIAL DATA')    
 1500 FORMAT(///5x,
     &          'DETONATION LINE            ',i10,/5x,
     &          '---------------            ',/5x,
     &          'X-COORDINATE FIRST POINT  =',1pg20.13,/5x,
     &          'y-coordinate first point  =',1PG20.13,/5X,
     &          'z-coordinate first point  =',1PG20.13,/5X,
     &          'x-coordinate second point =',1PG20.13,/5X,
     &          'y-coordinate second point =',1PG20.13,/5X,
     &          'z-coordinate second point =',1PG20.13,/5X,
     &          'lighting time             =',1PG20.13,/5X,
     &          'explosive material number =',I10 )
 1501 FORMAT(///5X,
     &          'detonation line            ',I10,/5X,
     &          '---------------            ',/5x,
     &          'FIRST NODE ID             =',i10     ,/5x,     
     &          ' X-COORDINATE FIRST POINT =',1pg20.13,/5x,
     &          ' Y-COORDINATE FIRST POINT =',1pg20.13,/5x,
     &          ' Z-COORDINATE FIRST POINT =',1pg20.13,/5x,
     &          'SECOND NODE ID            =',i10     ,/5x,        
     &          ' X-COORDINATE SECOND POINT=',1pg20.13,/5x,
     &          ' Y-COORDINATE SECOND POINT=',1pg20.13,/5x,
     &          ' Z-COORDINATE SECOND POINT=',1pg20.13,/5x,
     &          'LIGHTING TIME             =',1pg20.13,/5x,
     &          'EXPLOSIVE MATERIAL NUMBER =',i10 )     
C-----------------------------------------------