inter_trc_7.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/.
!||====================================================================
!||    inter_trc_7        ../engine/source/interfaces/int07/inter_trc_7.F
!||--- called by ------------------------------------------------------
!||    inttri             ../engine/source/interfaces/intsort/inttri.F
!||--- calls      -----------------------------------------------------
!||    i7trc              ../engine/source/interfaces/intsort/i7trc.F
!||--- uses       -----------------------------------------------------
!||    intbufdef_mod      ../common_source/modules/interfaces/intbufdef_mod.F90
!||    inter_struct_mod   ../engine/share/modules/inter_struct_mod.F
!||====================================================================
        SUBROUTINE inter_trc_7( ITASK,NIN,IPARI,IND_IMP,
     1                          INTBUF_TAB,NB_INTER_SORTED,LIST_INTER_SORTED,INTER_STRUCT)
!$COMMENT
!       INTER_TRC_7 description
!       INTER_TRC_7 organization :      
!$ENDCOMMENT
C============================================================================
C   M o d u l e s
C-----------------------------------------------
        USE intbufdef_mod  
        USE inter_struct_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER, INTENT(in) :: ITASK    ! thread id
        INTEGER, DIMENSION(NPARI,NINTER), INTENT(in) :: IPARI    !   interface data
        INTEGER, DIMENSION(*), INTENT(inout) :: IND_IMP    !   index for implicit option (pointer, cannot define the dimension)
        TYPE(intbuf_struct_),DIMENSION(NINTER) :: INTBUF_TAB    ! interface data
        INTEGER, INTENT(in) :: NB_INTER_SORTED        !   number of interfaces that need to be sorted
        INTEGER, DIMENSION(NB_INTER_SORTED), INTENT(in) :: LIST_INTER_SORTED   !   list of interfaces that need to be sorted
        TYPE(inter_struct_type), DIMENSION(NINTER), INTENT(inout) :: INTER_STRUCT   !   structure for interface
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER :: NIN
        INTEGER :: INACTI,IFQ,NSN,NSNROLD,INACIMP
        INTEGER :: I,N,K,NN,NI,ITIED
        INTEGER :: I_STOK,CAND_T,CAND_TF
        INTEGER :: NUM_IMP
        INTEGER :: KK
C-----------------------------------------------
 
        num_imp = 0
        DO kk=1,nb_inter_sorted
            nin = list_inter_sorted(kk)
            inacti = ipari(22,nin)
            ifq = ipari(31,nin)
            itied = ipari(85,nin)
 
            inter_struct(nin)%INACTI=inacti
            inter_struct(nin)%INACTII = inacti
 
            nsn     =ipari(5,nin)
            IF( inacti==5.OR.inacti==6.OR.inacti==7.OR.ifq>0.OR.
     .          num_imp>0.OR.itied/=0 )THEN
                nsnrold = ipari(24,nin)
            ELSE
                nsnrold = 0
            ENDIF
 
            IF( inacti==5.OR.inacti==6.OR.inacti==7.OR.ifq>0.OR.
     .          num_imp>0.OR.itied/=0)THEN
                inter_struct(nin)%INACTII=inacti
                IF (num_imp>0.AND.(inacti/=5.AND.inacti/=6.AND.ifq<=0)) THEN
                    inacimp = 0
                ELSE
                    inacimp = 1
                ENDIF
    
                inter_struct(nin)%I_SK_OLD = intbuf_tab(nin)%I_STOK(1)
    
                IF( inter_struct(nin)%SIZE_CAND_A<nsn+nsnrold+3) THEN
                    IF(ALLOCATED(inter_struct(nin)%CAND_A)) DEALLOCATE(inter_struct(nin)%CAND_A)
                    ALLOCATE(inter_struct(nin)%CAND_A(nsn+nsnrold+3) )
                    inter_struct(nin)%SIZE_CAND_A=nsn+nsnrold+3
                ENDIF
 
C=============================================
C     LES NOEUDS DELETES DEVIENNENT NSN+NSNROLD+1 
C=============================================
C------for travelling nodes detected in imp_solv----- 
                i_stok = intbuf_tab(nin)%I_STOK(1)
                CALL i7trc(
     1            nsn+nsnrold      ,inter_struct(nin)%I_SK_OLD,intbuf_tab(nin)%CAND_N,intbuf_tab(nin)%CAND_E,
     2            intbuf_tab(nin)%CAND_P,intbuf_tab(nin)%FTSAVX,intbuf_tab(nin)%FTSAVY,intbuf_tab(nin)%FTSAVZ,
     3            inter_struct(nin)%CAND_A,intbuf_tab(nin)%IFPEN ,inacti,ifq,
     4            num_imp          ,ind_imp          ,intbuf_tab(nin)%STFNS ,nin              ,
     5            nsn ,itied,intbuf_tab(nin)%CAND_F)
 
                IF(inter_struct(nin)%I_SK_OLD==0) THEN
                    inter_struct(nin)%INACTI=-abs(inacti)
                ENDIF
                intbuf_tab(nin)%I_STOK(1)=inter_struct(nin)%I_SK_OLD
            ELSE
                inter_struct(nin)%I_SK_OLD=0
                intbuf_tab(nin)%I_STOK(1)=0
                IF(.NOT.ALLOCATED(inter_struct(nin)%CAND_A)) THEN
                    ALLOCATE(inter_struct(nin)%CAND_A(0) )
                    inter_struct(nin)%SIZE_CAND_A=0
                ENDIF
            ENDIF
        ENDDO
 
        RETURN
        END SUBROUTINE inter_trc_7