alloc_surf_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_surf_str         ../engine/source/groups/alloc_surf_str.f
!||--- called by ------------------------------------------------------
!||    surf_ini               ../engine/source/groups/surf_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
!||    names_and_titles_mod   ../common_source/modules/names_and_titles_mod.F
!||====================================================================
      SUBROUTINE alloc_surf_str(IGRSURF)
C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE message_mod
      USE groupdef_mod
      USE names_and_titles_mod, ONLY: nchartitle
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(surf_)   ,DIMENSION(NSURF) :: IGRSURF
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER ISU,I,J,K,L_SURF,ID,NSEG,TYPE,ID_MADYMO,IAD_BUFR,
     .        NB_MADYMO,TYPE_MADYMO,LEVEL,TH_SURF,ISH4N3N,NSEG_R2R_ALL,
     .        NSEG_R2R_SHARE,NODE,ELTYP,ELEM,ERR,ITITLE(LTITR),LEN
      CHARACTER(LEN=NCHARTITLE) :: TITR
      INTEGER, ALLOCATABLE, DIMENSION (:)  ::  ISURFI
!
!=======================================================================
      DO isu=1,nsurf
        titr(1:nchartitle) = ' '
        CALL read_i_c(ititle,ltitr)
        CALL fretitl2(titr,ititle,ltitr)
        len = 0
        DO j=1,ltitr
          IF (titr(j:j) /= ' ') len = j
        ENDDO
        igrsurf(isu)%TITLE = titr
      ENDDO ! DO ISU=1,NSURF
!!------  Allocate and read subset data from restart
!
      ALLOCATE (isurfi(lenisurf) ,stat=err)
      CALL read_i_c(isurfi, lenisurf)
!=======================================================================
      l_surf = 0
!
      DO isu=1,nsurf
        id             = isurfi(l_surf+1)
          l_surf = l_surf+1
        nseg           = isurfi(l_surf+1)
          l_surf = l_surf+1
        TYPE           = isurfi(l_surf+1)
          l_surf = l_surf+1
        id_madymo      = isurfi(l_surf+1)
          l_surf = l_surf+1
        iad_bufr       = isurfi(l_surf+1)
          l_surf = l_surf+1
        nb_madymo      = isurfi(l_surf+1)
          l_surf = l_surf+1
        type_madymo    = isurfi(l_surf+1)
          l_surf = l_surf+1
        level          = isurfi(l_surf+1)
          l_surf = l_surf+1
        th_surf        = isurfi(l_surf+1)
          l_surf = l_surf+1
        ish4n3n        = isurfi(l_surf+1)
          l_surf = l_surf+1
        nseg_r2r_all   = isurfi(l_surf+1)
          l_surf = l_surf+1
        nseg_r2r_share = isurfi(l_surf+1)
          l_surf = l_surf+1
!
!        ICHAR(TITR)    = ISURFI(L_SURF+1)
!          L_SURF = L_SURF+1
!---
       igrsurf(isu)%ID             = id
       igrsurf(isu)%NSEG           = nseg
       igrsurf(isu)%TYPE           = TYPE
       igrsurf(isu)%ID_MADYMO      = id_madymo
       igrsurf(isu)%IAD_BUFR       = iad_bufr
       igrsurf(isu)%NB_MADYMO      = nb_madymo
       igrsurf(isu)%TYPE_MADYMO    = type_madymo
       igrsurf(isu)%LEVEL          = level
       igrsurf(isu)%TH_SURF        = th_surf
       igrsurf(isu)%ISH4N3N        = ish4n3n
       igrsurf(isu)%NSEG_R2R_ALL   = nseg_r2r_all
       igrsurf(isu)%NSEG_R2R_SHARE = nseg_r2r_share
!
! SURF ENTITIES (NODES, ELTYP, ELEM)
!
       ALLOCATE(igrsurf(isu)%NODES(nseg,4))
       ALLOCATE(igrsurf(isu)%ELTYP(nseg))
       ALLOCATE(igrsurf(isu)%ELEM(nseg))
!
        DO j=1,nseg
          DO k=1,4
            node = isurfi(l_surf+1)
              l_surf = l_surf+1
            igrsurf(isu)%NODES(j,k) = node
          ENDDO
          eltyp = isurfi(l_surf+1)
            l_surf = l_surf+1
          igrsurf(isu)%ELTYP(j) = eltyp
          elem = isurfi(l_surf+1)
            l_surf = l_surf+1
          igrsurf(isu)%ELEM(j) = elem
        ENDDO
      ENDDO ! DO ISU=1,NSURF
!---------
      DEALLOCATE (isurfi)
!---------
      RETURN
      END