We would like to use chiantipy for T_e=5000 to 10000 K

[Michihiro Takami]

Hi, I very much appreciate the development team for developing and maintaining this excellent software. Only the issue is chiantipy does not seem to work below T_e < 1e4 K. I am wondering if anyone could easily revise the software and allow us to make calculations for T_e=5000 to 10000 K. It would be tremendously useful for my student's MSc project. Many thanks for your consideration.

[Ken Dere]

you can do the following:

fe=ch.ioneq(26)

t=np.linspace(5000, 10000, 11)

fe.calculate(t)

fe.Ioneq.shape

Out[4]: (27, 11)

so that fe.Ioneq is the ionization equilibria for for the 27 fe ions (including the bare ion) for the 11 temperatures.

to use this you would probably need to create a some_name.ioneq files that can be read by ioneqRead().

Yes, it would be a bit of work and we should probably create  such a file ourselves.

hope this works,

Ken

[Michihiro Takami]

Tons of thanks for your quick advice, Ken! I have tested it and found that it works.

On the other hand. Here is our jupyter notebook.

https://www.dropbox.com/s/vm9x61c1iaxlsid/ChiantiPy%2C%20please%20check.ipynb?dl=0

We set the temperature in the 4th block of the commands.

Seems like we get an error in the 5th block of the temperature is set to be below 1e4 K.

I am wondering why. We very much appreciate your advice!

Thanks, Hiro

[Ken Dere]

the problem is that ionization equilibria are not available at that temperature.  the attribute .Intensity is a dict and here there is an errorMessage in it

I will see if I can put together a python script to run to create an ionization equilibria for a specified set of temperatures.

One thing that you can do is

o1.emiss() and then use o1.emissPlot()

that should at least get you started

regards,

Ken

[Michihiro Takami]

Many thanks for your advice, and sorry for my dalay. I finally had time to read your codes to find that o1.emiss stores the spectra in o1.Emiss['emiss'].
According you code, the unit is ergs s^-1 str^-1, and it does not include elemental abundance or ionization fraction. Does it mean the intensity/flux is given per atom/ion? I'd very much appreciate your advice. Thanks!

[Ken Dere]

Yes, that is true.

I do plan to create a jupyter ipython notebook for creating ionization equilibrium files for a specified set of temperatures but I have been busy

regards

[Michihiro Takami]

Thanks a lot for your reply, Ken! Please take your time for the jupyter notebook. My student still have a lot to do.