schlieren_buffer_gathering.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "param_c.inc"
#include "tabsiz_c.inc"
#include "spmd_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	schlieren_buffer_gathering (nercvois, nesdvois, lercvois, lesdvois, iparg, elbuf_tab, multi_fvm, itherm)

Function/Subroutine Documentation

schlieren_buffer_gathering()

subroutine schlieren_buffer_gathering	(	integer, dimension(snercvois), intent(inout)	nercvois,
		integer, dimension(snesdvois), intent(inout)	nesdvois,
		integer, dimension(slercvois), intent(inout)	lercvois,
		integer, dimension(slesdvois), intent(inout)	lesdvois,
		integer, dimension(nparg,ngroup), intent(in)	iparg,
		type (elbuf_struct_), dimension(ngroup), target	elbuf_tab,
		type(multi_fvm_struct), intent(in)	multi_fvm,
		integer, intent(in)	itherm )

Definition at line 41 of file schlieren_buffer_gathering.F.

C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C     This subroutine outputs is gathering density
C     from all SPMDdomains to calculate later 
C     numerical schlieren from schlieren subroutine
C-----------------------------------------------
C   P r e - C o n d i t i o n s
C-----------------------------------------------
C     ALE/EULER element groups only.
C       this can be checked with IALEL =IPARG(7,NG)+IPARG(11,NG) > 0
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE initbuf_mod
      USE elbufdef_mod  
      USE multi_fvm_mod       
      USE schlieren_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      "com04_c.inc"
#include      "param_c.inc"
#include      "tabsiz_c.inc"
#include      "spmd_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(INOUT) :: NERCVOIS(SNERCVOIS),NESDVOIS(SNESDVOIS),LERCVOIS(SLERCVOIS),LESDVOIS(SLESDVOIS)
      INTEGER,INTENT(IN) :: IPARG(NPARG,NGROUP)
      INTEGER, INTENT(IN) :: ITHERM
      TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
      TYPE(MULTI_FVM_STRUCT), INTENT(IN) :: MULTI_FVM      
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NG,MLW,IALEL,LENCOM,II,NEL,OFFSET,SWA_L
      TYPE(G_BUFEL_)  ,POINTER :: GBUF
      INTEGER :: ITY !< group value : elem type
      INTEGER :: NFT !< group value : shift (to inject local -> global)
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------      
      ii=0
      IF(n2d/=0)ii=1
      swa_l = max(iale,itherm,ieuler,ialelag)*( numels + numelq + ii*numeltg + nsvois+ nqvois + ii*ntgvois + nsegflu )        
      IF(.NOT.ALLOCATED(wa_l))ALLOCATE(wa_l(swa_l)) ! work array for /ANIM/ELEM/SCHLIEREN or /H3D/ELEM/SCHLIEREN (DENSITY GATHERING)            
      DO ng=1,ngroup
        mlw = iparg(01,ng)
        nel = iparg(02,ng)
        ity = iparg(05,ng)
       IF (mlw /= 13) THEN
        DO offset = 0,nel-1,nvsiz                                                                                    
          nft = iparg(3,ng) + offset
          nel = min(nvsiz,nel-offset)
          ialel = iparg(7,ng) + iparg(11,ng)                                                                            
          IF (ialel == 0)cycle
          IF (ity == 1 .OR. ity == 2 .OR. (ity == 7 .AND. n2d /= 0)) THEN
             IF (mlw == 151) THEN                                                                                       
              !-- COLLOCATED SCHEME                                                                                     
              wa_l(nft+1:nft+nel) =  multi_fvm%RHO(nft+1:nft+nel)
             ELSE                                                                                                       
              !-- STAGGERED SCHEME                                                                                      
              gbuf => elbuf_tab(ng)%GBUF                                                                                
              wa_l(nft+1:nft+nel) = gbuf%RHO(1:nel)
             ENDIF                                                                                                      
          ENDIF                                                                                                         
        ENDDO                                                                                                           
       endif!(MLW /= 13)                                                                                                
      enddo!next NG                                                                                               
      !--------------------                                                                                             
      !    SPMD EXCHANGE                                                                                                
      !--------------------    
      IF(iale+ieuler+itherm /= 0)THEN
        lencom = nercvois(nspmd+1)+nesdvois(nspmd+1)                                                                      
        IF (nspmd>1) THEN                                                                                                                                                                                     
          CALL spmd_e1vois( wa_l,nercvois,nesdvois,lercvois,lesdvois,lencom )
        END IF                                                                                                            
      ENDIF