@@ -701,10 +701,10 @@ def initialize(self, **kwargs):
701701 # This is necessary so that the PDep algorithm can identify the bath gas
702702 for spec in self .initial_species :
703703 if not spec .reactive :
704- self .reaction_model .enlarge (spec )
704+ self .reaction_model .enlarge (spec , requires_rms = requires_rms )
705705 for spec in self .initial_species :
706706 if spec .reactive :
707- self .reaction_model .enlarge (spec )
707+ self .reaction_model .enlarge (spec , requires_rms = requires_rms )
708708
709709 # chatelak: store constant SPC indices in the reactor attributes if any constant SPC provided in the input file
710710 # advantages to write it here: this is run only once (as species indexes does not change over the generation)
@@ -715,7 +715,7 @@ def initialize(self, **kwargs):
715715
716716 self .initialize_reaction_threshold_and_react_flags ()
717717 if self .filter_reactions and self .init_react_tuples :
718- self .react_init_tuples ()
718+ self .react_init_tuples (requires_rms = requires_rms )
719719 self .reaction_model .initialize_index_species_dict ()
720720
721721 self .initialize_seed_mech ()
@@ -828,6 +828,7 @@ def execute(self, initialize=True, **kwargs):
828828 unimolecular_react = self .unimolecular_react ,
829829 bimolecular_react = self .bimolecular_react ,
830830 trimolecular_react = self .trimolecular_react ,
831+ requires_rms = requires_rms ,
831832 )
832833
833834 if not np .isinf (self .model_settings_list [0 ].thermo_tol_keep_spc_in_edge ):
@@ -1899,7 +1900,7 @@ def initialize_reaction_threshold_and_react_flags(self):
18991900 if self .trimolecular :
19001901 self .trimolecular_react [:num_restart_spcs , :num_restart_spcs , :num_restart_spcs ] = False
19011902
1902- def react_init_tuples (self ):
1903+ def react_init_tuples (self , requires_rms = False ):
19031904 """
19041905 Reacts tuples given in the react block
19051906 """
@@ -1932,6 +1933,7 @@ def react_init_tuples(self):
19321933 unimolecular_react = self .unimolecular_react ,
19331934 bimolecular_react = self .bimolecular_react ,
19341935 trimolecular_react = self .trimolecular_react ,
1936+ requires_rms = requires_rms ,
19351937 )
19361938
19371939 def update_reaction_threshold_and_react_flags (
0 commit comments