1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Module name: VAVG_V_s(M) C 4 ! Purpose: Volume average V_s C 5 ! C 6 ! Author: M. Syamlal Date: 28-APR-94 C 7 ! Reviewer: Date: dd-mmm-yy C 8 ! C 9 ! Revision Number: C 10 ! Purpose: C 11 ! Author: Date: dd-mmm-yy C 12 ! Reviewer: Date: dd-mmm-yy C 13 ! C 14 ! Literature/Document References: C 15 ! C 16 ! Variables referenced: C 17 ! Variables modified: C 18 ! C 19 ! Local variables: C 20 ! C 21 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 22 ! 23 DOUBLE PRECISION FUNCTION VAVG_V_S (M) 24 !...Translated by Pacific-Sierra Research VAST-90 2.06G5 12:17:31 12/09/98 25 !...Switches: -xf 26 ! 27 !----------------------------------------------- 28 ! M o d u l e s 29 !----------------------------------------------- 30 USE param 31 USE param1 32 USE parallel 33 USE fldvar 34 USE bc 35 USE geometry 36 USE physprop 37 USE indices 38 USE compar 39 USE mpi_utility 40 USE functions 41 IMPLICIT NONE 42 !----------------------------------------------- 43 ! G l o b a l P a r a m e t e r s 44 !----------------------------------------------- 45 !----------------------------------------------- 46 ! D u m m y A r g u m e n t s 47 !----------------------------------------------- 48 ! 49 ! Indices 50 INTEGER IJK , M 51 ! 52 ! Integral of V_s*EP_s for entire volume 53 DOUBLE PRECISION SUM_V_s 54 ! 55 ! Total volume of computational cells 56 DOUBLE PRECISION SUM_VOL 57 58 ! Integrate the velocity values for the whole domain 59 ! 60 SUM_V_S = ZERO 61 SUM_VOL = ZERO 62 63 !!!!$omp parallel do private(IJK) reduction(+:SUM_VOL,SUM_V_S) 64 DO IJK = IJKSTART3, IJKEND3 65 IF(.NOT.IS_ON_myPE_wobnd(I_OF(IJK), J_OF(IJK), K_OF(IJK))) CYCLE 66 IF (FLUID_AT(IJK)) THEN 67 SUM_VOL = SUM_VOL + VOL_V(IJK) 68 SUM_V_S = SUM_V_S + V_S(IJK,M)*EP_S(IJK,M)*VOL_V(IJK) 69 ENDIF 70 END DO 71 72 CALL GLOBAL_ALL_SUM(SUM_VOL) 73 CALL GLOBAL_ALL_SUM(SUM_V_S) 74 VAVG_V_S = SUM_V_S/SUM_VOL 75 ! 76 RETURN 77 END FUNCTION VAVG_V_S 78