simple_node_box.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/.
!||====================================================================
!||    simple_node_box   ../starter/source/model/sets/simple_node_box.F
!||--- called by ------------------------------------------------------
!||    node_box          ../starter/source/model/sets/fill_clause_node_box.f
!||--- calls      -----------------------------------------------------
!||    checkcyl          ../starter/source/model/box/rdbox.F
!||    checkpara         ../starter/source/model/box/rdbox.F
!||    checksphere       ../starter/source/model/box/rdbox.F
!||--- uses       -----------------------------------------------------
!||====================================================================
      SUBROUTINE simple_node_box(
     .                    IBOX   ,X    ,SKEW  ,
     .                    IB, ND_ARRAY,ND_SIZE)
C-----------------------------------------------
C   ROUTINE DESCRIPTION :
C   ===================
C   create node list from BOX
C------------------------------------------------------------------
C   DUMMY ARGUMENTS DESCRIPTION:
C   ===================
C
C     NAME          DESCRIPTION                         
C
C     IBOX          BOX structure
C     X             Node position 
C     SKEW          SKEW array
C     NOD_ARRAY     Result list of nodes
C     ND_SIZE       number of stacked nodes
C     IB            Box to treat
C============================================================================
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE optiondef_mod
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 ND_ARRAY(*),IB,ND_SIZE
      my_real
     .      x(3,*),skew(lskew,*)
C-----------------------------------------------
      TYPE (BOX_) , DIMENSION(NBBOX)  :: IBOX
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,INSIDE,ISK,BOX_TYPE,NBOXBOX,IBX
      my_real
     .    xp1,yp1,zp1,xp2,yp2,zp2,diam,nodinb(3)
C-----------------------------------------------
       nd_size = 0
C------------
!      IBOX(IGS)%ISKBOX = ISK
!      IBOX(IGS)%NOD1 = J2(1)
!      IBOX(IGS)%NOD2 = J2(2)
C      IBOX(IGS)%NBOXBOX  : NUMBER OF SUB BOXES (BOXES OF BOXES)
!      IF(KEY(1:5) == 'RECTA')IBOX(IGS)%TYPE = 1
!      IF(KEY(1:5) == 'CYLIN')IBOX(IGS)%TYPE = 2
!      IF(KEY(1:5) == 'SPHER')IBOX(IGS)%TYPE = 3
!      IBOX(IGS)%DIAM = DIAM
!      IBOX(IGS)%X1 = XP1
!      IBOX(IGS)%Y1 = YP1
!      IBOX(IGS)%Z1 = ZP1
!      IBOX(IGS)%X2 = XP2
!      IBOX(IGS)%Y2 = YP2
!      IBOX(IGS)%Z2 = ZP2
C------------
        ibx = abs(ib) ! a box can be have negative user_ID if within Box of Box
        isk = ibox(ibx)%ISKBOX
        box_type = ibox(ibx)%TYPE
        xp1 = ibox(ibx)%X1
        yp1 = ibox(ibx)%Y1
        zp1 = ibox(ibx)%Z1
        xp2 = ibox(ibx)%X2
        yp2 = ibox(ibx)%Y2
        zp2 = ibox(ibx)%Z2
        diam = ibox(ibx)%DIAM
       !
       ! RECTA
        IF (box_type == 1) THEN
          DO j=1,numnod
            inside = 0
            nodinb(1) = x(1,j)
            nodinb(2) = x(2,j)
            nodinb(3) = x(3,j)
            CALL checkpara(xp1,yp1,zp1,xp2,yp2,zp2,
     .                     isk,nodinb,skew,inside)
 
            IF (inside == 1) THEN 
              nd_size = nd_size + 1
              nd_array(nd_size) = j ! add node
            ENDIF
 
          ENDDO ! DO J=1,NUMNOD
       !
       ! CYLIN
        ELSEIF (box_type == 2) THEN
          DO j=1,numnod
            inside = 0
            nodinb(1) = x(1,j)
            nodinb(2) = x(2,j)
            nodinb(3) = x(3,j)
            CALL checkcyl(xp1, yp1, zp1 , xp2, yp2, zp2,
     .                    nodinb  , diam, inside )
 
            IF (inside == 1) THEN 
              nd_size = nd_size + 1
              nd_array(nd_size) = j ! add node
            ENDIF
 
          ENDDO ! DO J=1,NUMNOD
       !
       ! SPHER
        ELSEIF (box_type == 3) THEN
          DO j=1,numnod
            inside = 0
            nodinb(1) = x(1,j)
            nodinb(2) = x(2,j)
            nodinb(3) = x(3,j)
            CALL checksphere(xp1, yp1, zp1, nodinb, diam, inside)
 
            IF (inside == 1) THEN 
              nd_size = nd_size + 1
              nd_array(nd_size) = j ! add node
            ENDIF
 
          ENDDO ! DO J=1,NUMNOD
        ENDIF ! IF (BOX_TYPE == 1)
!
C---------------
      RETURN
      END