File: N:\mfix\model\des\comp_mean_fields.f

1     !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
2     !                                                                      !
3     !  Subroutine: COMP_MEAN_FIELDS                                        !
4     !  Author: J.Musser                                   Date: 11-NOV-14  !
5     !                                                                      !
6     !  Purpose: Driver routine for calculating continuous field variables  !
7     !  corresponding to discrete data (ROP_s, u_s, v_s, w_s)               !
8     !                                                                      !
9     !  o The diffusion filter is only applied to the the solids bulk       !
10     !    density because DEM simulations do not utilize the other field    !
11     !    variables within a time loop.                                     !
12     !                                                                      !
13     !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
14           SUBROUTINE COMP_MEAN_FIELDS
15     
16     ! Modules
17     !---------------------------------------------------------------------//
18           use discretelement, only: DES_MMAX
19           use fldvar, only: rop_s
20     ! Flag: Diffuse DES field variables.
21           use particle_filter, only: DES_DIFFUSE_MEAN_FIELDS
22           use particle_filter, only: DES_INTERP_MEAN_FIELDS
23           use particle_filter, only: DES_INTERP_SCHEME_ENUM
24           use particle_filter, only: DES_INTERP_NONE
25           use particle_filter, only: DES_INTERP_GARG
26           use physprop, only: mmax
27           IMPLICIT NONE
28     
29     ! Local variables
30     !---------------------------------------------------------------------//
31     
32     ! Loop counter.
33           INTEGER :: M
34     
35     !......................................................................!
36     
37     ! Calculate field variables from particle data:
38           IF(DES_INTERP_MEAN_FIELDS) THEN
39              SELECT CASE(DES_INTERP_SCHEME_ENUM)
40              CASE(DES_INTERP_GARG) ; CALL COMP_MEAN_FIELDS0
41              CASE DEFAULT; CALL COMP_MEAN_FIELDS1
42              END SELECT
43           ELSE
44              CALL COMP_MEAN_FIELDS1
45           ENDIF
46     
47     ! Apply the diffusion filter.
48           IF(DES_DIFFUSE_MEAN_FIELDS) THEN
49              DO M=MMAX+1, MMAX+DES_MMAX
50                 CALL DIFFUSE_MEAN_FIELD(ROP_S(:,M),'ROP_S')
51              ENDDO
52           ENDIF
53     
54     ! Calculate the gas phase volume fraction from ROP_s.
55           CALL CALC_EPG_DES
56     
57           RETURN
58           END SUBROUTINE COMP_MEAN_FIELDS
59