Simon Jeffery's Software Store
TAP: cosmic abundances.
Quick Reference
| SUBROUTINE TAP_ABINI ( TAP_ABUND,
STATUS ) |
Initialise the elemental abundances.
|
| FUNCTION TAP_ABMET ( ZELEM
) |
Return the relative meteoritic abundance for a given element.
|
| FUNCTION TAP_ABSOL ( ZELEM
) |
Return the relative solar photospheric abundance for a given element.
|
| FUNCTION TAP_ABCOS ( ZELEM
) |
Return the relative cosmic abundance for a given element.
|
| FUNCTION TAP_ABNUC ( ZELEM,
AELEM ) |
Return the solar system abundance for a given nuclide.
|
Long Reference
SUBROUTINE TAP_ABINI ( TAP_ABUND, STATUS )
Purpose:
Initialise the elemental abundances.
Invocation:
CALL TAP_ABINI ( TAP_ABUND, STATUS )
Description:
The elemental abundances used by the Theoretical Astrophysics
library routines are initialised. This is done by providing a
DOUBLE PRECISION FUNCTION TAP_ABUND as an argument, which is then
used to initialise the internal array of abundances. This allows
the user to use their own abundance data, if required.
In order to use TAP_ABINI, the routine TAP_ABUND must be declared
EXTERNAL in the calling routine.
Arguments:
TAP_ABUND = REAL (Given)
The abundance initialisation function. The FUNCTION TAP_ABUND
has one argument of type INTEGER which is the atomic number, Z,
of a given element. The FUNCTION will then return a value for
the logarithm abundance relative to hydrogen, where the
logarithm hydrogen abundance is 12.00. Where no abundance data
are available, or the atomic number is out of range, TAP_ABUND
returns the value TAPC__DNULL. The FUNCTIONS TAP_ABMET and
TAP_ABSOL may be used as templates for user code.
STATUS = INTEGER (Given and Returned)
The global status.
DOUBLE PRECISION FUNCTION TAP_ABCOS ( ZELEM )
Purpose:
Return the relative cosmic abundance for a given element.
Invocation:
RESULT = TAP_ABCOS( ZELEM )
Description:
Given the atomic number of the element, the recommended value of
the cosmic element abundance relative to hydrogen is returned.
The abundances are based primarily on the compilation of solar
and meteorite data by Anders & Grevesse (1989), supplemented in a
few cases by more recent observations.
In general, the solar and meteorite data are in good agreement.
Some exceptions (e.g. Lithium) are obvious, whilst others (e.g.
Indium and Tungsten) are less so. The solar data are returned
where possible, meteoritic data are returned where there are no
solar data.
All solar and meteoritic abundances are given in logarithmic form
relative to the solar hydrogen abundance, which is defined to be
12.00. Cases where no data are available for a given element
from either source are indicated by a returned value of
TAP__DNULL.
Arguments:
ZELEM = INTEGER (Given)
The atomic number (Z) for which the abundance is required.
Returned Value:
TAP_ABSOL = DOUBLE PRECISION
The cosmic abundance (logarithm with H = 12.00).
Accuracy:
See references for full details.
Solar abundances: errors lie between +/- 0.03 and +/- 0.10 dex.
Exceptions: F, Cl, K, Ge, Rb, Rh, Ag, Cd, In, Sn, Sb, Ce,
Tb, Dy, Ho, Tm, Yb, Lu, W, Pt, Au, Tl.
Meteor abundances: errors lie between +/- 0.01 and +/- 0.05 dex.
Exceptions: F, Cl, Br, Kr, Sb, I, Xe, Au.
References:
Anders E. & Grevesse N. 1989. Geochim. Cosmochim. Acta 53, 197.
Holden N.E., Martin R.L., & Barnes I.L. 1984. Pure Appl. Chem. 56,
675.
Solar Fe:
Holweger H., Heise C. & Kock M. 1990. Astron.Astrophys. 232,
510.
Solar C:
Stuerenberg S. & Holweger H. 1990. Astron.Astrophys. 237, 125.
DOUBLE PRECISION FUNCTION TAP_ABSOL ( ZELEM )
Purpose:
Return the relative solar photospheric abundance for a given
element.
Invocation:
RESULT = TAP_ABSOL( ZELEM )
Description:
Given the atomic number of the element, the recommended value
of the solar photospheric element abundance relative to hydrogen
is returned. The abundances are based primarily on the
compilation of solar and meteorite data by Anders & Grevesse
(1989), supplemented in a few cases by more recent observations.
In general, the solar and meteorite data are in good agreement.
Some exceptions (e.g. Lithium) are obvious, whilst others (e.g.
Indium and Tungsten) are less so. The solar data will find
applications in stellar theory, whilst it is expected that the
meteorite data should be more useful for ISM studies.
All solar and meteoritic abundances are given in logarithmic form
relative to the solar hydrogen abundance, which is defined to be
12.00. Cases where no data are available for a given element
are indicated by a returned value of TAP__DNULL.
Arguments:
ZELEM = INTEGER (Given)
The atomic number (Z) for which the abundance is required.
Returned Value:
TAP_ABSOL = DOUBLE PRECISION
The solar photospheric abundance (logarithm with H = 12.00).
Accuracy:
See references for full details.
Solar abundances: errors lie between +/- 0.03 and +/- 0.10 dex.
Exceptions: F, Cl, K, Ge, Rb, Rh, Ag, Cd, In, Sn, Sb, Ce,
Tb, Dy, Ho, Tm, Yb, Lu, W, Pt, Au, Tl.
Meteor abundances: errors lie between +/- 0.01 and +/- 0.05 dex.
Exceptions: F, Cl, Br, Kr, Sb, I, Xe, Au.
References:
Anders E. & Grevesse N. 1989. Geochim. Cosmochim. Acta 53, 197.
Holden N.E., Martin R.L., & Barnes I.L. 1984. Pure Appl. Chem. 56,
675.
Solar Fe:
Holweger H., Heise C. & Kock M. 1990. Astron.Astrophys. 232,
510.
Solar C:
Stuerenberg S. & Holweger H. 1990. Astron.Astrophys. 237, 125.
DOUBLE PRECISION FUNCTION TAP_ABMET ( ZELEM )
Purpose:
Return the relative meteoritic abundance for a given element.
Invocation:
RESULT = TAP_ABMET ( ZELEM )
Description:
Given the atomic number of the element, the recommended value for
the meteoritic element abundance relative to hydrogen is
returned. The abundances are based primarily on the compilation
of solar and meteorite data by Anders & Grevesse (1989),
supplemented in a few cases by more recent observations.
In general, the solar and meteorite data are in good agreement.
Some exceptions (e.g. Lithium) are obvious, whilst others (e.g.
Indium and Tungsten) are less so. The solar data will find
applications in stellar theory, whilst it is expected that the
meteorite data should be more useful for ISM studies.
All solar and meteoritic abundances are given in logarithmic form
relative to the solar hydrogen abundance, which is defined to be
12.00. Cases where no data are available for a given element
are indicated by a returned value of TAP__DNULL.
Arguments:
ZELEM = INTEGER (Given)
The atomic number (Z) for which the abundance is required
Returned Value:
TAP_ABMET = DOUBLE PRECISION
The solar system meteorite abundance (logarithm with H =
12.00).
Accuracy:
See references for full details.
Solar abundances: errors lie between +/- 0.03 and +/- 0.10 dex.
Exceptions: F, Cl, K, Ge, Rb, Rh, Ag, Cd, In, Sn, Sb, Ce,
Tb, Dy, Ho, Tm, Yb, Lu, W, Pt, Au, Tl.
Meteor abundances: errors lie between +/- 0.01 and +/- 0.05 dex.
Exceptions: F, Cl, Br, Kr, Sb, I, Xe, Au.
References:
Anders E. & Grevesse N. 1989. Geochim. Cosmochim. Acta 53, 197.
Holden N.E., Martin R.L., & Barnes I.L. 1984. Pure Appl. Chem. 56,
675.
Solar Fe:
Holweger H., Heise C. & Kock M. 1990. Astron.Astrophys. 232,
510.
Solar C:
Stuerenberg S. & Holweger H. 1990. Astron.Astrophys. 237, 125.
DOUBLE PRECISION FUNCTION TAP_ABNUC ( ZELEM, AELEM )
Purpose:
Return the solar system abundance for a given nuclide.
Invocation:
RESULT = TAP_ABNUC ( ZELEM, AELEM )
Description:
The recommended values for solar system nuclide abundance,
relative to unit element abundance, are returned. In general, the
terrestrial representative isotopic compositions recommended
by IUPAC (Holden et al. 1984) are given, as cited by Anders and
Grevesse (1989).
The IUPAC compositions represent "the chemicals and/or materials
most commonly encountered in the laboratory" and NOT NECESSARILY
"the most abundant natural material". For light elements (Li, B,
C, O) and elements with a small radiogenic component (Sr, Nd, Hf,
Os) the differences are small, see Holden et al. (1984), and
references cited therein, for details. For elements that are
strongly depleted in the Earth (H, N, noble gases) or are largely
radiogenic (Pb), solar system values have been used.
Algorithm:
The nuclide abundances are stored and returned as percentages for
each element.
The nuclide abundances are stored in a REAL array with an INTEGER
pointer array indicating the end of the previous element's data.
The data in this array are ordered sequentially by element and
nuclide for each element.
Arguments:
ZELEM = INTEGER (Given)
The atomic number (Z) for which the abundance is required.
AELEM = INTEGER (Given)
The atomic (or nuclide) mass number (A) for which the nuclide
abundance is required.
Returned Value:
TAP_ABNUC = DOUBLE PRECISION
The percentage relative nuclide abundance.
References:
Anders E. & Grevesse N. 1989. Geochim. Cosmochim. Acta 53, 197.
Holden N.E., Martin R.L., & Barnes I.L. 1984. Pure Appl. Chem. 56, 675