scalingfactor = stone2000(latitude, pressure, fsp)
Calculates the geographic scaling factor for cosmogenic-nuclide production with varying latitude and altitude according to the scheme in:
Stone, J., 2000. Air pressure and cosmogenic isotope production. Journal of Geophysical Research 106, pp. 23,753-23,759.
The input arguments are
latitude, the site latitude (decimal degrees),
pressure, the atmospheric pressure at the site (hPa), and
fsp the fraction of production due to neutron spallation (nondimensional). Accepts vector arguments. The argument
fsp, if omitted, defaults to 0.978, which is the correct value for BeThe corresponding value for Al is 0.974. Note that using
fsp = 0.844 (Be) and
fsp = 0.844 (Al) will closely reproduce the scaling factors given in Lal (1991) as long as the standard atmosphere is used. Note also that this function can be made to yield the scaling factor for spallation only by specifying
fsp = 1 and the scaling factor for production by muons only by specifying
fsp = 0.
The geographic production rate scaling factor is:
where is the fraction of production due to spallation (input argument
fsp), is the scaling factor for production by neutron spallation:
and is the scaling factor for production by muons:
Where the constants and depend on the latitude and are defined for certain index latitudes:
This table duplicates Table 1 in Stone (2000).
This function actually calculates and for index latitudes that bound the actual latitude of the site, and then determines the value at the site by linear interpolation. The exact method of interpolation used is relatively unimportant, but we have chosen linear interpolation to avoid the problem of under- and over-shooting between index latitudes due to the polynomial form of the Lal equations. Basically, linear interpolation guarantees that the scaling factor at the site is not outside the range of the scaling factors at the bounding index latitudes.