read_detonators.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    read_detonators        ../starter/source/initial_conditions/detonation/read_detonators.F
!||--- called by ------------------------------------------------------
!||    lectur                 ../starter/source/starter/lectur.F
!||--- calls      -----------------------------------------------------
!||    ancmsg                 ../starter/source/output/message/message.F
!||    hm_option_count        ../starter/source/devtools/hm_reader/hm_option_count.F
!||    hm_option_read_key     ../starter/source/devtools/hm_reader/hm_option_read_key.F
!||    hm_option_start        ../starter/source/devtools/hm_reader/hm_option_start.F
!||    prelecdet              ../starter/source/initial_conditions/detonation/prelecdet.F
!||    read_dfs_detcord       ../starter/source/initial_conditions/detonation/read_dfs_detcord.F
!||    read_dfs_detline       ../starter/source/initial_conditions/detonation/read_dfs_detline.F
!||    read_dfs_detplan       ../starter/source/initial_conditions/detonation/read_dfs_detplan.F
!||    read_dfs_detpoint      ../starter/source/initial_conditions/detonation/read_dfs_detpoint.F
!||    read_dfs_wave_shaper   ../starter/source/initial_conditions/detonation/read_dfs_wave_shaper.F
!||--- uses       -----------------------------------------------------
!||    detonators_mod         ../starter/share/modules1/detonators_mod.F
!||    hm_option_read_mod     ../starter/share/modules1/hm_option_read_mod.F
!||    message_mod            ../starter/share/message_module/message_mod.F
!||    submodel_mod           ../starter/share/modules1/submodel_mod.F
!||====================================================================
      SUBROUTINE read_detonators(ITABM1     ,ITAB      ,IGRNOD ,
     .                           PM         ,IPM       ,X      ,UNITAB,
     .                           LSUBMODEL  ,DETONATORS)
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 , ONLY : detonators_struct_
      USE names_and_titles_mod , ONLY : 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),ITAB(NUMNOD)
      my_real,INTENT(INOUT)                         :: pm(npropm,nummat)
      INTEGER,INTENT(IN)                            :: IPM(NPROPMI,NUMMAT)
      my_real,INTENT(IN)                            :: x(3,numnod)
      TYPE(submodel_data),INTENT(IN)                :: LSUBMODEL(NSUBMOD)
      TYPE(detonators_struct_),TARGET,INTENT(INOUT) :: 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              :: ID,K
      INTEGER              :: IUNIT, UID
      INTEGER              :: IFLAGUNIT
      INTEGER              :: DET_ID,STAT
      CHARACTER*40         :: MESS
      CHARACTER(LEN=NCHARKEY) :: KEY, KEY2
      INTEGER              :: NDET,NDFS
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      DATA mess/'DETONATORS DEFINITION                   '/
C-----------------------------------------------
C   P r e - C o n d i t i o n
C-----------------------------------------------
      ndet = detonators%N_DET
      detonators%IS_SHADOWING_REQUIRED = .false.
      IF (ndet <= 0) RETURN
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
 
      !---READING DFS OPTIONS & CHECK UNIT ID
      CALL hm_option_count('/DFS', ndfs)
      CALL hm_option_start('/DFS')
      WRITE(iout,1000)
      DO k=1,ndfs
        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)
        !CHECK UNIT SYSTEM ID (IF DEFINED)
        IF(uid > 0)THEN
          iflagunit=0
          DO iunit=1,unitab%NUNITS
            IF (unitab%UNIT_ID(iunit) == uid) THEN
              iflagunit = 1
              EXIT
            ENDIF
          ENDDO
          IF (iflagunit == 0) THEN
            CALL ancmsg(msgid=659,anmode=aninfo,msgtype=msgerror,
     .                  i2=uid,i1=det_id,c1='DETONATORS',c2='DETONATORS',c3='/DFS')
          ENDIF
        ENDIF
      ENDDO !K=1,NDFS
 
      !---SIZING & ALLOCATION
      CALL prelecdet(igrnod, lsubmodel, detonators)
 
      ! --- READING ALL DETONATOR TYPES
      CALL read_dfs_detpoint(detonators,x,ipm,pm,itabm1,unitab,lsubmodel,igrnod)
      CALL read_dfs_detline(detonators,x,ipm,itabm1,unitab,lsubmodel)
      CALL read_dfs_wave_shaper(detonators,igrnod,ipm,itabm1,unitab,lsubmodel,itab)
      CALL read_dfs_detplan(detonators,x,ipm,itabm1,unitab,lsubmodel)
      CALL read_dfs_detcord(detonators,x,igrnod,ipm,itabm1,unitab,lsubmodel,itab)
      WRITE(iout,2000)
C-----------------------------------------------
C   O u t p u t   F o r m a t
C-----------------------------------------------
 
 1000 FORMAT(////'    DETONATORS' /
     .           '    ----------'/)
 2000 FORMAT(//)
      END