139 subroutine initialize(this, chem_spec_data, aero_phase, aero_rep, n_cells)
150 integer(kind=i_kind),
intent(in) :: n_cells
152 type(
property_t),
pointer :: spec_props, reactants, products
153 character(len=:),
allocatable :: key_name, spec_name, string_val
154 integer(kind=i_kind) :: i_spec, i_qty
156 integer(kind=i_kind) :: temp_int
157 real(kind=
dp) :: temp_real
160 if (.not.
associated(this%property_set))
call die_msg(510658779, &
161 "Missing property set needed to initialize reaction")
162 key_name =
"reactants"
164 this%property_set%get_property_t(key_name, reactants), &
165 "Troe reaction is missing reactants")
166 key_name =
"products"
168 this%property_set%get_property_t(key_name, products), &
169 "Troe reaction is missing products")
172 call reactants%iter_reset()
176 call assert(844081716, reactants%get_property_t(val=spec_props))
178 if (spec_props%get_int(key_name, temp_int)) i_spec = i_spec+temp_int-1
179 call reactants%iter_next()
187 allocate(this%condensed_data_int(num_int_prop_ + &
188 (i_spec + 2) * (i_spec + products%size())))
189 allocate(this%condensed_data_real(num_real_prop_ + products%size()))
190 this%condensed_data_int(:) = int(0, kind=
i_kind)
191 this%condensed_data_real(:) = real(0.0, kind=
dp)
194 this%num_env_params = num_env_param_
199 num_prod_ = products%size()
202 conv_ =
const%avagadro /
const%univ_gas_const * 10.0d0**(-12.0d0)
208 if (.not. this%property_set%get_real(key_name, k0_a_))
then
212 if (.not. this%property_set%get_real(key_name, k0_b_))
then
216 if (.not. this%property_set%get_real(key_name, k0_c_))
then
220 if (.not. this%property_set%get_real(key_name, kinf_a_))
then
224 if (.not. this%property_set%get_real(key_name, kinf_b_))
then
228 if (.not. this%property_set%get_real(key_name, kinf_c_))
then
232 if (.not. this%property_set%get_real(key_name, fc_))
then
236 if (.not. this%property_set%get_real(key_name,
n_))
then
239 key_name =
"time unit"
240 scaling_ = real(1.0, kind=
dp)
241 if (this%property_set%get_string(key_name, string_val))
then
242 if (trim(string_val).eq.
"MIN")
then
243 scaling_ = real(1.0d0/60.0d0, kind=
dp)
248 k0_a_ = k0_a_ * real(1.0d6, kind=
dp)
251 call reactants%iter_reset()
256 react_(i_spec) = chem_spec_data%gas_state_id(
spec_name)
259 call assert_msg(595701751, react_(i_spec).gt.0, &
263 call assert(221292659, reactants%get_property_t(val=spec_props))
265 if (spec_props%get_int(key_name, temp_int))
then
266 do i_qty = 1, temp_int - 1
267 react_(i_spec + i_qty) = react_(i_spec)
269 i_spec = i_spec + temp_int - 1
272 call reactants%iter_next()
277 call products%iter_reset()
282 prod_(i_spec) = chem_spec_data%gas_state_id(
spec_name)
285 call assert_msg(480024633, prod_(i_spec).gt.0, &
289 call assert(393355097, products%get_property_t(val=spec_props))
291 if (spec_props%get_real(key_name, temp_real))
then
292 yield_(i_spec) = temp_real
297 call products%iter_next()