bhol2.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/.
!||====================================================================
!||    bhol2         ../engine/source/ale/ale2d/bhol2.F
!||--- called by ------------------------------------------------------
!||    alethe        ../engine/source/ale/alethe.F
!||--- uses       -----------------------------------------------------
!||    element_mod   ../common_source/modules/elements/element_mod.F90
!||====================================================================
      SUBROUTINE bhol2(ITASK ,IPARG ,PM ,BHOLE ,IXQ   )
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------        
      use element_mod , only : nixq        
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      "com01_c.inc"
#include      "com04_c.inc"
#include      "vect01_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 ITASK,IPARG(NPARG,NGROUP),BHOLE(NMULT,*),IXQ(NIXQ,NUMELQ)
      my_real pm(npropm,nummat)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NM, NG, JMUL, I, J, MX
C-----------------------------------------------
C   S o u r c e   L i n e s
C------------------------------------------------
      DO nm=1,nmult
       DO ng=itask+1,ngroup,nthread
C       ALE ON / OFF
        IF (iparg(76, ng)  ==  1) cycle ! --> OFF
        jale=iparg(7,ng)
        jeul=iparg(11,ng)
        IF(jale+jeul == 0)   cycle
        IF(iparg(8,ng) == 1) cycle
        jmul=iparg(20,ng)
        llt=iparg(2,ng)
        nft=iparg(3,ng)
        lft=1
        IF(jmul == 0)THEN
          DO i=lft,llt
           j=i+nft
           mx=ixq(1,j)
           bhole(nm,j)=mx
          ENDDO
C       ELSEIF(JMUL<NM)THEN
C         DO I=LFT,LLT
C          J=I+NFT
C          MX=IXQ(1,J)
C          BHOLE(NM,J)=0
C         ENDDO
        ELSE
          DO i=lft,llt
           j=i+nft
           mx=ixq(1,j)
           bhole(nm,j)=nint(pm( 20+nm,mx))
           IF(nint(pm(185+nm,mx)) /= 0)bhole(nm,j)=-bhole(nm,j)
          ENDDO
        ENDIF
       enddo! next NG
      ENDDO !next NM=1,NMULT
C-----------------------------------------------      
      RETURN
      END