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 ! External functions 39 !----------------------------------------------- 40 DOUBLE PRECISION, EXTERNAL :: CALC_MW 41 !----------------------------------------------- 42 43 IF (MW_AVG /= UNDEFINED) RETURN 44 45 !!$omp parallel do private(ijk) & 46 !!$omp schedule(dynamic,chunk_size) 47 DO IJK = ijkstart3, ijkend3 48 ! calculate mw_mix_g in all fluid and flow boundary cells 49 ! set_bc0 will have already defined mw_mix_g in MI and PI boundary cells 50 ! (redundant-remove in set_bc0?) 51 IF (.NOT.WALL_AT(IJK)) MW_MIX_G(IJK) = & 52 CALC_MW(X_G,DIMENSION_3,IJK,NMAX(0),MW_G) 53 ENDDO 54 55 RETURN 56 END SUBROUTINE SET_MW_MIX_G 57 58 59