DH_t DH_t_0*exp(DH_t_opt) "Heat of reaction wood -> tar"
DH_g DH_t "Heat of reaction wood -> gas"
DH_is DH_t "Heat of reaction wood -> intermediate solid"
DH_c2 -42[kJ/kg] "Heat of reaction tar -> char"
DH_g2 DH_c2 "Heat of reaction tar -> gas"
DH_c DH_c_0*exp(DH_c_opt) "Heat of reaction intermediate solid -> char"
A_t 1.08e10[1/s] "Frequency factor w (wood) -> t (tar)"
E_t 148[kJ/mol] "Activation energy w -> t"
A_is A_is_0*exp(A_is_opt) "Frequency factor w -> is (intermediate solid)"
E_is 111.7[kJ/mol] "Activation energy w -> is"
A_g 4.38e9[1/s] "Frequency factor w -> g (primary gas)"
E_g 152.7[kJ/mol] "Activation energy w -> g"
A_c 1.38e10[1/s] "Frequency factor is-> c (char)"
E_c 161[kJ/mol] "Activation energy is -> c"
A_c2 1.00e5[1/s] "Frequency factor t -> c"
E_c2 108[kJ/mol] "Activation energy t  -> c"
A_g2 4.28e6[1/s] "Frequency factor t -> g"
E_g2 108[kJ/mol] "Activation energy t -> g"

A_is_0 1e7[1/s] "Initial value for physical parameter A_is"
A_is_lower 1e4[1/s] "Lower bound for physical parameter A_is"
A_is_upper 1e9[1/s] "Upper bound for physical parameter A_is"
A_is_opt 0 "Control parameter to be optimized"
A_is_opt_0 log(A_is_0/A_is_0) "Initial value for A_is_opt"
A_is_opt_lower log(A_is_lower/A_is_0) "Lower bound for A_is_opt"
A_is_opt_upper log(A_is_upper/A_is_0) "Upper bound for A_is_opt"

DH_t_0 50[kJ/kg] "Initial value for physical parameter DH_t"
DH_t_lower 1[kJ/kg] "Lower bound for physical parameter DH_t"
DH_t_upper 400[kJ/kg] "Upper bound for physical parameter DH_t"
DH_t_opt 0 "Control parameter to be optimized"
DH_t_opt_0_ log(DH_t_0/DH_t_0) "Initial value for DH_t_opt"
DH_t_opt_lower log(DH_t_lower/DH_t_0) "Lower bound for DH_t_opt"
DH_t_opt_upper log(DH_t_upper/DH_t_0) "Upper bound for DH_t_opt"

DH_c_0 -200[kJ/kg] "Initial value for physical parameter DH_c"
DH_c_lower -400[kJ/kg] "Lower bound for physical parameter DH_c"
DH_c_upper -10[kJ/kg] "Upper bound for physical parameter DH_c"
DH_c_opt 0 "Control parameter to be optimized"
DH_c_opt_0_ log(DH_c_0/DH_c_0) "Initial value for DH_c_opt"
DH_c_opt_lower log(DH_c_upper/DH_c_0) "Lower bound for DH_c_opt"
DH_c_opt_upper log(DH_c_lower/DH_c_0) "Upper bound for DH_c_opt"