w_secbuf.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/.
C
!||====================================================================
!||    w_secbuf   ../starter/source/restart/ddsplit/w_secbuf.F
!||--- called by ------------------------------------------------------
!||    ddsplit    ../starter/source/restart/ddsplit/ddsplit.F
!||--- calls      -----------------------------------------------------
!||    nlocal     ../starter/source/spmd/node/ddtools.F
!||====================================================================
      SUBROUTINE w_secbuf(SECBUF ,LSECBUF_L, NSTRF, PROC,
     2                    LEN_AM )
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------  
 
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  PROC, LSECBUF_L, LEN_AM,
     .          nstrf(*)
      my_real
     .          secbuf(*)
C-----------------------------------------------
C   F u n c t i o n
C-----------------------------------------------
      INTEGER  NLOCAL
      EXTERNAL nlocal    
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NNOD_S, NNOD_S_L, TYP_S, NSINT_S,
     .        n, ip, j, k0, k2, kr_l, kr, k
      my_real
     .        secbuf_l(lsecbuf_l)
C
      kr_l = 0
      kr = 0
C header commun
      DO j = 1, 10
        secbuf_l(kr_l+j) = secbuf(kr+j)
      END DO
      kr_l = kr_l + 10
      kr = kr + 10
C
      k0 = nstrf(25)
      DO n = 1, nsect
        typ_s  = nstrf(k0)
        nnod_s = nstrf(k0+6)
        nsint_s = nstrf(k0+14)
        k2=k0+30+nsint_s
C noeuds locaux
        nnod_s_l = 0
        DO j = 1, nnod_s
          k = nstrf(k2 + j -1)
          IF(nlocal(k,proc+1)==1)THEN
            nnod_s_l = nnod_s_l + 1
          END IF
        END DO
        DO j = 1, 10
          secbuf_l(kr_l+j) = secbuf(kr+j)
        END DO
        kr_l = kr_l + 10
        kr = kr + 10
        IF(typ_s>=100) THEN
          DO j = 1, 12*nnod_s_l
            secbuf_l(kr_l+j) = zero
          END DO
          kr_l = kr_l + 12*nnod_s_l
          kr = kr + 12*nnod_s
        END IF
        IF(typ_s>=101) THEN
          DO j = 1, 12*nnod_s_l
            secbuf_l(kr_l+j) = zero
          END DO
          kr_l = kr_l + 12*nnod_s_l
          kr = kr + 12*nnod_s
        END IF
        IF(typ_s>=102) THEN
          DO j = 1, 6*nnod_s_l
            secbuf_l(kr_l+j) = zero
          END DO
          kr_l = kr_l + 6*nnod_s_l
          kr = kr + 6*nnod_s
        END IF
        k0  = nstrf(k0+24)
      ENDDO
C
      CALL write_db(secbuf_l,lsecbuf_l)
      len_am = len_am + lsecbuf_l
      if(lsecbuf_l/=kr_l)print*,'error secbuf',lsecbuf_l,kr_l,nnod_s_l
C
      RETURN
      END