simple_rbody_box.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	simple_rbody_box (ibox, x, skew, ib, nd_array, nd_size, rby_msn)

Function/Subroutine Documentation

simple_rbody_box()

subroutine simple_rbody_box	(	type (box_), dimension(nbbox)	ibox,
			x,
			skew,
		integer	ib,
		integer, dimension(*)	nd_array,
		integer	nd_size,
		integer, dimension(2,nrbody), intent(in)	rby_msn )

Definition at line 33 of file simple_rbody_box.F.

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
      INTEGER, INTENT(IN), DIMENSION(2,NRBODY) :: RBY_MSN
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
          DO i=1,nrbody
            j =  rby_msn(2,i)
            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 I=1,NRBODY
       !
       ! CYLIN
        ELSEIF (box_type == 2) THEN
!          DO J=1,NUMNOD
          DO i=1,nrbody
            j =  rby_msn(2,i)
            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 I=1,NRBODY
       !
       ! SPHER
        ELSEIF (box_type == 3) THEN
!          DO J=1,NUMNOD
          DO i=1,nrbody
            j =  rby_msn(2,i)
            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 I=1,NRBODY
        ENDIF ! IF (BOX_TYPE == 1)
!
C---------------
      RETURN