w_elmx.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_elmx      ../starter/source/restart/ddsplit/w_elmx.F
!||--- called by ------------------------------------------------------
!||    ddsplit     ../starter/source/restart/ddsplit/ddsplit.F
!||--- calls      -----------------------------------------------------
!||====================================================================
      SUBROUTINE w_elmx(KXX,NUMELX_L,ISUMNX_L,IXX,CEP,
     +                  PROC,NODLOCAL,LEN_IA)
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      "scr23_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NUMELX_L, PROC, ISUMNX_L,LEN_IA,
     .        kxx(nixx,*), cep(*), nodlocal(*), ixx(*)
C-----------------------------------------------
C   L o c a l  V a r i a b l e s
C-----------------------------------------------
      INTEGER I, NN, IAD, IAD_L, J, NB_L, NNOD, IE_L, OFF
      INTEGER, DIMENSION(:), ALLOCATABLE :: IXX_L
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: KXX_L
C-----------------------------------------------
C
! ----------------------------
!     allocate 1d array
      ALLOCATE( ixx_l(isumnx_l) ) 
!     allocate 2d array
      ALLOCATE( kxx_l(nixx,numelx_l) )
! ----------------------------
      off = numelc+numeltg+numels+numelq+numelt+numelp+numelr
      ie_l = 0
      nb_l = 0
C
      DO i = 1, numelx
        IF(cep(off+i)==proc) THEN
          nnod  = kxx(3,i)
          iad   = kxx(4,i)
C
          ie_l = ie_l + 1
          DO j = 1, nixx
            kxx_l(j,ie_l) = kxx(j,i)
          ENDDO
          iad_l = nb_l+1
          kxx_l(4,ie_l) = iad_l
          nb_l = nb_l + nnod
          DO j = 1, nnod
            ixx_l(iad_l+j-1) = nodlocal(ixx(iad+j-1))
          ENDDO
        ENDIF
      ENDDO
C
      CALL write_i_c(kxx_l,numelx_l*nixx)
      CALL write_i_c(ixx_l,isumnx_l)
      len_ia = len_ia + numelx_l*nixx + isumnx_l
C
! ----------------------------
!     deallocate 1d array
      DEALLOCATE( ixx_l ) 
!     deallocate 2d array
      DEALLOCATE( kxx_l )
! ----------------------------
      RETURN
      END