insurfigeo.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/.
!||====================================================================
!||    insurfigeo             ../starter/source/interfaces/interf1/insurfigeo.F
!||--- called by ------------------------------------------------------
!||    lecins                 ../starter/source/interfaces/interf1/lecins.F
!||    lecint                 ../starter/source/interfaces/interf1/lecint.F
!||--- calls      -----------------------------------------------------
!||    countigeo              ../starter/source/interfaces/interf1/countigeo.F
!||--- uses       -----------------------------------------------------
!||    format_mod             ../starter/share/modules1/format_mod.F90
!||    message_mod            ../starter/share/message_module/message_mod.F
!||====================================================================
      SUBROUTINE insurfigeo(NRT_IGE ,OFFSET_SEG   ,MSN_IGE,OFFSET_NODE,IADTABIGE,IR    ,IRECT  ,NOINT  ,
     .                 SURF_NODES_IGE,ITAB  ,MSV   ,ID,TITR,
     .                 NIGE,RIGE,XIGE,VIGE,NIGE_TMP,RIGE_TMP,XIGE_TMP,VIGE_TMP)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE names_and_titles_mod , ONLY : nchartitle
      USE format_mod , ONLY : fmw_4i
C-----------------------------------------------------------------------
C     READING SURFACES AND COUNTING SEGMENTS
C           ENTREE : 
C                    NRT_IGE NUMBER OF ITEMS TO READ
C           SORTIE : 
C                    IRECT 
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      "units_c.inc"
#include      "scr03_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NRT_IGE,OFFSET_SEG,NOINT,IR,MSN_IGE,OFFSET_NODE,IADTABIGE
      INTEGER IRECT(4,*),
     .        itab(*), msv(*),nige_tmp(*),nige(*),surf_nodes_ige(nrt_ige,4)
      INTEGER ID
      CHARACTER(LEN=NCHARTITLE) :: TITR
      my_real
     .   rige(3,*),xige(3,*),vige(3,*),
     .   rige_tmp(3,*),xige_tmp(3,*),vige_tmp(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, N23, J, K, IR1, K4
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
C=======================================================================   
CC in intbuf_tab%nige, we write all the fictitious points second then the fictitious points main
c NIGE is offset at the routine call
 
      DO j=1,msn_ige
        nige(j)=nige_tmp(j+iadtabige) 
        rige(:,j)=rige_tmp(:,j+iadtabige)       
        xige(:,j)=xige_tmp(:,j+iadtabige)                                     
        vige(:,j)=vige_tmp(:,j+iadtabige)  
        xige(:,j)=xige_tmp(:,j+iadtabige) 
      ENDDO
      DO j=offset_seg+1,offset_seg+nrt_ige                                            
        DO k=1,4                                          
          irect(k,j)=surf_nodes_ige(j,k)+offset_node
        ENDDO    
      ENDDO                                                 
C-------    
      IF(ir/=0.AND.ir/=2)THEN              
        DO j=offset_seg+1,offset_seg+nrt_ige                 
          ir1=irect(1,j)           
          irect(1,j)=irect(2,j)    
          irect(2,j)=ir1           
          ir1=irect(3,j)           
          irect(3,j)=irect(4,j)    
          irect(4,j)=ir1            
        ENDDO                      
      ENDIF                        
      n23=4                        
C-------          
      IF(ipri>=1) THEN                                              
        WRITE(iout,'(/,A,/)')' SEGMENTS USED FOR ISOGEOMETRIC SURFACE DEFINITION'  
        k=1
        k4=4
        DO i=offset_seg+1,offset_seg+nrt_ige  
          WRITE(iout,fmt=fmw_4i)(irect(k,i),k=1,k4)  !  ECRITURE NON DEFINITIVE
        ENDDO                                                         
      ENDIF                                                           
      IF (nrt_ige/=0) CALL countigeo(irect,msv,msn_ige,nrt_ige,offset_seg)                    
C------------------------------------------------------------
      RETURN
      END