surface_min_max_computation.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/.
!||====================================================================
!||    surface_min_max_computation   ../starter/source/initial_conditions/inivol/surface_min_max_computation.F
!||--- called by ------------------------------------------------------
!||    init_inivol                   ../starter/source/initial_conditions/inivol/init_inivol.F90
!||--- uses       -----------------------------------------------------
!||    inivol_def_mod                ../starter/share/modules1/inivol_mod.F
!||====================================================================
      SUBROUTINE surface_min_max_computation(
     .                                       NSURF, NUMNOD,
     .                                       SURFACE_NUMBER,MIN_MAX_POSITION,X,IGRSURF,INIVOL)
 
        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
        END SUBROUTINE surface_min_max_computation