3
X_                 @   s   d Z dddZdS )a  
************************************ 
*      Solar number fraction       *
************************************
This code is part of the LIFELINE program.

Copyright (C) 2020-2021 University of Liege (Belgium)
Enmanuelle Mossoux (STAR Institute)
This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public
License as published by the Free Software Foundation, either version 3 of the License, or any later version.
This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied
warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
details.
You should have received a copy of the GNU General Public License along with this program.
If not, see <http://www.gnu.org/licenses/>.

Run: 
  from sun_abund import abund
  number_fraction=abund(ref="name")

Abundances taken from: https://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/node117.html#abund

# Parametres du modele 
# ====================
ref - String. The reference for the abundances among:
	wilm - Wilms, Allen & McCray (2000, ApJ 542, 914)
	angr - Anders E. & Grevesse N. (1989, Geochimica et Cosmochimica Acta 53, 197)
	aspl - Asplund M., Grevesse N., Sauval A.J. & Scott P. (2009, ARAA, 47, 481)
	feld - Feldman U.(1992, Physica Scripta 46, 202)
	aneb - Anders E. & Ebihara (1982, Geochimica et Cosmochimica Acta 46, 2363)
	grsa - Grevesse, N. & Sauval, A.J. (1998, Space Science Reviews 85, 161)
	lodd - Lodders, K (2003, ApJ 591, 1220)

Nc             C   s  | d krd} | dkrd}d}d}d}d}d}d}d}d}	d}
d	}d
}d}d}d}d}d}d}d}d}d}d}d}d}d}d}d}d}d}d}n,| dkrd}d}d}d}d}d}d}d}d}	d }
d}d!}d"}d#}d$}d%}d&}d'}d}d(}d)}d*}d+}d,}d-}d.}d}d/}d0}d1}n| d2krd}d3}d4}d5}d6}d7}d8}d}d}	d9}
d:}d;}d<}d=}d>}d?}d&}d@}dA}dB}dC}dD}dE}dF}dG}dH}d*}dI}dJ}d}n | dKkr"d}d}dL}dM}dN}dO}dP}d}d}	dQ}
d}d!}d"}d#}d$}d%}d&}dR}d}d(}dS}dT}d+}d,}d-}d=}d}d/}d0}d1}n| dUkrd}dV}dW}dX}dY}dZ}d[}d\}d]}	d^}
d_}d`}da}db}dc}dd}de}df}dg}dh}di}dj}dk}dl}dm}dn}do}dp}dq}dr}n| dskr,d}d3}dL}dM}dN}dt}du}dv}d}	dw}
d}d!}d"}d#}d$}dx}d&}d@}d}d(}dS}dT}d+}d,}d-}dH}d}d/}d0}d1}n| dykrd}dz}d{}d|}d}}d~}d8}d}d}	d}
d}d#}d}d}d}d}d}d}d}dB}d}d}d+}d}d&}d}d}dI}d}d}nt d g S |||||||||	|
||||||||||||||||||||gS )Nwilmg      ?g ?    giUMu/?gZ1?gXj@?g.?g:ܳS>gݺQ>g8q>g[mD>gs;N>gpP >gIhG|>gƼq>g3_>g17GWq>g'G>g"s2_d>g/D4>g(JuUv>gz]ʲ>angrg?8#:=g&=gaզY=g !7?g_>Y1\?gVbK?gu}c>gFn1 ?g
F?g>gI?g`Ⱥ>gf녧>g4>gֺ]s>gpmUм5>gብ>g9z>g:0yE>gH g-h>g\+-eq>g8xpc?g 
ݽ>gڧeQ>g|V^e>asplgӼɵ?gG@
=gAfc=gm6>g*31?g&(?gNN?gy;>1>gAH?gxz)e>gf녧 ?gk9?>gÔ/>g}/>g |>g;^>g͸_>9>gVxw>gWlFB>gLLJSS>gzb>g/tG ?gVMٻ>g创P>feldgP^+=g1@=g"9e=gRT:?g-C6?ghY ?gmby>>gM\m>g߈/|>anebgh o?ge">g[T=g7o'N>g=\r)=?gT?g_vj.7H?g创`>g$D"?gӋ5>g:`?g*#&,>g7_6K?g2iĵ>ghM=n>g
m>g8>g- >g?̔>gqmM>gR"w>g\1G>ge	w֊>g_eG|>gQWZ?g{.3v>gJ]>gGC1KT>g~mlh>grsag&:,B5?gkNu?g30&F?giUMu?gŸ́p>loddgׁsF?gX$R >g*A=g+I>gXj0?gr{^>g	yNl?g~>gv(>gZ1?gg+eS>g!>@>g9xm>gtR#>gB$gP>g>g[f>gmO>gbu>g`$ǊS>gl-f>zThis reference is undefined)print)refnhZnheZnliZnbenbncnnnonfZnneZnnaZnmgZnalZnsinpnsZnclZnarnkZncaZnscZntinvZncrZnmnZnfeZncoZnniZncuZnzn r   2/local/network/htdocs/Lifeline/codes3/sun_abund.pyabund$   s    





r   )N)__doc__r   r   r   r   r   <module>#   s   