@@ -99,7 +99,10 @@ def generate_concentrations(self) -> None:
9999 self .time_track (start , 'Concentrations' )
100100 return
101101
102- def _evaluate_conc (self , cur_conc : float , cur_p_conc : float , lam_p : float , lam : float , ti : int , dt : list [float ], fission_rates : list [float ], concs : list [float ], p_concs : list [float ], y_p : float , y : float ) -> tuple [float , float ]:
102+ def _evaluate_conc (self , cur_conc : float , cur_p_conc : float , lam_p : float ,
103+ lam : float , ti : int , dt : list [float ],
104+ fission_rates : list [float ], concs : list [float ],
105+ p_concs : list [float ], y_p : float , y : float ) -> tuple [float , float ]:
103106 """
104107 Evaluates the concentration of a nuclide and its decay parent.
105108 Yield of parent is assumed to be scaled by branching ratio.
@@ -197,7 +200,10 @@ def _add_debug_dnp_data(self, data: list[dict[str, float]]) -> list[dict[str, fl
197200 y_p = 0
198201 lam_p = 1
199202 for ti , t in enumerate (times [:- 1 ]):
200- cur_conc , cur_p_conc = self ._evaluate_conc (cur_conc , cur_p_conc , lam_p , lam , ti , dt , fission_rates , concs , p_concs , y_p , y )
203+ cur_conc , cur_p_conc = self ._evaluate_conc (cur_conc , cur_p_conc ,
204+ lam_p , lam , ti , dt ,
205+ fission_rates , concs ,
206+ p_concs , y_p , y )
201207 p_concs .append (cur_p_conc )
202208 concs .append (cur_conc )
203209 data .append (
@@ -274,8 +280,10 @@ def OMC_concentrations(self) -> list[dict[str, float]]:
274280 env = Environment (loader = PackageLoader ('mosden' ))
275281 file = self .openmc_settings ['omc_file' ]
276282 template = env .get_template (file )
277- chain_file = os .path .join (self .unprocessed_data_dir , self .openmc_settings ['chain' ])
278- cross_sections = os .path .join (self .unprocessed_data_dir , self .openmc_settings ['x_sections' ])
283+ chain_file = os .path .join (self .unprocessed_data_dir ,
284+ self .openmc_settings ['chain' ])
285+ cross_sections = os .path .join (self .unprocessed_data_dir ,
286+ self .openmc_settings ['x_sections' ])
279287 omc_dir = self .openmc_settings ['omc_dir' ]
280288 time_rate_data = self ._get_times_and_rates (self .f_in )
281289 render_data = {
@@ -407,7 +415,10 @@ def _collect_omc_fissions(self) -> tuple[dict[str, np.ndarray], np.ndarray]:
407415 for tally in sp .tallies .keys ():
408416 tally_data = sp .get_tally (id = tally )
409417 if 'fissionrate' in tally_data .name :
410- df = tally_data .get_pandas_dataframe (filters = False , scores = False , derivative = False , paths = False )
418+ df = tally_data .get_pandas_dataframe (filters = False ,
419+ scores = False ,
420+ derivative = False ,
421+ paths = False )
411422 try :
412423 df ['mean' ] = df ['mean' ] * self .openmc_settings ['source' ]
413424 except KeyError :
@@ -485,7 +496,10 @@ def _collect_omc_nuyield(self) -> dict[str, dict[str, np.ndarray]]:
485496 tally_data = sp .get_tally (id = tally )
486497 if 'delnuyield' in tally_data .name :
487498 use_dict = nuyield ['d' ]
488- df = tally_data .get_pandas_dataframe (filters = False , scores = False , derivative = False , paths = False )
499+ df = tally_data .get_pandas_dataframe (filters = False ,
500+ scores = False ,
501+ derivative = False ,
502+ paths = False )
489503 try :
490504 df ['mean' ] = (df ['mean' ] * self .openmc_settings ['source' ] / avg_fission_rate )
491505 except KeyError :
@@ -504,7 +518,10 @@ def _collect_omc_nuyield(self) -> dict[str, dict[str, np.ndarray]]:
504518
505519 if 'pmtnuyield' in tally_data .name :
506520 use_dict = nuyield ['p' ]
507- df = tally_data .get_pandas_dataframe (filters = False , scores = False , derivative = False , paths = False )
521+ df = tally_data .get_pandas_dataframe (filters = False ,
522+ scores = False ,
523+ derivative = False ,
524+ paths = False )
508525 try :
509526 df ['mean' ] = (df ['mean' ] * self .openmc_settings ['source' ] / avg_fission_rate )
510527 except KeyError :
0 commit comments