scalingfactor = desilets2006sp(h,Rc)
Calculates the geographic scaling factor for cosmogenic-nuclide production for particular cutoff rigidity and atmospheric pressure according to the scheme in:
Desilets D., Zreda M., Prabu T., 2006. Extended scaling factors for in situ cosmogenic nuclides: New measurements at low latitude. Earth and Planetary Science Letters, v. 246, pp. 265-276.
The input arguments are
h, atmospheric pressure (hPa), and
Rc, cutoff rigidity (GV). Accepts vector arguments.
Converts atmospheric pressure (hPa) to atmospheric depth (g cm) by .
Assigns cutoff rigidities below 2 GV a value of 2 GV.
Obtains the effective attenuation length in air via Equation (4) in the source paper:
where is the cutoff rigidity (GV) and the function is:
Obtains the altitude scaling factor :
Obtains the latitude scaling factor via Equation (6) in the source paper:
where is the cutoff rigidity, , and .
Finally, obtains the total scaling factor .
Accepts either scalars or vectors of equal sizes for all the input arguments. Returns either a scalar or a vector of the appropriate size.