delnumbc_ply.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/.
!||====================================================================
!||    delnumbc_ply          ../engine/source/output/anim/generate/delnumbc_ply.F
!||--- called by ------------------------------------------------------
!||    genani                ../engine/source/output/anim/generate/genani.F
!||--- calls      -----------------------------------------------------
!||    get_ply_elem          ../engine/source/output/anim/generate/anin_ply_init.F
!||    spmd_iget_partn_ply   ../engine/source/mpi/anim/spmd_iget_partn_ply.F
!||    write_i_c             ../common_source/tools/input_output/write_routines.c
!||--- uses       -----------------------------------------------------
!||    element_mod           ../common_source/modules/elements/element_mod.F90
!||    plyxfem_mod           ../engine/share/modules/plyxfem_mod.F
!||====================================================================
      SUBROUTINE delnumbc_ply(IPLY,NEL_PLY,IPARG,
     2                  IXC ,IXTG ,INVERT,
     3                  EL2FA,NBF  ,INUM ,NELCUT,
     4                  DD_IAD  ,IADD ,NBF_L,
     5                  NBPART,IADG,NODGLOB,IDCMAX,
     6                  NBF_PXFEMG)
      USE plyxfem_mod
      use element_mod , only : nixc,nixtg
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
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "param_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXC(NIXC,*),IXTG(NIXTG,*),IPARG(NPARG,*),
     .         EL2FA(*),INUM(*),
     .        INVERT(*), NBF,NELCUT,NBPART,IADG(NSPMD,*)
      INTEGER DD_IAD(NSPMD+1,*), IADD(*), NBF_L,IPLY,
     .        NEL_PLY,NBF_PXFEMG,NODGLOB(*)
C-----------------------------------------------
C     REAL
      INTEGER II, NG, ITY, LFT, LLT, N, I, 
     . nel, iad, nft,
     . idcmax,ishplyxfem,idply,
     .        ipt,elc,plyelems(numelc),plys,rbuf,nelcount,elem
      INTEGER GET_PLY_ELEM
      EXTERNAL get_ply_elem
C-----------------------------------------------
      nelcount = 0
      DO plys = 1,nplypart
         iply = indx_ply(plys)
 
         plyelems=0
         DO i=1,plyshell(iply)%PLYNUMSHELL
           ipt = plyshell(iply)%SHELLIPT(i)
           elc = plyshell(iply)%SHID(i)
           plyelems(elc)=ipt
         ENDDO
C
        idply  = idcmax + nel_ply
        ii = 0
        DO 490 ng=1,ngroup
          nel   =iparg(2,ng)
          nft   =iparg(3,ng)
          iad   =iparg(4,ng)
          ity   =iparg(5,ng)
          ishplyxfem  = iparg(50,ng)
          lft=1
          llt=nel
C-----------------------------------------------
C       QUAD
C-----------------------------------------------
          IF(ishplyxfem > 0) THEN
C-----------------------------------------------
C       SHELLS
C-----------------------------------------------
            IF(ity==3)THEN
              DO i=lft,llt
                n = i + nft
                IF(plyelems(n) > 0) THEN
                  ii = ii + 1
                    elem = get_ply_elem(iply,n)
                    inum(el2fa(nel_ply + ii)) =
     *                             plyshell(iply)%SHELLID(elem)
                ENDIF
              ENDDO
C-----------------------------------------------
C       3-NODE SHELLS
C-----------------------------------------------
            ELSEIF(ity==7)THEN
              DO i=lft,llt
cc             N = I + NFT
cc             INUM(EL2FA(NN5+N)) = IXTG(NIXTG,N)
              ENDDO
C-----------------------------------------------
C       RNUR
C-----------------------------------------------
            ENDIF
          ENDIF
 490    CONTINUE
C-----------------------------------------------
        nel_ply = nel_ply + plyshell(iply)%PLYNUMSHELL
      ENDDO
 
      IF (nspmd == 1) THEN
        DO i=1,nbf
          inum(i) = inum(i) + idcmax
        ENDDO
        CALL write_i_c(inum,nbf)
 
        idcmax=0
        DO i=1,nbf
           idcmax=max(idcmax,inum(i))
        ENDDO
 
      ELSE
C
        IF (ispmd==0) THEN
          rbuf = nbf_pxfemg
          CALL spmd_iget_partn_ply(1,nbf_l,inum,nplypart,iadg,rbuf,1,
     *                             idcmax)
        ELSE
          rbuf = 1
          CALL spmd_iget_partn_ply(1,nbf_l,inum,nplypart,iadg,rbuf,1,
     *                             idcmax)
        END IF
      ENDIF
 
C-----------------------------------------------
C
      RETURN
      END