1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Subroutine: SET_MW_MIX_g C 4 ! Purpose: calculate gas mixture molecular weights C 5 ! C 6 ! Author: M. Syamlal Date: 19-OCT-92 C 7 ! Reviewer: S. Venkatesan Date: 11-DEC-92 C 8 ! C 9 ! Literature/Document References: C 10 ! C 11 ! Variables referenced: IJKMAX2, X_g C 12 ! Variables modified: MW_MIX_g C 13 ! Local variables: NONE C 14 ! C 15 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 16 17 SUBROUTINE SET_MW_MIX_G 18 19 !----------------------------------------------- 20 ! Modules 21 !----------------------------------------------- 22 USE param 23 USE param1 24 USE parallel 25 USE physprop 26 USE geometry 27 USE fldvar 28 USE constant 29 USE indices 30 USE compar 31 USE functions 32 IMPLICIT NONE 33 !----------------------------------------------- 34 ! Local variables 35 !----------------------------------------------- 36 INTEGER :: IJK 37 !----------------------------------------------- 38 39 IF (MW_AVG /= UNDEFINED) RETURN 40 41 !!$omp parallel do private(ijk) & 42 !!$omp schedule(dynamic,chunk_size) 43 DO IJK = ijkstart3, ijkend3 44 ! calculate mw_mix_g in all fluid and flow boundary cells 45 ! set_bc0 will have already defined mw_mix_g in MI and PI boundary cells 46 ! (redundant-remove in set_bc0?) 47 IF (.NOT.WALL_AT(IJK)) MW_MIX_G(IJK) = & 48 CALC_MW(X_G,DIMENSION_3,IJK,NMAX(0),MW_G) 49 ENDDO 50 51 RETURN 52 END SUBROUTINE SET_MW_MIX_G 53 54 55