search

agnwinds / python

This is the repository for Python, the radiative transfer code used to winds in AGN and other syatems
GNU General Public License v3.0
24 stars 25 forks source link

Atomic data - census #1062

Closed kslong closed 1 month ago

kslong commented 1 month ago

With an increasing number of macro atom datasets, we need a robust way to determine how many levels, collision x-sections, etc are involved. At present, python exits if one provides it with an atomic dataset with too many values for various things, but one is left to guess at how many values are needed, and python has to be recompiled after each guess.

Either Python needs a mode for counting the inputs of various types that have limits, or we need to write a helper script to do this assessesment.

kslong commented 1 month ago

A problem with doing this in atomic data is that there are many, many exit commands.

kslong commented 1 month ago

This problem is now addressed by a new routine called CheckAtomic.py, which describes the number of lements, ions, levels, and macro atoms for the files that are read by a masterfile. Additional information, for a file that can be run with Python is documented in the diag files.