seggetv.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/.
C remplit le tableau segvar avec les valeurs des variables de  l ments voisins
C sert pour les interfaces et pour les ebc
!||====================================================================
!||    seggetv                ../engine/source/interfaces/interf/seggetv.F
!||--- called by ------------------------------------------------------
!||    alemain                ../engine/source/ale/alemain.F
!||--- calls      -----------------------------------------------------
!||    initbuf                ../engine/share/resol/initbuf.F
!||    my_barrier             ../engine/source/system/machine.F
!||    varcondec              ../engine/source/interfaces/interf/varcondec.F
!||--- uses       -----------------------------------------------------
!||    ale_connectivity_mod   ../common_source/modules/ale/ale_connectivity_mod.F
!||    ale_mod                ../common_source/modules/ale/ale_mod.F
!||    elbufdef_mod           ../common_source/modules/mat_elem/elbufdef_mod.F90
!||    initbuf_mod            ../engine/share/resol/initbuf.F
!||    segvar_mod             ../engine/share/modules/segvar_mod.F
!||====================================================================
      SUBROUTINE seggetv(IPARG,ELBUF_TAB,ALE_CONNECTIVITY,ITASK,SEGVAR)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE initbuf_mod
      USE elbufdef_mod    
      USE segvar_mod              
      USE ale_connectivity_mod
      USE ale_mod
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      "com01_c.inc"
#include      "vect01_c.inc"
#include      "param_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IPARG(NPARG,ngroup), ITASK
      TYPE(t_segvar),TARGET :: SEGVAR
      TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
      TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECTIVITY
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NG, I, J, IV, KVAR, KK, JCODV(ALE%GLOBAL%LCONV),CODTOT, SEGAD, IAD2, LGTH
      TYPE(g_bufel_)  ,POINTER :: GBUF     
      TYPE(buf_mat_)  ,POINTER :: MBUF  
C=======================================================================
C
      CALL my_barrier
C
      DO 100 ng=itask+1,ngroup,nthread
         IF (iparg(8,ng) == 1) GO TO 100
c
         CALL varcondec(jcodv,iparg(34,ng),codtot)
         IF (codtot == 0) GOTO 100
         CALL initbuf(iparg    ,ng      ,                     
     2         mtn     ,llt     ,nft     ,iad     ,ity     ,   
     3         npt     ,jale    ,ismstr  ,jeul    ,jtur    ,   
     4         jthe    ,jlag    ,jmult   ,jhbe    ,jivf    ,   
     5         nvaux   ,jpor    ,jcvt    ,jclose  ,jplasol ,   
     6         irep    ,iint    ,igtyp   ,israt   ,isrot   ,   
     7         icsen   ,isorth  ,isorthg ,ifailure,jsms    )            
         lft=1
c         NBB(1) =  NB4
c         NBB(2) =  NB3
c         NBB(3) =  NB10
c         NBB(4) =  NB12
         gbuf => elbuf_tab(ng)%GBUF
         mbuf => elbuf_tab(ng)%BUFLY(1)%MAT(1,1,1)
 
         DO  i=lft,llt
           j=i+nft
           iad2 = ale_connectivity%ee_connect%iad_connect(j)
           lgth = ale_connectivity%ee_connect%iad_connect(j+1)-
     .          ale_connectivity%ee_connect%iad_connect(j)
           DO iv=1,lgth
             IF (ale_connectivity%ee_connect%connected(iad2 + iv - 1) < 0) THEN 
               kk=-ale_connectivity%ee_connect%connected(iad2 + iv - 1)
c
               kvar = 1
               IF (jcodv(kvar) == 1) THEN            
                 segvar%RHO(kk) = gbuf%RHO(i)
               ENDIF                              
c
               kvar = 2
               IF (jcodv(kvar) == 1) THEN           
                 segvar%EINT(kk) = gbuf%EINT(i)
               ENDIF                              
c
               kvar = 3
               IF (jcodv(kvar) == 1) THEN           
                 segvar%RK(kk) = gbuf%RK(i)
               ENDIF                              
c
               kvar = 4
               IF (jcodv(kvar) == 1) THEN            
                 segvar%RE(kk) = gbuf%RE(i)
               ENDIF                              
c
               kvar = 5         
               IF (jcodv(kvar) == 1) THEN    !  UVAR(1)  
                 segvar%UVAR(kk) = mbuf%VAR(llt*(i-1)+1)
               ENDIF                              
 
             ENDIF
           ENDDO
   50    CONTINUE
         ENDDO
  100  CONTINUE
c
c      write(iout,*)'NVCONV =',NVCONV,' NSEGFLU=',NSEGFLU
c      write(iout,'(2(1pE10.3))')(segvar(i),i=1,NVCONV*NSEGFLU)
C
      CALL my_barrier
c-----------
      RETURN
      END