fill_gr_surf_ellipse.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/.
!||====================================================================
!||    fill_surf_ellipse   ../starter/source/model/sets/fill_gr_surf_ellipse.F
!||--- called by ------------------------------------------------------
!||    fill_igr            ../starter/source/model/sets/fill_igr.F
!||--- calls      -----------------------------------------------------
!||--- uses       -----------------------------------------------------
!||    message_mod         ../starter/share/message_module/message_mod.F
!||====================================================================
      SUBROUTINE fill_surf_ellipse(SET,IGRSURF,IGRS,BUFSF,LISURF1,NSURF)
C-----------------------------------------------
C   ROUTINE DESCRIPTION :
C   ===================
C   Merge SET%SURFACE into Radioss Surface 
C-----------------------------------------------
C   DUMMY ARGUMENTS DESCRIPTION:
C   ===================
C
C     NAME          DESCRIPTION                         
C
C     SET           Set Structure - Current SET
C     IGRSURF       SURFACES 
C============================================================================
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE my_alloc_mod
      USE message_mod
      USE groupdef_mod
      USE setdef_mod
      USE qa_out_mod
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(INOUT) :: IGRS
      TYPE (SURF_)  , TARGET ,INTENT(INOUT):: IGRSURF(*)
      TYPE (SET_)   , INTENT(INOUT) :: SET
      INTEGER, INTENT(IN) :: LISURF1,NSURF
      my_real, INTENT(INOUT) :: bufsf(lisurf1*(nsurf+nsets))
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NSEG,J,IAD
      CHARACTER MESS*40
      DATA mess/'SET SURF GROUP DEFINITION               '/
C-----------------------------------------------
!
!     create new  (IGRSURF, etc) from elems of /SET
!
      nseg = set%NB_ELLIPSE
!      IF (NSEG == 0) RETURN      ! create a Surface if empty
!---
      igrs = igrs + 1 ! increment NSURF = IGRS + 1
!
      igrsurf(igrs)%ID = set%SET_ID
      igrsurf(igrs)%TITLE = set%TITLE
      igrsurf(igrs)%NSEG = nseg
!
      igrsurf(igrs)%TYPE = 0
      igrsurf(igrs)%ID_MADYMO = 0
      igrsurf(igrs)%IAD_BUFR = 0
      igrsurf(igrs)%NB_MADYMO = 0
      igrsurf(igrs)%TYPE_MADYMO = 0
      igrsurf(igrs)%LEVEL = 1
      igrsurf(igrs)%TH_SURF = 0
      igrsurf(igrs)%ISH4N3N = 0
      igrsurf(igrs)%NSEG_R2R_ALL = 0
      igrsurf(igrs)%NSEG_R2R_SHARE = 0
      igrsurf(igrs)%IAD_IGE = 0
      igrsurf(igrs)%NSEG_IGE = 0
!
!     not printout empty group
!
      IF (nseg == 0) igrsurf(igrs)%SET_GROUP = 1
!
!
      IF (nseg > 0) THEN
!
        CALL my_alloc(igrsurf(igrs)%NODES,nseg,4)
        igrsurf(igrs)%NODES(1:nseg,1:4) = 0
        CALL my_alloc(igrsurf(igrs)%ELTYP,nseg)
        igrsurf(igrs)%ELTYP(1:nseg) = 0
        CALL my_alloc(igrsurf(igrs)%ELEM,nseg)
        igrsurf(igrs)%ELEM(1:nseg) = 0
!
        igrsurf(igrs)%TYPE = 101                                          
        igrsurf(igrs)%IAD_BUFR = set%ELLIPSE_IAD_BUFR
        iad=set%ELLIPSE_IAD_BUFR
        igrsurf(igrs)%ID_MADYMO = set%ELLIPSE_ID_MADYMO
        DO j=1,9
          bufsf(iad+7+j-1)=set%ELLIPSE_SKEW(j)
        ENDDO
 
        bufsf(iad+1)=set%ELLIPSE_A
        bufsf(iad+2)=set%ELLIPSE_B
        bufsf(iad+3)=set%ELLIPSE_C
        bufsf(iad+4)=set%ELLIPSE_XC                                                  
        bufsf(iad+5)=set%ELLIPSE_YC                                                  
        bufsf(iad+6)=set%ELLIPSE_ZC                                                 
        !Init application point for force and momentum     
        !/* ellipsoides : defining center ! */      
        bufsf(iad+16)=set%ELLIPSE_XC                                               
        bufsf(iad+17)=set%ELLIPSE_YC                                                
        bufsf(iad+18)=set%ELLIPSE_ZC 
                                                
        bufsf(iad+36)=set%ELLIPSE_N
 
!
      ENDIF ! IF (NSEG > 0)
 
      set%SET_NSURF_ID = igrs
      set%HAS_SURF_SEG = nseg
C-----
      RETURN
      END