read_pcyl.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_pcyl   ../engine/source/output/restart/read_pcyl.F
!||--- called by ------------------------------------------------------
!||    rdresb      ../engine/source/output/restart/rdresb.F
!||--- calls      -----------------------------------------------------
!||    read_db     ../common_source/tools/input_output/read_db.F
!||    read_i_c    ../common_source/tools/input_output/write_routtines.c
!||--- uses       -----------------------------------------------------
!||    loads_mod   ../common_source/modules/loads/loads_mod.F90
!||====================================================================
      SUBROUTINE read_pcyl(LOADS)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE loads_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      "units_c.inc"
#include      "tabsiz_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE (LOADS_)  :: LOADS
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I,J,NSEG,LLINT,IAD,NLOAD_CYL
      INTEGER :: MY_SIZE
      INTEGER ,DIMENSION(:) ,ALLOCATABLE :: ITMP
      my_real, DIMENSION(3) :: rtmp
C-----------------------------------------------
c         -> PCYL_ID
c         -> NSEG
c         -> SEGNOD(NSEG,4)  (N1,N2,NB3,N4) by segment
c         -> FRAME_ID
c         -> SENS_ID
c         -> TABLE_ID
c         -> XSCALE_R
c         -> XSCALE_T
c         -> YSCALE_P
c=======================================================================
      nload_cyl = loads%NLOAD_CYL
      ALLOCATE(loads%LOAD_CYL(nload_cyl))
c      
      DO i=1,nload_cyl
        CALL read_i_c(nseg,1)
        llint = 4 + nseg*4
        ALLOCATE (itmp(llint))
        CALL read_i_c(itmp,llint)
        
        loads%LOAD_CYL(i)%ID       = itmp(1)
        loads%LOAD_CYL(i)%ISENS    = itmp(2)
        loads%LOAD_CYL(i)%ITABLE   = itmp(3)
        loads%LOAD_CYL(i)%IFRAME   = itmp(4)
c
        ALLOCATE(loads%LOAD_CYL(i)%SEGNOD(nseg,4))
        iad = 4
        loads%LOAD_CYL(i)%NSEG = nseg
        DO j=1,nseg
          loads%LOAD_CYL(i)%SEGNOD(j,1) = itmp(iad+1)
          loads%LOAD_CYL(i)%SEGNOD(j,2) = itmp(iad+2)
          loads%LOAD_CYL(i)%SEGNOD(j,3) = itmp(iad+3)
          loads%LOAD_CYL(i)%SEGNOD(j,4) = itmp(iad+4)
          iad = iad + 4
        END DO
c
        CALL read_db(rtmp,3)
        loads%LOAD_CYL(i)%XSCALE_R = rtmp(1) 
        loads%LOAD_CYL(i)%XSCALE_T = rtmp(2) 
        loads%LOAD_CYL(i)%YSCALE   = rtmp(3) 
 
        ! ---------------
        ! load/pcyl option : adress of FSKY array
        my_size = loads%LOAD_CYL(i)%NSEG
        ALLOCATE(loads%LOAD_CYL(i)%SEGMENT_ADRESS(4,my_size))
        CALL read_i_c(loads%LOAD_CYL(i)%SEGMENT_ADRESS,4*my_size)
        ! ---------------
c        
        DEALLOCATE(itmp)
      ENDDO
c-----------
      RETURN
      END