alloc_subset_str.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/.
!||====================================================================
!||    alloc_subset_str   ../engine/source/groups/alloc_subset_str.F
!||--- called by ------------------------------------------------------
!||    subset_ini         ../engine/source/groups/subset_ini.F
!||--- calls      -----------------------------------------------------
!||    fretitl2           ../engine/source/input/freform.F
!||    read_i_c           ../common_source/tools/input_output/write_routtines.c
!||--- uses       -----------------------------------------------------
!||    groupdef_mod       ../common_source/modules/groupdef_mod.F
!||    message_mod        ../engine/share/message_module/message_mod.F
!||====================================================================
      SUBROUTINE alloc_subset_str(SUBSET)
C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE message_mod
      USE groupdef_mod
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      "scr17_c.inc"
#include      "tabsiz_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE (SUBSET_) , DIMENSION(NSUBS) :: SUBSET
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,ID,ISU,NCHILD,NBPART,NBTPART,ICHILD,IPART,IPARTTOT,
     .        PARENT,LEVEL,TH_PART,L_SUB,ERR,ITITLE(LTITR),THIAD,IVARTH,
     .        LEN
      CHARACTER(LEN=LTITR) ::  TITR,NEW_TITLE(NSUBS)
      INTEGER, ALLOCATABLE, DIMENSION (:)  ::  ISUB
!=======================================================================
      DO isu=1,nsubs
        CALL read_i_c(ititle,ltitr)
        CALL fretitl2(titr,ititle,ltitr)
        len = 0
        DO j=1,ltitr
          IF (titr(j:j) /= ' ') len = j
        ENDDO
        new_title(isu) = titr(1:len)
      ENDDO ! DO ISU=1,NSUBS
!!------  Allocate and read subset data from restart
      ALLOCATE (isub(lensubs) ,stat=err)
      CALL read_i_c(isub, lensubs)
!=======================================================================
!
      l_sub = 0
!
      DO isu=1,nsubs
        id       = isub(l_sub+1)
          l_sub = l_sub+1
        nchild   = isub(l_sub+1)
          l_sub = l_sub+1
        nbpart   = isub(l_sub+1)
          l_sub = l_sub+1
        nbtpart  = isub(l_sub+1)
          l_sub = l_sub+1
        parent   = isub(l_sub+1)
          l_sub = l_sub+1
        level   = isub(l_sub+1)
          l_sub = l_sub+1
        thiad   = isub(l_sub+1)
          l_sub = l_sub+1
!!        TH_FLAG = ISUB(L_SUB+1)
!!          L_SUB = L_SUB+1
!!        ICHAR(TITR)    = ISUB(L_SUB+1)
!!          L_SUB = L_SUB+1
!---
       subset(isu)%ID     = id
       subset(isu)%NCHILD = nchild
       subset(isu)%NPART  = nbpart
       subset(isu)%NTPART = nbtpart
       subset(isu)%PARENT = parent
       subset(isu)%LEVEL  = level
       subset(isu)%THIAD  = thiad
!!       SUBSET(ISU)%TH_FLAG = TH_FLAG
       subset(isu)%TITLE = new_title(isu)
       ALLOCATE(subset(isu)%CHILD(nchild))
       ALLOCATE(subset(isu)%PART(nbpart))
       ALLOCATE(subset(isu)%TPART(nbtpart))
       ALLOCATE(subset(isu)%NVARTH(10))  !  /iTH (9) + /TH (1) -> SUBSET(I)%NVARTH(ITHFLAG)
!---
! NOMBRE DE SUBSETS FILS
        DO j=1,nchild
          ichild = isub(l_sub+1)
            l_sub = l_sub+1
          subset(isu)%CHILD(j) = ichild
        ENDDO
! NOMBRE DE PARTS FILLES
        DO j=1,nbpart
          ipart = isub(l_sub+1)
            l_sub = l_sub+1
          subset(isu)%PART(j) = ipart
        ENDDO
! NOMBRE TOTAL DE PARTS (DESCENDANCE)
        DO j=1,nbtpart
          iparttot = isub(l_sub+1)
            l_sub = l_sub+1
          subset(isu)%TPART(j) = iparttot
        ENDDO
!---
!   - TIME HISTORY - /iTH (9) + /TH (1)
!---
        DO j=1,10
          ivarth =  isub(l_sub+1)
            l_sub = l_sub+1
          subset(isu)%NVARTH(j) = ivarth
        ENDDO
      ENDDO ! DO ISU=1,NSUBS
!---------
      DEALLOCATE (isub)
!---------
      RETURN
      END