w_igrav.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_igrav     ../starter/source/restart/ddsplit/w_igrav.F
!||--- called by ------------------------------------------------------
!||    ddsplit     ../starter/source/restart/ddsplit/ddsplit.F
!||--- calls      -----------------------------------------------------
!||====================================================================
      SUBROUTINE w_igrav(
     .  IGRAV,LGRAV,NODLOCAL,PROC,LLGRAV_L,LEN_IA,NUMNOD_L)
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"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER PROC, LEN_IA, LLGRAV_L,NUMNOD_L,
     .        igrav(nigrv,*), lgrav(*), nodlocal(*)
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 I, J, IAD, IADL, NL, N, NLL, NN, IGRTMP(NIGRV,NGRAV)
      INTEGER, DIMENSION(:), ALLOCATABLE :: LGRTMP ! LLGRAV_L
C
      ALLOCATE(lgrtmp(llgrav_l))
      iadl = 1
      DO i = 1, ngrav
        nl = igrav(1,i)
        iad= igrav(4,i)
        DO j = 1, nigrv
          igrtmp(j,i) = igrav(j,i)
        ENDDO
        nll = 0
        DO j = 1, nl
          nn = lgrav(iad+j-1)
          n = abs(nn)
          IF (nodlocal(n)/=0.AND.nodlocal(n)<=numnod_l)THEN
            IF(nn>0) THEN
              lgrtmp(iadl+nll) = nodlocal(n)
            ELSE
              lgrtmp(iadl+nll) = -nodlocal(n)
            ENDIF
            nll = nll + 1
          ENDIF
        ENDDO
        igrtmp(1,i) = nll
        igrtmp(4,i) = iadl
        iadl = iadl+nll
      ENDDO
C
      CALL write_i_c(igrtmp,nigrv*ngrav)
      CALL write_i_c(lgrtmp,iadl-1)
      len_ia = len_ia + nigrv*ngrav + iadl-1
 
      DEALLOCATE(lgrtmp)
C
      RETURN
      END