surface_min_max_computation.F File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	surface_min_max_computation (nsurf, numnod, surface_number, min_max_position, x, igrsurf, inivol)

Function/Subroutine Documentation

surface_min_max_computation()

subroutine surface_min_max_computation	(	integer, intent(in)	nsurf,
		integer, intent(in)	numnod,
		integer, intent(in)	surface_number,
		intent(inout)	min_max_position,
		intent(in)	x,
		type (surf_), dimension(nsurf), intent(in)	igrsurf,
		type (inivol_struct_), intent(in)	inivol )

Definition at line 30 of file surface_min_max_computation.F.

 
        USE inivol_def_mod , only : inivol_struct_, container_
        USE groupdef_mod , only : surf_
!$COMMENT
!       SURFACE_MIN_MAX_COMPUTATION description
!       SURFACE_MIN_MAX_COMPUTATION computes the min/max position of surface's nodes
!       
!       SURFACE_MIN_MAX_COMPUTATION organization :
!$ENDCOMMENT
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER,INTENT(IN) :: NSURF,NUMNOD
        INTEGER, INTENT(IN) :: SURFACE_NUMBER ! number of surface
        my_real, DIMENSION(6), INTENT(INOUT) :: min_max_position ! min/max position
        my_real, DIMENSION(3,NUMNOD), INTENT(IN) :: x ! position
        TYPE (SURF_), DIMENSION(NSURF), INTENT(IN) :: IGRSURF ! surface data
        TYPE (INIVOL_STRUCT_), INTENT(IN) :: INIVOL ! inivol data
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER :: I,J,K
        INTEGER :: SURFACE_ID,SURFACE_NODE_NUMBER,SURFACE_TYPE
        INTEGER :: NODE_ID
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
        DO i=1,surface_number
          surface_id = inivol%CONTAINER(i)%SURF_ID ! surface id
          surface_node_number = igrsurf(surface_id)%NSEG ! number of segment of the surface
          surface_type = igrsurf(surface_id)%TYPE ! type of surface
          IF(surface_type /= 200 .AND. surface_type /= 101) THEN
            DO j=1,4
              DO k=1,surface_node_number
                node_id = igrsurf(surface_id)%NODES(k,j)
                min_max_position(1:3) = min(min_max_position(1:3),x(1:3,node_id))
                min_max_position(4:6) = max(min_max_position(4:6),x(1:3,node_id))
              ENDDO
            ENDDO
          ENDIF
        ENDDO
        ! --------------------------
 
        RETURN