README.md 1.54 KB
Newer Older
Matthias Steffen's avatar
Matthias Steffen committed
1
**Python 3 script for evaluating fitting fuctions FFI, FFII, and FFII**
Matthias Steffen's avatar
Matthias Steffen committed
2
3
4
5
6
7
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, A&A ???, ???)

The three functions are defined in ff_all.py
Matthias Steffen's avatar
Matthias Steffen committed
8
They are called in the example script test_ff123.py that produces the output file 'test_ff123_py.out'
Matthias Steffen's avatar
Matthias Steffen committed
9
10

With Teff = effective tempaerature, log g = surface gravity, Z = metallicity [Fe/H], qiso = lithium isotopic ratio 6Li/7Li,
Matthias Steffen's avatar
Matthias Steffen committed
11
the call to evaluate fitting function **FFI ** is:
Matthias Steffen's avatar
Matthias Steffen committed
12
13
14
15

`ff1_out=ff1(Teff,logg,Z,ALi_LTE,qiso)`

where ff1_out is a two-element array returing the (logarithmic) lithium abundance corrections
Matthias Steffen's avatar
Matthias Steffen committed
16

Matthias Steffen's avatar
Matthias Steffen committed
17
18
`ff1_out[0] = 1D_NLTE correction`
`ff1_out[1] = 3D_NLTE correction`
Matthias Steffen's avatar
Matthias Steffen committed
19

Matthias Steffen's avatar
Matthias Steffen committed
20
21
for a given LTE lithium abundance ALi_LTE

Matthias Steffen's avatar
Matthias Steffen committed
22
The call to evaluate fitting function **FFII** is:
Matthias Steffen's avatar
Matthias Steffen committed
23
24
25
26

`ff2_out=ff2(Teff,logg,Z,ALi,qiso)`

where ff2_out is a two-element array returing the (logarithmic) non-LTE equivalent widths (EW)
Matthias Steffen's avatar
Matthias Steffen committed
27
28
29
30

`ff2_out[0] = 1D_NLTE EW`
`ff2_out[1] = 3D_NLTE EW`

Matthias Steffen's avatar
Matthias Steffen committed
31
32
for a given lithium abundance ALi

Matthias Steffen's avatar
Matthias Steffen committed
33
The call to evaluate fitting function **FFIII** is:
Matthias Steffen's avatar
Matthias Steffen committed
34
35
36

`ff3_out=ff3(Teff,logg,Z,EW,qiso)`

Matthias Steffen's avatar
Matthias Steffen committed
37
where ff3_out is a two-element array returning the (logarithmic) 1D NLTE and 3D NLTE lithium abundances
Matthias Steffen's avatar
Matthias Steffen committed
38
39
40
41

`ff3_out[0] = ALI_1DNLTE`
`ff3_out[1] = ALi_3DNLTE`

Matthias Steffen's avatar
Matthias Steffen committed
42
43
for a given (logarithmic) equivalent width EW.

Matthias Steffen's avatar
Matthias Steffen committed
44
The IDL script test_ff123.pro, provided for the convenience of IDL users, produces the output file 
Matthias Steffen's avatar
Matthias Steffen committed
45
46
test_ff123_idl.out with idenitcal content as test_ff123_py.out.

Matthias Steffen's avatar
Matthias Steffen committed
47