madamow / pymoogi

python wrapper for MOOG. All plots are handled with matplotlib.
10 stars 4 forks source link

Phoenix file type support #8

Closed WLNORGREN closed 1 year ago

WLNORGREN commented 3 years ago

I am using pymoogi on the following file types, here is an example of what I am trying to do. I am using the abfind algorithm with the following file types, the atmosphere is in phoenix format. Here is the .par file I am using:

abfind standard_out 'b' summary_out 'bb' model_in 'GJ3313.phx' lines_in 'GJ3313.lin' atmosphere 1 molecules 1 lines 1 flux/int 0 plot 0 terminal xterm abundances 1 1 29 0.0 damping 2

The name of this file is abgj3313.par, and the model atmosphere and line files are GJ3313.phx and GJ3313.lin, respectively. When I run pymoogi abfind_example.par, there are no issues. However, when I run pymoogi abgj3313.par , I get the following error message:

Screenshot from 2021-04-02 10-59-39

In addition, I have modified Inmodel.f to support the phoenix model atmosphere type as follows: Inmodelf.txt

Is it possible to use the phoenix model atmosphere type with pymoogi, or have I configured something incorrectly?

madamow commented 3 years ago

Form the quick look on your error, it seems like the problem is with the format of your output files. Could you send me an example of your output files? I can work on adjusting pymoogi to your case.

WLNORGREN commented 3 years ago

Hello Ms. Adamow,

Do you mean b and bb?

Thank you,

Wolf Norgren

On Fri, Apr 2, 2021 at 12:01 PM Monika Adamow @.***> wrote:

Form the quick look on your error, it seems like the problem is with the format of your output files. Could you send me an example of your output files? I can work on adjusting pymoogi to your case.

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/madamow/pymoogi/issues/8#issuecomment-812617160, or unsubscribe https://github.com/notifications/unsubscribe-auth/ATP4Q7GCIDAD4DRVEHU2XPLTGXZ45ANCNFSM42JCLHBQ .

WLNORGREN commented 3 years ago

Hello Ms. Adamow,

Thank you for helping me with this software!

Both me and my professor greatly appreciate pymoogi!

Thank you,

Wolf Norgren

On Fri, Apr 2, 2021 at 5:03 PM W Norgren @.***> wrote:

Hello Ms. Adamow,

Do you mean b and bb?

Thank you,

Wolf Norgren

On Fri, Apr 2, 2021 at 12:01 PM Monika Adamow @.***> wrote:

Form the quick look on your error, it seems like the problem is with the format of your output files. Could you send me an example of your output files? I can work on adjusting pymoogi to your case.

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/madamow/pymoogi/issues/8#issuecomment-812617160, or unsubscribe https://github.com/notifications/unsubscribe-auth/ATP4Q7GCIDAD4DRVEHU2XPLTGXZ45ANCNFSM42JCLHBQ .

madamow commented 3 years ago

Yes, b and bb. I will start from there. It is Friday today, so I will probably contact you back early next week :) Thank you for using pymoogi!

WLNORGREN commented 3 years ago

Hello Ms. Adamow,

I have changed b and bb to eqw1 and eqw2. That seems to work!

Thank you for the help!

Is it OK if i contact you in the future about further pymoogi issues?

I am in an ongoing research project relating to stellar atmospheres.

Th

On Fri, Apr 2, 2021 at 12:07 PM Monika Adamow @.***> wrote:

Yes, b and bb. I will start from there. It is Friday today, so I will probably contact you back early next week :) Thank you for using pymoogi!

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/madamow/pymoogi/issues/8#issuecomment-812620899, or unsubscribe https://github.com/notifications/unsubscribe-auth/ATP4Q7E7G535SZEWEFH6XT3TGX2V5ANCNFSM42JCLHBQ .

WLNORGREN commented 3 years ago

Thank you,

Wolf Norgren

On Fri, Apr 2, 2021 at 12:12 PM W Norgren @.***> wrote:

Hello Ms. Adamow,

I have changed b and bb to eqw1 and eqw2. That seems to work!

Thank you for the help!

Is it OK if i contact you in the future about further pymoogi issues?

I am in an ongoing research project relating to stellar atmospheres.

Th

On Fri, Apr 2, 2021 at 12:07 PM Monika Adamow @.***> wrote:

Yes, b and bb. I will start from there. It is Friday today, so I will probably contact you back early next week :) Thank you for using pymoogi!

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/madamow/pymoogi/issues/8#issuecomment-812620899, or unsubscribe https://github.com/notifications/unsubscribe-auth/ATP4Q7E7G535SZEWEFH6XT3TGX2V5ANCNFSM42JCLHBQ .

madamow commented 3 years ago

I am happy that it worked. That is weird though that something depends on the name of the output files, so it is something I need to check. Please contact me about any issues with pymoogi, always happy to help!

WLNORGREN commented 3 years ago

Hello Ms. Adamow,

I had previously thought that changing the output names worked, but instead, it seems like I keep seeing the data for the example star, instead of the one I am trying to calculate the abundances for, when the output files are set to eqw1 and eqw2. What do you think I am doing wrong? I am still trying to use Phoenix model atmospheres.

Thank you,

Wolf Norgren

On Fri, Apr 2, 2021 at 12:16 PM Monika Adamow @.***> wrote:

I am happy that it worked. That is weird though that something depends on the name of the output files, so it is something I need to check. Please contact me about any issues with pymoogi, always happy to help!

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/madamow/pymoogi/issues/8#issuecomment-812624720, or unsubscribe https://github.com/notifications/unsubscribe-auth/ATP4Q7BGGDRCWUYTKCIMUKLTGX3XXANCNFSM42JCLHBQ .

madamow commented 3 years ago

Ok, I will work on that. Could you please send me an example file with phoenix model?

WLNORGREN commented 3 years ago

Sure, here it is.

Hello Ms. Adamow,

Sorry for not getting back sooner.

The first screenshot is what I see when I first log in.

The second is the data for carbon for abfind_example.par The third is the data for abGJ3313.par, a star I am trying to examine.

In both cases, I run pymoogi, and the output is identical.

Is there anything I haven't configured correctly?

I am using eqw1 and eqw2 for the output files in both cases.

On Tue, Apr 27, 2021 at 12:00 PM Monika Adamow @.***> wrote:

Ok, I will work on that. Could you please send me an example file with phoenix model?

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/madamow/pymoogi/issues/8#issuecomment-827763257, or unsubscribe https://github.com/notifications/unsubscribe-auth/ATP4Q7B2DKAIYYXKAKHKCWTTK3US3ANCNFSM42JCLHBQ .

WLNORGREN commented 3 years ago

Hello Ms. Adamow,

Here are the other two attachments.

Thank you,

Wolf Norgren

On Tue, Apr 27, 2021 at 2:24 PM W Norgren @.***> wrote:

Sure, here it is.

Hello Ms. Adamow,

Sorry for not getting back sooner.

The first screenshot is what I see when I first log in.

The second is the data for carbon for abfind_example.par The third is the data for abGJ3313.par, a star I am trying to examine.

In both cases, I run pymoogi, and the output is identical.

Is there anything I haven't configured correctly?

I am using eqw1 and eqw2 for the output files in both cases.

On Tue, Apr 27, 2021 at 12:00 PM Monika Adamow @.***> wrote:

Ok, I will work on that. Could you please send me an example file with phoenix model?

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/madamow/pymoogi/issues/8#issuecomment-827763257, or unsubscribe https://github.com/notifications/unsubscribe-auth/ATP4Q7B2DKAIYYXKAKHKCWTTK3US3ANCNFSM42JCLHBQ .

madamow commented 3 years ago

I do not see attached files. Could you send them directly to my email? madamow[at]icloud.com

madamow commented 3 years ago

My guess is the new output files are not created, old eqw1 and eqw2 are read in instead. But to recreate the error I need input files, at least the phoenix model file.

WLNORGREN commented 3 years ago

Hello Ms. Adamow,

Here are the input files I used.

Thank you,

Wolf Norgren

On Tue, Apr 27, 2021 at 3:23 PM Monika Adamow @.***> wrote:

My guess is the new output files are not created, old eqw1 and eqw2 are read in instead. But to recreate the error I need input files, at least the phoenix model file.

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/madamow/pymoogi/issues/8#issuecomment-827904670, or unsubscribe https://github.com/notifications/unsubscribe-auth/ATP4Q7EIAMO4PYAYHROPB53TK4MLDANCNFSM42JCLHBQ .