madamow / pymoogi

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

Problem when using "blends" driver #26

Open flaulhabe opened 8 months ago

flaulhabe commented 8 months ago

Hello,

I'm not being able to successfully use the "blends" driver in pymoogi. It seems that the lines are not getting recognized and the output file don't give the elemental abundances. I've tried this same process with the same files in MOOG 2019 in another computer and it worked, and I still haven't figured out why this is not working in pymoogi.

I'm attaching the inp file used to run pymoogi, the input data and the output.

Thank you!

IC2714-28_La_hf_list.txt IC2714-28_atm.txt inp_la_hf.txt IC2714-28_zr_out2_La_ab.txt

madamow commented 8 months ago

Hi! I am on winter break now. I will take a look at this issue on Monday. It may be a bug in blends.

madamow commented 7 months ago

Can you post the correct output you expect?

flaulhabe commented 7 months ago

Hello, I'm attaching the expected output.

The input is expected to be a list of spectral lines according to the elemental hyperfine structure. In this case, the input I've attached earlier is a list of Lanthanum spectral lines. As result, I expect an output giving the Lanthanum abundance, in a similar way to the abfind output.

IC2714-28_zr_out2_La_ab.txt

madamow commented 7 months ago

This file has 999.990 and one NaN in abundance, so I assume this is an output that is not correct. You mention running MOOG on your input files on a different machine and getting the correct output. Could you post that file, please?

flaulhabe commented 7 months ago

Hey, Sorry, I've mistaken the files. I'm attaching the right one. In this case, MOOG provided a mean lanthanum abundance of 1.279.

IC2714-5_zr_out2_La_ab.txt

madamow commented 6 months ago

Here is an output I am getting. It looks correct. IC2714-28_zr_out2_La_ab.txt

If you see differences between outputs from different machines, this may indicate that Fortran is not correctly installed on one of them. I cannot recreate that error, though.