prelecdet.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/.
!||====================================================================
!||    prelecdet              ../starter/source/initial_conditions/detonation/prelecdet.F
!||--- called by ------------------------------------------------------
!||    read_detonators        ../starter/source/initial_conditions/detonation/read_detonators.F
!||--- calls      -----------------------------------------------------
!||    hm_get_intv            ../starter/source/devtools/hm_reader/hm_get_intv.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
!||    ngr2usrn               ../starter/source/system/nintrr.F
!||--- uses       -----------------------------------------------------
!||    detonators_mod         ../starter/share/modules1/detonators_mod.F
!||    hm_option_read_mod     ../starter/share/modules1/hm_option_read_mod.F
!||    submodel_mod           ../starter/share/modules1/submodel_mod.F
!||====================================================================
      SUBROUTINE prelecdet(IGRNOD, LSUBMODEL, DETONATORS)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE groupdef_mod
      USE detonators_mod
      USE submodel_mod         
      USE hm_option_read_mod
      USE detonators_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"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE(submodel_data),INTENT(IN) :: LSUBMODEL(NSUBMOD)
      TYPE(detonators_struct_) :: DETONATORS
C-----------------------------------------------
      TYPE (GROUP_)  , DIMENSION(NGRNOD)  :: IGRNOD
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER K, ID, USER_ID, INTERNAL_ID
      INTEGER NNOD_DETCORD, NNOD_WAV_SHA,NUM
      CHARACTER(LEN=NCHARKEY) :: KEY,KEY2,KEY3
      CHARACTER(LEN=NCHARTITLE) :: TITR
      LOGICAL :: IS_AVAILABLE
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER  NGR2USRN
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------      
 
! Allocation for /DFS/WAV_SHA & /DFS/DETCORD
 
      !ALLOCATE NUMBER OF DETONATORS IN DATA STRUCTURE
      detonators%N_DET = 0 !total number of detonators
      !
      num = detonators%N_DET_POINT
      detonators%N_DET = detonators%N_DET + num
      IF(num > 0) ALLOCATE(detonators%POINT(num))
      !
      num = detonators%N_DET_LINE
      detonators%N_DET = detonators%N_DET + num
      IF(num > 0) ALLOCATE(detonators%LINE(num))
      !
      num = detonators%N_DET_PLANE
      detonators%N_DET = detonators%N_DET + num
      IF(num > 0) ALLOCATE(detonators%PLANE(num))
      !
      num = detonators%N_DET_WAVE_SHAPER
      detonators%N_DET = detonators%N_DET + num
      IF(num > 0) ALLOCATE(detonators%WAVE_SHAPER(num))
      !
      num = detonators%N_DET_CORD
      detonators%N_DET = detonators%N_DET + num
      IF(num > 0) ALLOCATE(detonators%CORD(num))
 
      !---SPECIFIC ALLOCATIONS (WAV_SHA & DETCORD)
      !- /dfs/wav_sha
      CALL hm_option_start('/DFS/WAV_SHA')
      DO k=1,detonators%N_DET_WAVE_SHAPER
        CALL hm_option_read_key(lsubmodel,option_id=id,keyword2= key)
        CALL hm_get_intv('entityid', user_id, is_available, lsubmodel)
        internal_id = ngr2usrn(user_id,igrnod,ngrnod,nnod_wav_sha)
        detonators%WAVE_SHAPER(k)%NUMNOD = nnod_wav_sha
        ALLOCATE(detonators%WAVE_SHAPER(k)%NODES(nnod_wav_sha))
        ALLOCATE(detonators%WAVE_SHAPER(k)%IORDR(nnod_wav_sha))
        ALLOCATE(detonators%WAVE_SHAPER(k)%FLAG(nnod_wav_sha))
        ALLOCATE(detonators%WAVE_SHAPER(k)%TIME(nnod_wav_sha))
        detonators%WAVE_SHAPER(k)%NODES(:) = 0
        detonators%WAVE_SHAPER(k)%IORDR(:) = 0
        detonators%WAVE_SHAPER(k)%FLAG(:) = 0 
        detonators%WAVE_SHAPER(k)%TIME(:) = zero
      END DO
 
      CALL hm_option_start('/DFS/DETCORD')
      DO k=1,detonators%N_DET_CORD
        CALL hm_option_read_key(lsubmodel,option_id=id,keyword2= key)
        CALL hm_get_intv('entityid', user_id, is_available, lsubmodel)
        internal_id = ngr2usrn(user_id,igrnod,ngrnod,nnod_detcord)
        detonators%CORD(k)%NUMNOD = nnod_detcord
        ALLOCATE(detonators%CORD(k)%NODES(nnod_detcord))
        ALLOCATE(detonators%CORD(k)%TDET_PATH(nnod_detcord))
        detonators%CORD(k)%NODES(:) = 0
        detonators%CORD(k)%TDET_PATH(:) = zero
      ENDDO
 
 
 
 
C-------------------------------------------
      RETURN
      END SUBROUTINE prelecdet