Greetings,
I run into the errors below about undeclared duplicate reactions. As a test, if I comment out the first reaction (H2 + M <=> 2 H + M), other reactions appear.
Any thoughts on this? C3MechV3.3 is a massive mechanism, so a script to fix this would be appreciated.
ck2yaml --permissive --input C3MechV3.3.MECH --thermo C3MechV3.3.THERM --transport C3MechV3.3.TRAN
Wrote YAML mechanism file to 'C3MechV3.3.yaml'.
Mechanism contains 3761 species and 16532 reactions.
Validating mechanism...
/opt/anaconda3/envs/ct-env/lib/python3.11/site-packages/cantera/ck2yaml.py:2249: UserWarning: NasaPoly2::validate:
For species OHV, discontinuity in h/RT detected at Tmid = 1000
Value computed using low-temperature polynomial: 53.62056162666667
Value computed using high-temperature polynomial: 53.5841554314
gas = ct.Solution(out_name)
/opt/anaconda3/envs/ct-env/lib/python3.11/site-packages/cantera/ck2yaml.py:2249: UserWarning: NasaPoly2::validate:
For species CHV, discontinuity in h/RT detected at Tmid = 1000
Value computed using low-temperature polynomial: 107.5046684
Value computed using high-temperature polynomial: 107.34847808033334
gas = ct.Solution(out_name)
/opt/anaconda3/envs/ct-env/lib/python3.11/site-packages/cantera/ck2yaml.py:2249: UserWarning: NasaPoly2::validate:
For species KCHOCH2CO3H, discontinuity in cp/R detected at Tmid = 1000
Value computed using low-temperature polynomial: 22.934254100000004
Value computed using high-temperature polynomial: 23.200262099999996
gas = ct.Solution(out_name)
FAILED
Undeclared duplicate reaction H2 + M <=> 2 H + M
found on lines 1303 and 1314 of the kinetics input file.