From fcaa1f60aa67c4c716dd311155cf36a8c7af1cd4 Mon Sep 17 00:00:00 2001 From: Atharva Arya Date: Fri, 27 Dec 2024 01:28:53 -0500 Subject: [PATCH] standard workflow debugging 10 --- tardis/io/model/parse_radiation_field_configuration.py | 6 ++++++ tardis/plasma/standard_plasmas.py | 9 ++++++++- 2 files changed, 14 insertions(+), 1 deletion(-) diff --git a/tardis/io/model/parse_radiation_field_configuration.py b/tardis/io/model/parse_radiation_field_configuration.py index d0ff62badc9..1020457e7fa 100644 --- a/tardis/io/model/parse_radiation_field_configuration.py +++ b/tardis/io/model/parse_radiation_field_configuration.py @@ -126,6 +126,12 @@ def parse_radiation_field_state_from_csvy( dilution_factor = csvy_model_data["dilution_factor"].iloc[1:].values else: dilution_factor = calculate_geometric_dilution_factor(geometry) + + print("\n=== DilutePlanckianRadiationField in parse_radiation_field_configuration/parse_radiation_field_state_from_csvy===") + print(f"t_radiative = {t_radiative}") + print(f"dilution_factor = {dilution_factor}") + print(f"geometry = {geometry}") + print("="*50) return DilutePlanckianRadiationField(t_radiative, dilution_factor, geometry) diff --git a/tardis/plasma/standard_plasmas.py b/tardis/plasma/standard_plasmas.py index 2f9267d80ea..42bc6e6bf53 100644 --- a/tardis/plasma/standard_plasmas.py +++ b/tardis/plasma/standard_plasmas.py @@ -114,7 +114,11 @@ def assemble_plasma(config, simulation_state, atom_data=None): species_string_to_tuple(species) for species in config.plasma.nlte_excitation_species ] - + print("\n=== DilutePlanckianRadiationField in standard_plasmas===") + print(f"simulation_state.t_radiative = {simulation_state.t_radiative}") + print(f"simulation_state.dilution_factor = {simulation_state.dilution_factor}") + print("exiting standard_plasmas") + print("="*50) dilute_planckian_radiation_field = DilutePlanckianRadiationField( simulation_state.t_radiative, simulation_state.dilution_factor ) @@ -211,6 +215,9 @@ def assemble_plasma(config, simulation_state, atom_data=None): config.plasma.link_t_rad_t_electron * dilute_planckian_radiation_field.temperature.to(u.K).value ) + print("\n=== t_electrons in standard_plasmas===") + print(f"t_electrons = {t_electrons}") + print("="*50) initial_continuum_solver = DiluteBlackBodyContinuumPropertiesSolver( atom_data )