1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Module name: Scalar_PROP(IER) C 4 ! Purpose: Calculate diffusion coefficeint and sources for user-defined 5 ! scalars 6 ! C 7 ! Author: Date: C 8 ! Reviewer: Date: C 9 ! C 10 ! C 11 ! Literature/Document References: C 12 ! C 13 ! Variables referenced: None C 14 ! Variables modified: None C 15 ! C 16 ! Local variables: None C 17 ! C 18 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 19 ! 20 SUBROUTINE SCALAR_PROP() 21 !...Translated by Pacific-Sierra Research VAST-90 2.06G5 12:17:31 12/09/98 22 !...Switches: -xf 23 ! 24 !----------------------------------------------- 25 ! M o d u l e s 26 !----------------------------------------------- 27 USE param 28 USE param1 29 USE fldvar 30 USE physprop 31 USE geometry 32 USE indices 33 USE run 34 USE scalars 35 USE toleranc 36 USE compar 37 USE sendrecv 38 USE functions 39 IMPLICIT NONE 40 !----------------------------------------------- 41 ! G l o b a l P a r a m e t e r s 42 !----------------------------------------------- 43 !----------------------------------------------- 44 ! D u m m y A r g u m e n t s 45 !----------------------------------------------- 46 ! 47 INTEGER L,IJK 48 ! 49 !----------------------------------------------- 50 51 IF(NScalar == 0) RETURN 52 ! 53 ! --- Remember to include all the local variables here for parallel 54 ! ---- processing 55 !!!$omp parallel do private(ijk, L) 56 DO IJK = IJKSTART3, IJKEND3 57 IF (FLUID_AT(IJK)) THEN 58 DO L = 1, NScalar 59 60 ! d (Scalar)/dt = S 61 ! S is linearized as S = Scalar_c - Scalar_p * Scalar 62 ! Scalar_c and Scalar_p must be >= 0 63 ! *** Uncomment next two lines *** 64 Scalar_c (IJK, L) = ZERO 65 Scalar_p (IJK, L) = ZERO 66 ! 67 ! Diffusion coefficient for User-defined Scalars 68 ! *** Uncomment next one line *** 69 Dif_Scalar(IJK, L) =ZERO 70 END DO 71 ! 72 ENDIF 73 END DO 74 !\\Sendrecv operations - just to make sure all the variables computed are 75 ! are passed and updated locally - fool-proof approach - Sreekanth - 102199 76 77 ! call send_recv(Scalar_c,2) 78 ! call send_recv(Scalar_p,2) 79 ! call send_recv(Dif_Scalar,2) 80 RETURN 81 END SUBROUTINE SCALAR_PROP 82