1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Module name: CALC_VTC_SS(VxTC_ss, IER) C 4 ! Purpose: multiply solid-solid granular energy transfer coefficient C 5 ! with cell volume 6 ! C 7 ! Author: J. Galvin Date: C 8 ! Reviewer: S. benyahia Date: C 9 ! C 10 ! Literature/Document References: C 11 ! C 12 ! Variables referenced: C 13 ! Variables modified: C 14 ! C 15 ! Local variables: C 16 ! C 17 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 18 ! 19 SUBROUTINE CALC_VTC_SS(VXTC_SS) 20 ! 21 !...Translated by Pacific-Sierra Research VAST-90 2.06G5 12:17:31 12/09/98 22 !...Switches: -xf 23 ! 24 !----------------------------------------------- 25 ! Modules 26 !----------------------------------------------- 27 USE param 28 USE param1 29 USE geometry 30 USE indices 31 USE physprop 32 USE compar 33 USE kintheory 34 USE functions 35 IMPLICIT NONE 36 !----------------------------------------------- 37 ! Local variables 38 !----------------------------------------------- 39 40 ! Indices 41 INTEGER IJK 42 ! 43 ! Index of both solids phases 44 INTEGER L, M, LM 45 ! 46 ! Volume x interphase transfer coefficient 47 DOUBLE PRECISION VxTC_ss(DIMENSION_3, DIMENSION_LM) 48 !----------------------------------------------- 49 ! 50 DO M = 1, MMAX 51 DO L = 1, MMAX 52 LM = FUNLM(L,M) 53 IF (L .NE. M) THEN 54 !!!$omp parallel do private(IJK) 55 DO IJK = ijkstart3, ijkend3 56 IF (FLUID_AT(IJK)) THEN 57 VXTC_SS(IJK,LM) = ED_ss_ip(IJK,LM)*VOL(IJK) 58 ELSE 59 VXTC_SS(IJK,LM) = ZERO 60 ENDIF 61 ENDDO 62 ENDIF 63 ENDDO 64 ENDDO 65 ! 66 RETURN 67 END SUBROUTINE CALC_VTC_SS 68 !----------------------------------------------- 69 70 !// Comments on the modifications for DMP version implementation 71 !// 001 Include header file and common declarations for parallelization 72 !// 350 Changed do loop limits: 1,ijkmax2-> ijkstart3, ijkend3 73