Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

PlasmaIonizationError: n_electron just turned "nan" - aborting #2484

Open
erins184 opened this issue Nov 30, 2023 · 1 comment
Open

PlasmaIonizationError: n_electron just turned "nan" - aborting #2484

erins184 opened this issue Nov 30, 2023 · 1 comment

Comments

@erins184
Copy link

erins184 commented Nov 30, 2023

I am attempting tp run TARDIS on a custom ejecta model from HESMA. I replaced the model section in the Tardis configuration file with csvy_model: . I created a sample csvy file from my HESMA file and have been attempting to run Tardis with it. However, when I execute the run_tardis command, I am getting this error: PlasmaIonizationError: n_electron just turned "nan" - aborting. I am getting the same error when I attempt to run TARDIS with the example csvy configuration file from the TARDIS website.

this is the line of code raising the error
sim = run_tardis("tardis_config_tester.yml",
virtual_packet_logging=True,
show_convergence_plots=True,
export_convergence_plots=True,
log_level="INFO")

csvy:
Screenshot 2023-11-29 at 8 28 09 PM

main Tardis config file:
Screenshot 2023-11-29 at 8 29 06 PM

Mac OS

Any help is greatly appreciated!!

@DeerWhale
Copy link
Contributor

I just tried with the example csvy model on the TARDIS website with the current version of TARDIS, and it seems working fine, could you try with the newest version of TARDIS?

If the issue still there, my guess would be the model setup, some of the reasons that could lead to the issue could be:

  • The inner boundary velocity (indicating where the photosphere is in the code) is not set properly, by default it uses the first shells available I think, and the first shell in your csvy is really deep in, it might be not enough energy packets escaping from the inner shells.
  • Helium is the only element in the composition but the helium_treatment is not set in the plasma configuration, so it could be no ionized atoms in the plasma, try add helium_treatment: recomb-nlte in the plasma part of the .yml file
  • Double check the all the values in the csvy model, the density value seems really high.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

3 participants