comp_mean_fields.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
!                                                                      !
!  Subroutine: COMP_MEAN_FIELDS                                        !
!  Author: J.Musser                                   Date: 11-NOV-14  !
!                                                                      !
!  Purpose: Driver routine for calculating continuous field variables  !
!  corresponding to discrete data (ROP_s, u_s, v_s, w_s)               !
!                                                                      !
!  o The diffusion filter is only applied to the the solids bulk       !
!    density because DEM simulations do not utilize the other field    !
!    variables within a time loop.                                     !
!                                                                      !
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
      SUBROUTINE comp_mean_fields

! Modules
!---------------------------------------------------------------------//
      use discretelement, only: des_mmax
      use fldvar, only: rop_s
! Flag: Diffuse DES field variables.
      use particle_filter, only: des_diffuse_mean_fields
      use particle_filter, only: des_interp_mean_fields
      use particle_filter, only: des_interp_scheme_enum
      use particle_filter, only: des_interp_none
      use particle_filter, only: des_interp_garg
      use physprop, only: mmax
      IMPLICIT NONE

! Local variables
!---------------------------------------------------------------------//

! Loop counter.
      INTEGER :: M

!......................................................................!

! Calculate field variables from particle data:
      IF(des_interp_mean_fields) THEN
         SELECT CASE(des_interp_scheme_enum)
         CASE(des_interp_garg) ; CALL comp_mean_fields0
         CASE DEFAULT; CALL comp_mean_fields1
         END SELECT
      ELSE
         CALL comp_mean_fields1
      ENDIF

! Apply the diffusion filter.
      IF(des_diffuse_mean_fields) THEN
         DO m=mmax+1, mmax+des_mmax
            CALL diffuse_mean_field(rop_s(:,m),'ROP_S')
         ENDDO
      ENDIF

! Calculate the gas phase volume fraction from ROP_s.
      CALL calc_epg_des

      RETURN
      END SUBROUTINE comp_mean_fields