vavg_w_g.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: VAVG_W_g                                               C
!  Purpose: Calculate volume averaged W_g                              C
!                                                                      C
!  Author: M. Syamlal                                 Date: 28-APR-94  C
!  Reviewer:                                          Date: dd-mmm-yy  C
!                                                                      C
!  Revision Number:                                                    C
!  Purpose:                                                            C
!  Author:                                            Date: dd-mmm-yy  C
!  Reviewer:                                          Date: dd-mmm-yy  C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced:                                               C
!  Variables modified:                                                 C
!                                                                      C
!  Local variables:                                                    C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
!
      DOUBLE PRECISION FUNCTION vavg_w_g ()
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE bc
      USE compar
      USE fldvar
      USE functions
      USE geometry
      USE indices
      USE mpi_utility
      USE parallel
      USE param
      USE param1
      USE physprop
      USE run
      USE utilities, ONLY: mfix_isnan
      IMPLICIT NONE
!-----------------------------------------------
!   G l o b a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   L o c a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   L o c a l   V a r i a b l e s
!-----------------------------------------------
!
!                      Indices
      INTEGER          IJK
!
!                      Integral of W_g*EP_g for entire volume
      DOUBLE PRECISION SUM_W_g
!
!                      Total volume of computational cells
      DOUBLE PRECISION SUM_VOL

!  Integrate the velocity values for the whole domain,
!
      sum_w_g = zero
      sum_vol = zero

!!!!$omp   parallel do private(IJK) reduction(+:SUM_VOL,SUM_W_G)
      DO ijk = ijkstart3, ijkend3
      IF(.NOT.is_on_mype_wobnd(i_of(ijk), j_of(ijk), k_of(ijk))) cycle
         IF (fluid_at(ijk)) THEN
            sum_vol = sum_vol + vol_w(ijk)
            sum_w_g = sum_w_g + w_g(ijk)*ep_g(ijk)*vol_w(ijk)
         ENDIF
      END DO

      CALL global_all_sum(sum_vol)
      CALL global_all_sum(sum_w_g)
      vavg_w_g = sum_w_g/sum_vol
!
! uncomment the following lines to enable trapping NaN's.
!      IF( mfix_isnan(VAVG_W_G) ) THEN
!        write(*,*) VAVG_W_G,  ' NaN being caught in VAVG_W_G.f '
!        AUTOMATIC_RESTART = .TRUE.
!      ENDIF
!
      RETURN
      END FUNCTION vavg_w_g



      DOUBLE PRECISION FUNCTION vavg_flux_w_g ()
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE parallel
      USE run
      USE fldvar
      USE bc
      USE geometry
      USE physprop
      USE indices
      USE compar
      USE mpi_utility
      USE mflux
      USE functions
      IMPLICIT NONE
!-----------------------------------------------
!   G l o b a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   L o c a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   L o c a l   V a r i a b l e s
!-----------------------------------------------
!
!                      Indices
      INTEGER          IJK
!
!                      Integral of W_g*EP_g for entire volume
      DOUBLE PRECISION SUM_W_g
!
!                      Total volume of computational cells
      DOUBLE PRECISION SUM_AREA
!
!-----------------------------------------------
!
!  Integrate the velocity values for the whole domain,
!
      sum_w_g = zero
      sum_area = zero

!!!!$omp   parallel do private(IJK) reduction(+:SUM_AREA,SUM_W_G)
      DO ijk = ijkstart3, ijkend3
      IF(.NOT.is_on_mype_wobnd(i_of(ijk), j_of(ijk), k_of(ijk))) cycle
         IF (fluid_at(ijk)) THEN
           IF(.NOT.added_mass) THEN
              sum_w_g = sum_w_g + flux_gt(ijk)
           ELSE
              sum_w_g = sum_w_g + flux_gst(ijk)
           ENDIF
           sum_area = sum_area + axy(ijk)
         ENDIF
      END DO

      CALL global_all_sum(sum_area)
      CALL global_all_sum(sum_w_g)
      vavg_flux_w_g = sum_w_g/sum_area
!
      RETURN
      END FUNCTION vavg_flux_w_g