Skip to content

Commit 72f8896

Browse files
committed
Chemkin read: cope with reactions with 4 reactants
It needs fifth order units to be specified also to prepare for pressure-dependent reactions.
1 parent ec63dfb commit 72f8896

1 file changed

Lines changed: 2 additions & 0 deletions

File tree

rmgpy/chemkin.pyx

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1405,13 +1405,15 @@ def read_reactions_block(f, species_dict, read_comments=True):
14051405
'{0}^3/({1}*{2})'.format(volume_units, molecule_units, time_units), # Second-order
14061406
'{0}^6/({1}^2*{2})'.format(volume_units, molecule_units, time_units), # Third-order
14071407
'{0}^9/({1}^3*{2})'.format(volume_units, molecule_units, time_units), # Fourth-order
1408+
f"{volume_units}^12/({molecule_units}^4*{time_units})", # Fifth-order
14081409
]
14091410

14101411
Aunits_surf = [
14111412
'', # Zeroth-order
14121413
's^-1'.format(time_units), # First-order
14131414
'{0}^2/({1}*{2})'.format(area_units, molecule_units, time_units), # Second-order
14141415
'{0}^4/({1}^2*{2})'.format(area_units, molecule_units, time_units), # Third-order
1416+
'{0}^6/({1}^3*{2})'.format(area_units, molecule_units, time_units), # Fourth-order
14151417
]
14161418
Eunits = energy_units
14171419

0 commit comments

Comments
 (0)