File: N:\mfix\model\bodyforce_mod.f

1     !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
2     !                                                                      C
3     !  Module:  BODYFORCE                                                  C
4     !  Purpose: Include file for all body force statement functions        C
5     !                                                                      C
6     !  Author: M. Syamlal                                 Date:  6-MAR-92  C
7     !                                                                      C
8     !                                                                      C
9     !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
10     
11           MODULE bodyforce
12     
13           USE constant, only: gravity_x, gravity_y, gravity_z
14           IMPLICIT NONE
15     
16           CONTAINS
17     
18     ! Body force on gas at i+1/2, j, k
19           DOUBLE PRECISION FUNCTION BFX_g(IJK)
20           IMPLICIT NONE
21           INTEGER, INTENT(IN) :: ijk
22           BFX_g = GRAVITY_X
23           END FUNCTION BFX_g
24     
25     ! Body force on gas at i, j+1/2, k
26           DOUBLE PRECISION FUNCTION BFY_g(IJK)
27           IMPLICIT NONE
28           INTEGER, INTENT(IN) :: ijk
29           BFY_g = GRAVITY_Y
30           END FUNCTION BFY_g
31     
32     ! Body force on gas at i, j, k+1/2
33           DOUBLE PRECISION FUNCTION BFZ_g(IJK)
34           IMPLICIT NONE
35           INTEGER, INTENT(IN) :: ijk
36           BFZ_g = GRAVITY_Z
37           END FUNCTION BFZ_g
38     
39     ! Body force on solids m at i+1/2, j, k
40           DOUBLE PRECISION FUNCTION BFX_s(IJK,M)
41           IMPLICIT NONE
42           INTEGER, INTENT(IN) :: ijk,m
43           BFX_s = GRAVITY_X
44           END FUNCTION BFX_s
45     
46     ! Body force on solids m at i, j+1/2, k
47           DOUBLE PRECISION FUNCTION BFY_s(IJK,M)
48           IMPLICIT NONE
49           INTEGER, INTENT(IN) :: ijk,m
50           BFY_s = GRAVITY_Y
51           END FUNCTION BFY_s
52     
53     ! Body force on solids m at i, j, k+1/2
54           DOUBLE PRECISION FUNCTION BFZ_s(IJK,M)
55           IMPLICIT NONE
56           INTEGER, INTENT(IN) :: ijk,m
57           BFZ_s = GRAVITY_Z
58           END FUNCTION BFZ_s
59     
60           END MODULE bodyforce
61