hm_prelecgrns.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/.
!||====================================================================
!||    hm_prelecgrns            ../starter/source/groups/hm_prelecgrns.f
!||--- called by ------------------------------------------------------
!||    lectur                   ../starter/source/starter/lectur.F
!||--- calls      -----------------------------------------------------
!||    hm_get_int_array_index   ../starter/source/devtools/hm_reader/hm_get_int_array_index.F
!||    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
!||    usr2sys                  ../starter/source/system/sysfus.F
!||--- uses       -----------------------------------------------------
!||    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 hm_prelecgrns(ITABM1 ,IGRNOD ,LSUBMODEL)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE my_alloc_mod
      USE message_mod
      USE groupdef_mod
      USE submodel_mod
      USE hm_option_read_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-----------------------------------------------
      INTEGER ITABM1(*)
C-----------------------------------------------
      TYPE (GROUP_)  , DIMENSION(NGRNOD) :: IGRNOD
      TYPE(submodel_data) LSUBMODEL(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER J10(10)
      INTEGER I,J,K,ID,NNOD,NL,IGS,KK,JJ,NENTITY,UID,NN
      CHARACTER(LEN=NCHARTITLE) :: TITR,TITR1
      CHARACTER(LEN=NCHARKEY) :: KEY,KEY2
      CHARACTER MESS*40
      LOGICAL IS_AVAILABLE
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER USR2SYS
      DATA mess/'NODENS GROUP DEFINITION                 '/
C-----------------------------------------------
C    IGRNOD(IGS)%ID : GROUP identifier
C    IGRNOD(IGS)%TITLE : GROUP title
C    IGRNOD(IGS)%NENTITY : Entities (nodes) number of the GROUP
C    IGRNOD(IGS)%GRTYPE : TYPE ( 0-NOEUDS, 1-BRIC, 2-QUAD,   3-SHELL_4N,
!                                4-TRUSS,  5-BEAM, 6-SPRINGS,7-SHELL_3N )
!!                   GRTYPE --- > OBSOLETE
C    IGRNOD(IGS)%SORTED  : FLAG for sorted/unsorted nodes
!                    = 0 -> sorted
!                    = 1 -> unsorted
C    IGRNOD(IGS)%GRPGRP : TYPE of nodal GROUP
!                                 = 1 FOR /GRNOD/NOD
!                                 = 2 FOR /GRNOD/GNROD
!                                 = 0 ALL THE REST
C    IGRNOD(IGS)%LEVEL  : Hierarchy level
!                             (FLAG 'SUBLEVEL DONE' FOR GROUPS OF GROUPS)
!                            = 0 ---> not yet initialized
!                            = 1 ---> done
!  - R2R -
!                              ! R2R_ALL <--- IGROUP2(8,I) = IGROUP2(2,I) -- before splitting
!                              ! R2R_SHARE <--- IGROUP2(9,I)   (shared on boundary subdomain)
C    IGRNOD(IGS)%R2R_ALL   ! Multidomaines -> number of node(elems, parts) before split
C    IGRNOD(IGS)%R2R_SHARE ! shared on boundary subdomain
C-----------------------------------------------
      is_available = .false.
C=======================================================================
C groupes de noeuds(type NODE,NODENS) + init IGRNOD(IGS)%ID
C=======================================================================
      igs=0
      titr1='NODENS GROUP'
!
      CALL hm_option_start('/GRNOD')
!
C      
      ! Loop over all GRNOD
      DO i=1,ngrnod
C
        ! Read keys and uid
        CALL hm_option_read_key(lsubmodel,
     .                          option_id   = id,
     .                          option_titr = titr  ,
     .                          unit_id     = uid,
     .                          keyword2    = key   ,
     .                          keyword3    = key2)   
C
        ! Increasing the counter
        igs=igs+1
C
        ! Initialization of data structure
        igrnod(igs)%ID=id
        igrnod(igs)%GRTYPE=0  !!! obsolete
        igrnod(igs)%SORTED=0
        igrnod(igs)%LEVEL=1
        igrnod(igs)%TITLE=titr
        nn=0
        nnod=0
        !-----------------------------------------------------------------
        ! GRNOD from NODENS        
        IF(key(1:6) == 'NODENS') THEN
!         groupe de noeuds (non sorted)
          CALL hm_get_intv('idsmax' ,nentity,is_available,lsubmodel)
          DO kk = 1,nentity
            CALL hm_get_int_array_index ('ids',jj  ,kk,is_available,lsubmodel)
            IF (jj /= 0) THEN
              nnod = nnod + 1
            ENDIF
          ENDDO
          igrnod(igs)%NENTITY = nnod
          igrnod(igs)%GRPGRP  = 1
          CALL my_alloc(igrnod(igs)%ENTITY,nnod)
          igrnod(igs)%SORTED  = 1
          DO kk = 1,nentity
            CALL hm_get_int_array_index ('ids',jj  ,kk,is_available,lsubmodel)
            IF (jj /= 0) THEN
              nn = nn + 1
              igrnod(igs)%ENTITY(nn)=usr2sys(jj,itabm1,mess,id)
            ENDIF
          ENDDO          
        ENDIF ! IF(LINE(1:1) == '/')
      ENDDO
!-----
      RETURN
      END