Commit b5231bee authored by Matthias Steffen's avatar Matthias Steffen
Browse files

Update README.md

parent 4e1d161c
......@@ -4,11 +4,12 @@ developed in the paper
**Improving spectroscopic lithium abundances. Fitting functions for 3D non-LTE corrections in FGK stars of different metallicity**
by A. Mott et al. (2020, Astronomy & Astrophysics ???, ???)
The three functions are defined in ff_all.py
They are called in the example script test_ff123.py that produces the output file 'test_ff123_py.out' (see description below)
The three python functions are defined in ff_all.py, the IDL functions in ff1.pro, ff2.pro, ff3.pro
They are called in the example script test_ff123.py (test_ff123.pro) that produces the output file 'test_ff123_py.out' (`test_ff123_idl.out`)
(see description of the contents of the output file below)
With Teff = effective temperature, log g = surface gravity, Z = metallicity [Fe/H], qiso = lithium isotopic ratio 6Li/7Li,
the call to evaluate fitting function **FFI ** is:
the call to evaluate fitting function **FFI** is:
`ff1_out=ff1(Teff,logg,Z,ALi_LTE,qiso)`
......@@ -23,12 +24,12 @@ The call to evaluate fitting function **FFII** is:
`ff2_out=ff2(Teff,logg,Z,ALi,qiso)`
where ff2_out is a two-element array returning the (logarithmic) non-LTE equivalent widths (EW)
where ff2_out is a two-element array returning the (logarithmic) non-LTE equivalent widths, log_10(EW/mA)
`ff2_out[0] = 1D_NLTE EW`
`ff2_out[1] = 3D_NLTE EW`
`ff2_out[0] = 1D_NLTE log_10(EW/mA)`
`ff2_out[1] = 3D_NLTE log_10(EW/mA)`
for a given lithium abundance ALi
for a given lithium (logarithmic) abundance ALi, where A(Li)=log_10(N(Li)/N(H)) + 12
The call to evaluate fitting function **FFIII** is:
......@@ -36,13 +37,12 @@ The call to evaluate fitting function **FFIII** is:
where ff3_out is a two-element array returning the (logarithmic) 1D NLTE and 3D NLTE lithium abundances
`ff3_out[0] = ALI_1DNLTE`
`ff3_out[0] = ALi_1DNLTE`
`ff3_out[1] = ALi_3DNLTE`
for a given (logarithmic) equivalent width EW.
for a given (logarithmic) equivalent width log_10(EW/mA).
The IDL script test_ff123.pro, provided for the convenience of IDL users, produces the output file
test_ff123_idl.out with identical content as test_ff123_py.out.
The output files produced by `test_ff123.py` and `test_ff123.pro` are identical.
**Contents of the output file**
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment