fill_gr_surf_ellipse.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	fill_surf_ellipse (set, igrsurf, igrs, bufsf, lisurf1, nsurf)

Function/Subroutine Documentation

fill_surf_ellipse()

subroutine fill_surf_ellipse	(	type (set_), intent(inout)	set,
		type (surf_), dimension(*), intent(inout), target	igrsurf,
		integer, intent(inout)	igrs,
		dimension(lisurf1*(nsurf+nsets)), intent(inout)	bufsf,
		integer, intent(in)	lisurf1,
		integer, intent(in)	nsurf )

Definition at line 31 of file fill_gr_surf_ellipse.F.

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