write_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/.
!||====================================================================
!||    write_pcyl      ../starter/source/loads/general/load_pcyl/write_pcyl.F
!||--- called by ------------------------------------------------------
!||    ddsplit         ../starter/source/restart/ddsplit/ddsplit.F
!||--- calls      -----------------------------------------------------
!||--- uses       -----------------------------------------------------
!||    message_mod     ../starter/share/message_module/message_mod.F
!||====================================================================
      SUBROUTINE write_pcyl(LOAD_CYL ,NLOAD_CYL ,NODLOCAL)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE pload_cyl_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER :: NLOAD_CYL
      TYPE (PRESS_CYL_) ,DIMENSION(NLOAD_CYL) ,INTENT(IN) :: LOAD_CYL
      INTEGER , DIMENSION(*) ,INTENT(IN) :: NODLOCAL
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I,J,N1,N2,N3,N4,NSEG,LINT,IAD,OFF
      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=======================================================================
 
      ! --------------
      ! loop over the /load/pcyl
      DO i=1,nload_cyl        
        nseg = load_cyl(i)%NSEG ! get the number of segment 
        lint = 5 + nseg*4
        ALLOCATE(itmp(lint))
        rtmp(1) = load_cyl(i)%XSCALE_R
        rtmp(2) = load_cyl(i)%XSCALE_T
        rtmp(3) = load_cyl(i)%YSCALE
c
        itmp(1) = nseg
        itmp(2) = load_cyl(i)%ID
        itmp(3) = load_cyl(i)%ISENS
        itmp(4) = load_cyl(i)%ITABLE
        itmp(5) = load_cyl(i)%IFRAME
        iad     = 5
        ! --------------
        ! loop over the segment to convert global node id into local node id
        DO j=1,load_cyl(i)%NSEG
            n1 = load_cyl(i)%SEGNOD(j,1) ! get the global node id N1
            n2 = load_cyl(i)%SEGNOD(j,2) ! get the global node id N2
            n3 = load_cyl(i)%SEGNOD(j,3) ! get the global node id N3
            n4 = load_cyl(i)%SEGNOD(j,4) ! get the global node id N4
            itmp(iad+1) = nodlocal(n1) ! global --> local for N1
            itmp(iad+2) = nodlocal(n2) ! global --> local for N2
            itmp(iad+3) = nodlocal(n3) ! global --> local for N3
            ! global --> local for N4 : for triangle, N4 = 0, for shell N4/=0
            IF (n4 == 0) THEN
                itmp(iad+4) = 0
            ELSE
                itmp(iad+4) = nodlocal(n4)
            ENDIF
            iad = iad + 4
        ENDDO
        ! --------------
 
        ! --------------        
        ! write the data related to the load/pcyl
        CALL write_i_c(itmp,iad)
        CALL write_db(rtmp,3)
        ! /LOAD/PCYL option : adress for parith/on 
        CALL write_i_c(load_cyl(i)%SEGMENT_ADRESS,4*load_cyl(i)%S_SEGMENT_ADRESS)
        ! --------------
                
        DEALLOCATE(itmp)
      ENDDO
      ! --------------
 
 
      RETURN
      END