1 | include(dom.inc)
2 |
3 | ! ----------!
4 | ! Constants !
5 | ! ----------!
6 |
7 | DOM_REAL, PARAMETER :: sigma = 5.667E-08 ! [WK4/m2]
8 | DOM_REAL, PARAMETER :: pi = 3.141592654
9 | DOM_REAL, PARAMETER :: PASCAL = 9.8692E-6 ! [Atm]
10 | DOM_REAL, PARAMETER :: c0 = 2.9979E8 ! [m/s]
11 | DOM_REAL, PARAMETER :: Boltzmann= 1.3805E-23 ! [J/K]
12 | DOM_REAL, PARAMETER :: Planckcst= 6.626E-34 ! [J.s]
13 |
14 | ! ------------------!
15 | ! Generic constants !
16 | ! ------------------!
17 |
18 | DOM_REAL, PARAMETER :: zero = 0.
19 | DOM_REAL, PARAMETER :: half = 0.5
20 | DOM_REAL, PARAMETER :: small = 1.0d-12
21 | DOM_REAL, PARAMETER :: big = 1.0d12
22 |
23 | ! Memory needed for 10000 cells
24 | DOM_REAL, PARAMETER :: Memtot = 50 ! [Mb]
25 |
26 | ! -------------------!
27 | ! Species properties !
28 | ! -------------------!
29 | DOM_REAL, PARAMETER :: W_H2O = 0.01801530
30 | DOM_REAL, PARAMETER :: W_CO2 = 0.04400980
31 | DOM_REAL, PARAMETER :: W_CO = 0.02801400
32 | DOM_REAL, PARAMETER :: W_O2 = 0.03199880
33 | DOM_REAL, PARAMETER :: W_N2 = 0.02801340
34 | DOM_REAL, PARAMETER :: W_C3H8 = 0.04409721
35 | DOM_REAL, PARAMETER :: W_CH4 = 0.01604230
36 |
37 | ! --------------!
38 | ! Element types !
39 | ! --------------!
40 | DOM_INT, PARAMETER :: EL_TRI = 1
41 | DOM_INT, PARAMETER :: EL_QUAD = 2
42 | DOM_INT, PARAMETER :: EL_TETRA = 3
43 | DOM_INT, PARAMETER :: EL_PYRAM = 4
44 | DOM_INT, PARAMETER :: EL_PRISM = 5
45 | DOM_INT, PARAMETER :: EL_HEXA = 6
46 |
47 | ! -------------------------!
48 | ! Element types for Gambit !
49 | ! -------------------------!
50 | DOM_INT, PARAMETER :: GB_QUAD = 2
51 | DOM_INT, PARAMETER :: GB_TRI = 3
52 | DOM_INT, PARAMETER :: GB_HEXA = 4
53 | DOM_INT, PARAMETER :: GB_PRISM = 5
54 | DOM_INT, PARAMETER :: GB_TETRA = 6
55 | DOM_INT, PARAMETER :: GB_PYRAM = 7
56 |
57 | ! -----------!
58 | ! MAX values !
59 | ! -----------!
60 |
61 | DOM_INT, PARAMETER :: MAX_NFACES_CELL = 8
62 | DOM_INT, PARAMETER :: MAX_NNODES_CELL = 10
63 |
64 | ! ---------------------------------!
65 | ! Quadrature discretization method !
66 | ! ---------------------------------!
67 |
68 | CHARACTER*5, PARAMETER :: SNDOM='SNDOM'
69 | CHARACTER*5, PARAMETER :: TNDOM='TNDOM'
70 | CHARACTER*3, PARAMETER :: FVM ='FVM'
71 |
72 | CHARACTER*4 , PARAMETER :: EXPON='EXPO'
73 | CHARACTER*4 , PARAMETER :: DMFS ='DMFS'
74 | CHARACTER*4 , PARAMETER :: SMFS ='STEP'
75 |
76 | ! ------------!
77 | ! IN/OUT file !
78 | ! ------------!
79 | DOM_INT, PARAMETER :: FILE_NODES = 100
80 | DOM_INT, PARAMETER :: FILE_CELLS = 101
81 | DOM_INT, PARAMETER :: FILE_FACES = 102
82 | DOM_INT, PARAMETER :: FILE_C2C = 103
83 | DOM_INT, PARAMETER :: FILE_CCELL = 104
84 | DOM_INT, PARAMETER :: FILE_CFACE = 105
85 | DOM_INT, PARAMETER :: FILE_CLNOD = 106
86 | DOM_INT, PARAMETER :: FILE_NORM = 107
87 | DOM_INT, PARAMETER :: FILE_VOLA = 108
88 | DOM_INT, PARAMETER :: FILE_EMIS = 109
89 | DOM_INT, PARAMETER :: FILE_KSCA = 110
90 | DOM_INT, PARAMETER :: FILE_EXTRN = 111
91 | DOM_INT, PARAMETER :: FILE_QUADR = 112
92 | DOM_INT, PARAMETER :: FILE_CLPRO = 113
93 | DOM_INT, PARAMETER :: FILE_CLFAC = 114
94 | DOM_INT, PARAMETER :: FILE_PROG = 115
95 | DOM_INT, PARAMETER :: FILE_SNBWN = 116
96 | DOM_INT, PARAMETER :: FILE_LSPEC = 117
97 | DOM_INT, PARAMETER :: FILE_SSPEC = 118
98 | DOM_INT, PARAMETER :: FILE_G = 119
99 | DOM_INT, PARAMETER :: FILE_Sr = 120
100 | DOM_INT, PARAMETER :: FILE_H = 121
101 | DOM_INT, PARAMETER :: FILE_Qw = 122
102 | DOM_INT, PARAMETER :: FILE_Qr = 123
103 | DOM_INT, PARAMETER :: FILE_C2FAC = 124
104 | DOM_INT, PARAMETER :: FILE_FAC = 125
105 | DOM_INT, PARAMETER :: FILE_T = 126
106 | DOM_INT, PARAMETER :: FILE_E = 127
107 | DOM_INT, PARAMETER :: FILE_GDATA = 128
108 | DOM_INT, PARAMETER :: FILE_PROP = 129
109 | DOM_INT, PARAMETER :: FILE_CHCS = 130
110 | DOM_INT, PARAMETER :: FILE_Tf = 131
111 | DOM_INT, PARAMETER :: FILE_EMIS2 = 132
112 | DOM_INT, PARAMETER :: FILE_G2 = 133
113 | DOM_INT, PARAMETER :: FILE_Sr2 = 134
114 | DOM_INT, PARAMETER :: FILE_H2 = 135
115 | DOM_INT, PARAMETER :: FILE_Qw2 = 136
116 | DOM_INT, PARAMETER :: FILE_Qr2 = 137
117 | DOM_INT, PARAMETER :: FILE_GDAT2 = 138
118 | DOM_INT, PARAMETER :: FILE_MSH = 139
119 | DOM_INT, PARAMETER :: FILE_KEXT = 140
120 | DOM_INT, PARAMETER :: FILE_Kp = 141
121 | DOM_INT, PARAMETER :: FILE_Kp2 = 142
122 | DOM_INT, PARAMETER :: FILE_Wsgg = 143
123 | DOM_INT, PARAMETER :: FILE_CLDAT = 144
124 | DOM_INT, PARAMETER :: FILE_FC = 145
125 | DOM_INT, PARAMETER :: FILE_TRACK = 146
126 | DOM_INT, PARAMETER :: FILE_QRP = 147
127 | DOM_INT, PARAMETER :: FILE_TAB = 148
128 |
129 | ! ---------------!
130 | ! Gaz identifier !
131 | ! ---------------!
132 | DOM_INT, PARAMETER :: GAZ_CO = 1
133 | DOM_INT, PARAMETER :: GAZ_C = 2
134 | DOM_INT, PARAMETER :: GAZ_H = 3