set_mw_mix_g.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Subroutine: SET_MW_MIX_g                                            C
!  Purpose: calculate gas mixture molecular weights                    C
!                                                                      C
!  Author: M. Syamlal                                 Date: 19-OCT-92  C
!  Reviewer: S. Venkatesan                            Date: 11-DEC-92  C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced: IJKMAX2, X_g                                  C
!  Variables modified: MW_MIX_g                                        C
!  Local variables: NONE                                               C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

      SUBROUTINE set_mw_mix_g

!-----------------------------------------------
! Modules
!-----------------------------------------------
      USE param
      USE param1
      USE parallel
      USE physprop
      USE geometry
      USE fldvar
      USE constant
      USE indices
      USE compar
      USE functions
      IMPLICIT NONE
!-----------------------------------------------
! Local variables
!-----------------------------------------------
      INTEGER :: IJK
!-----------------------------------------------

      IF (mw_avg /= undefined) RETURN

!!$omp parallel do private(ijk) &
!!$omp schedule(dynamic,chunk_size)
      DO ijk = ijkstart3, ijkend3
! calculate mw_mix_g in all fluid and flow boundary cells
! set_bc0 will have already defined mw_mix_g in MI and PI boundary cells
! (redundant-remove in set_bc0?)
         IF (.NOT.wall_at(ijk)) mw_mix_g(ijk) = &
            calc_mw(x_g,dimension_3,ijk,nmax(0),mw_g)
      ENDDO

      RETURN
      END SUBROUTINE set_mw_mix_g