File: /nfs/home/0/users/jenkins/mfix.git/model/energy_mod.f
1 MODULE energy
2
3 USE fldvar
4 USE param
5 USE param1
6
7
8 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: HOR_g
9
10
11 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: HOR_s
12
13
14 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: GAMA_gs
15
16
17 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: GAMA_Rg
18
19
20 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: GAMA_Rs
21
22
23 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: T_Rg
24
25
26 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: T_Rs
27
28 CONTAINS
29
30
31
32
33
34
35
36
37
38 DOUBLE PRECISION FUNCTION S_Rpg(IJK)
39 IMPLICIT NONE
40 INTEGER IJK
41 S_Rpg = 4.d0 * GAMA_Rg(IJK) * T_g(IJK)**3
42 END FUNCTION S_Rpg
43
44
45 DOUBLE PRECISION FUNCTION S_Rcg(IJK)
46 IMPLICIT NONE
47 INTEGER IJK
48 S_Rcg = GAMA_Rg(IJK) * ( T_Rg(IJK)**4 + 3.d0 * T_g(IJK)**4 )
49 END FUNCTION S_Rcg
50
51
52 DOUBLE PRECISION FUNCTION S_Rps(IJK, M)
53 IMPLICIT NONE
54 INTEGER IJK, M
55 S_Rps = 4.d0 * GAMA_Rs(IJK, M) * T_s(IJK, M)**3
56 END FUNCTION S_Rps
57
58
59 DOUBLE PRECISION FUNCTION S_Rcs(IJK, M)
60 IMPLICIT NONE
61 INTEGER IJK, M
62 S_Rcs = GAMA_Rs(IJK, M) * ( T_Rs(IJK, M)**4 &
63 + 3.d0 * T_s(IJK, M)**4 )
64 END FUNCTION S_Rcs
65
66 END MODULE energy
67