-
Notifications
You must be signed in to change notification settings - Fork 33
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
unit of the radiative loss rate #56
Comments
Yes I think you are correct in that the label on this plot is not correct as this gives units of erg s^-1 cm^5 The radiative loss rate is typically reported in units of erg cm^3 s^-1 such that I believe the free-bound loss should also have units of erg cm^3 s^-1 and so you would be correct on that count as well. @kdere can you confirm this? |
I agree that the label is not correct. What is plotted is Lambda and Lambdan_en_H gives the radiative loss rate per unit volume erg cm^-3 s^-1. The units for Lambda are erg cm^3 s^-1, and the same goes for the free-bound loss rate |
Ok so it seems that the units listed in the plot labels and docstrings in several places need to be fixed. @smoon24 would you mind listing all of the places where you've found these inconsistencies and adding to it as you find more? These could all be fixed in a single PR. |
I used Chianti for coursework by now, so I'm afraid whether if I can find more. I just have checked functions related to RadLoss and found that boundBoundLoss method in Ion.py seems to has similar error. Docstring says it returns loss rate in unit of erg cm^-3 s^-1, which should be erg cm^3 s^-1. And it might be better to add unit in the docstring of freeFreeLoss and twoPhotonLoss for consistency. Thanks! |
Thanks for your feedback @smoon24. If you happen to come across any more of these inconsistencies, feel free to list them here. We'll keep this issue open until these get fixed. |
Yes, this is one of the things I need to get to. I think it is also on the QuickStart page and notebook |
The unit of the radiative loss rate in the document looks strange to me. I think it should be something like$\Lambda/(n_e n_H)$ where \Lambda is usual volumetric cooling function.
Also, the docstring in freeBoundLoss method says it returns the cooling rate in units of ergs cm^-2 s^-1 which is again conflicting with the y-label on above plot.
The text was updated successfully, but these errors were encountered: