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

Parameters of XS_ismabs #171

Open
rwiller opened this issue Oct 7, 2021 · 3 comments
Open

Parameters of XS_ismabs #171

rwiller opened this issue Oct 7, 2021 · 3 comments
Assignees

Comments

@rwiller
Copy link
Contributor

rwiller commented Oct 7, 2021

The model XS_ismabs as build by astromodels only has 29 out of the required 31 parameters. I believ the reason is that the abundances of S II and S III are missing due to their names being taken by Si I (sii in the code) and Si II (siii). I suggest to use the name su instead of s to denote the ionization states of sulphur. Maybe consider naming silicon differently as well to avoid confusion with the indication of the ionization state.

@grburgess
Copy link
Contributor

Thanks @rwiller I'm glad you found a work around! I think that we will wait for the newer interface to XSPEC to go online rather than looking for a fix to this. I will keep open just in case anyone else runs into the same problem . Feel free to post your workaround here if you think it would be helpful!

@rwiller
Copy link
Contributor Author

rwiller commented Oct 8, 2021

Right, in the file XS_ismabs.py, which is created by astromodels (located in /home//.astromodels/data) I just changed the name of the parameter describing the abundance of neutral Sulfur from si to sui and added the missing parameters for singly and doubly ionized sulphur as suii and suiii respectively. I then modified every function that uses the parameter list to include the renamed/ added names for sulphur.
These are the new parameters with the names and values i used. The initial values and boundaries can/should be changed via threeML according to the individual use. I did not verify that the model gives correct numbers with this workaround, but it can be used without throwing an error.

sui:
delta: 0.001
desc: (see https://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/XspecModels.html)
free: true
initial value: 2.14
max: 1E6
min: '0.'
unit: 10^16
suii:
delta: 0.001
desc: (see https://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/XspecModels.html)
free: false
initial value: 0.0
max: 1E6
min: '0.'
unit: 10^16
suiii:
delta: 0.001
desc: (see https://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/XspecModels.html)
free: false
initial value: 0.0
max: 1E6
min: '0.'
unit: 10^16

@omodei
Copy link
Contributor

omodei commented Nov 17, 2021

Can this be closed?

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