@@ -20,7 +20,6 @@ def write_rmg_input_file(rmg: dict,
2020 t3 : dict ,
2121 path : str ,
2222 iteration : int = 1 ,
23- run_sa : bool = False ,
2423 walltime : str = '00:00:00:00' ,
2524 ):
2625 """
@@ -33,7 +32,6 @@ def write_rmg_input_file(rmg: dict,
3332 path (str): The path where the RMG input file should be saved.
3433 iteration (int, optional): The T3 iteration, used to determine ``core_tolerance`` and
3534 ``tolerance_interrupt_simulation``. Does not matter for simulating or computing SA.
36- run_sa (bool, optional): Whether to add sensitivity analysis information into the RMG input file.
3735 walltime (str, optional): The time cap for an RMG run. Should pass here t3['options']['max_RMG_walltime']
3836 """
3937 rmg = rmg .copy ()
@@ -103,7 +101,7 @@ def write_rmg_input_file(rmg: dict,
103101 pressure=${pressure},
104102 initialMoleFractions={${concentrations()} },
105103 ${termination}
106- nSims=${conditions_per_iteration},${balance}${constant}${sensitivity}
104+ nSims=${conditions_per_iteration},${balance}${constant}
107105)
108106<%def name="concentrations()">
109107% for spc in species_list:
@@ -121,7 +119,7 @@ def write_rmg_input_file(rmg: dict,
121119 temperature=${temperature},
122120 initialConcentrations={${concentrations()} },
123121 ${termination}
124- nSims=${conditions_per_iteration},${constant}${sensitivity}
122+ nSims=${conditions_per_iteration},${constant}
125123)
126124<%def name="concentrations()">
127125% for spc in species_list:
@@ -169,16 +167,6 @@ def write_rmg_input_file(rmg: dict,
169167 constant = '\n constantSpecies=['
170168 constant += f"'{ spc ['label' ]} ', "
171169 constant += '],' if constant else ''
172- sensitivity = ''
173- observables = [spc ['label' ] for spc in species if spc ['observable' ]]
174- if run_sa and len (observables ):
175- sensitivity = '\n sensitivity=['
176- for i , observable in enumerate (observables ):
177- sensitivity += f"'{ observable } '"
178- if i < len (observables ) - 1 :
179- sensitivity += ', '
180- sensitivity += f"],\n sensitivityThreshold={ t3 ['sensitivity' ]['SA_threshold' ]} "
181-
182170
183171 if reactor ['type' ] == 'gas batch constant T P' :
184172 if isinstance (reactor ['P' ], float ):
@@ -201,7 +189,6 @@ def write_rmg_input_file(rmg: dict,
201189 conditions_per_iteration = reactor ['conditions_per_iteration' ],
202190 balance = balance ,
203191 constant = constant ,
204- sensitivity = sensitivity ,
205192 )
206193
207194 elif reactor ['type' ] == 'liquid batch constant T V' :
0 commit comments