vtc_scalar.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: CALC_VTC_SS(VxTC_ss, IER)                              C
!  Purpose: multiply solid-solid granular energy transfer coefficient  C
!           with cell volume
!                                                                      C
!  Author: J. Galvin                                 Date:             C
!  Reviewer: S. benyahia                             Date:             C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced:                                               C
!  Variables modified:                                                 C
!                                                                      C
!  Local variables:                                                    C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
!
      SUBROUTINE calc_vtc_ss(VXTC_SS)
!
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!-----------------------------------------------
!     Modules
!-----------------------------------------------
      USE param
      USE param1
      USE geometry
      USE indices
      USE physprop
      USE compar
      USE kintheory
      USE functions
      IMPLICIT NONE
!-----------------------------------------------
!     Local variables
!-----------------------------------------------

!                      Indices
      INTEGER          IJK
!
!                      Index of both solids phases
      INTEGER          L, M, LM
!
!                      Volume x interphase transfer coefficient
      DOUBLE PRECISION VxTC_ss(dimension_3, dimension_lm)
!-----------------------------------------------
!
      DO m = 1, mmax
          DO l = 1, mmax
               lm = funlm(l,m)
               IF (l .NE. m) THEN
!!!$omp  parallel do private(IJK)
                    DO ijk = ijkstart3, ijkend3
                         IF (fluid_at(ijk)) THEN
                              vxtc_ss(ijk,lm) = ed_ss_ip(ijk,lm)*vol(ijk)
                         ELSE
                              vxtc_ss(ijk,lm) = zero
                         ENDIF
                    ENDDO
               ENDIF
          ENDDO
      ENDDO
!
      RETURN
      END SUBROUTINE calc_vtc_ss
!-----------------------------------------------

!// Comments on the modifications for DMP version implementation
!// 001 Include header file and common declarations for parallelization
!// 350 Changed do loop limits: 1,ijkmax2-> ijkstart3, ijkend3