anim_nodal_vector_fvmbags.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "vect01_c.inc"
#include "param_c.inc"
#include "inter22.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	anim_nodal_vector_fvmbags (key, wa4, monvol, volmon, fvdata, nfvbag, smonvol, svolmon, airbags_total_fvm_in_h3d, is_written_node_fvm, airbags_node_id_shift)

Function/Subroutine Documentation

anim_nodal_vector_fvmbags()

subroutine anim_nodal_vector_fvmbags	(	character*4, intent(in)	key,
		dimension(3,airbags_total_fvm_in_h3d), intent(inout)	wa4,
		integer, dimension(smonvol), intent(in)	monvol,
		dimension(svolmon), intent(in)	volmon,
		type(fvbag_data), dimension(nfvbag), intent(in)	fvdata,
		integer, intent(in)	nfvbag,
		integer, intent(in)	smonvol,
		integer, intent(in)	svolmon,
		integer, intent(in)	airbags_total_fvm_in_h3d,
		integer, dimension(airbags_total_fvm_in_h3d), intent(inout)	is_written_node_fvm,
		integer, intent(in)	airbags_node_id_shift )

Definition at line 31 of file anim_nodal_vector_fvmbags.F.

C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C     This suroutine computes nodal vectors from fvmbag polyhedra (fvm case)
C     Done when requested by Engine keyword
C        velocity    : /ANIM/NODA/VEL or /H3D/NODA/VEL
C-----------------------------------------------
C   P r e - C o n d i t i o n s
C-----------------------------------------------
C     none
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE groupdef_mod , only:group_
      USE fvbag_mod , only:fvbag_data !data structure definition
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      "vect01_c.inc"
#include      "param_c.inc"
#include      "inter22.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      CHARACTER*4, INTENT(IN) :: KEY
      INTEGER,INTENT(IN) :: NFVBAG, SMONVOL,SVOLMON,AIRBAGS_TOTAL_FVM_IN_H3D, AIRBAGS_NODE_ID_SHIFT
      my_real,INTENT(INOUT) :: wa4(3,airbags_total_fvm_in_h3d)
      my_real,INTENT(IN) :: volmon(svolmon)
      INTEGER,INTENT(IN) :: MONVOL(SMONVOL)
      TYPE(FVBAG_DATA), INTENT(IN) :: FVDATA(NFVBAG) 
      INTEGER,INTENT(INOUT) :: IS_WRITTEN_NODE_FVM(AIRBAGS_TOTAL_FVM_IN_H3D)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: K1        !< index shift for MONVOL ARRAY
      INTEGER :: NN, II    !< loop
      INTEGER :: ITYP      !< monvol type (refer to read_monvol.F)
      INTEGER :: IH3D_FLAG !< node group identifier
      INTEGER :: NNODES    !< number of nodes in group
      INTEGER :: IFV       !< FVMBAG identifier in [1, NFVBAG] where NFVBAG <= NVOLU
      INTEGER :: nodeID    !< node identifier (internal)
      my_real :: value(3)  !< polyhedron value
      INTEGER :: K_SHIFT   !< shift value for index : sum of previous NPOLH
      INTEGER :: INDX
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------    
 
 
      !loop over all airbags, skip ones which are not FVMBAG (no internal mesh)
      ! then if grnod_id is provided get corresponding value in the output buffer WA4 depending on KEY value (pressure, temperature, density ...)  
      k1=1
      k_shift=0
      DO nn=1,nvolu
         ityp=monvol(k1-1+2)
         IF (ityp == 6.OR.ityp == 8 .OR. ityp == 11) THEN   ! /MONVOL/FVMBAG, or FVMBAG1, or FVMBAG2
            ih3d_flag = monvol(k1-1 +75)  !GRNOD internal identifier in [1:NGRNOD]
            ifv = monvol(k1-1 +45)
            !check if grnod_id was provided, otherwise skip
            IF(ih3d_flag == 1 .AND. ifv /= 0)THEN
              
              IF(fvdata(ifv)%NPOLH > 0)THEN !if GRNOD is not empty and if there are any polyhedrong
                                 
                 SELECT CASE (trim(key))
                 
                   CASE('VEL')
                     DO ii=1,fvdata(ifv)%NPOLH !loop over polyhedra composing the airbag mesh
                       value(1:3) = zero
                       IF(fvdata(ifv)%MPOLH(ii) /= zero)THEN
                         value(1) = fvdata(ifv)%QPOLH(1,ii) / fvdata(ifv)%MPOLH(ii)
                         value(2) = fvdata(ifv)%QPOLH(2,ii) / fvdata(ifv)%MPOLH(ii)
                         value(3) = fvdata(ifv)%QPOLH(3,ii) / fvdata(ifv)%MPOLH(ii)                                                  
                       ENDIF
                       wa4(1:3,k_shift+ii) = value(1:3)
                       is_written_node_fvm(k_shift+ii) = 1
                    enddo! next I (next polyhedron)                                        
 
                END SELECT
 
              ENDIF !NPOLH>0>0 
            ENDIF !IH3D_FLAG>0
            k_shift = k_shift+fvdata(ifv)%NPOLH
         ENDIF ! ITYP
         k1=k1+nimv
 
      ENDDO !next NN
      
      ! set minimum value to other nodes in the group (otherwise visualization may lead to unsuitable color legend/gradient)          
      IF(k_shift < airbags_total_fvm_in_h3d)THEN                                                                       
        DO ii= k_shift+1, airbags_total_fvm_in_h3d
          wa4(1:3, ii) = zero
          is_written_node_fvm(ii) = 0                                                                                                 
        ENDDO                                                                                                                         
      ENDIF        
 
      
      RETURN