make_al_be_consts

This function has no arguments or outputs. It creates a structure called al_be_consts that contains all the constants relevant to the exposure age and erosion rate calculations, and saves it as a .mat file. The fields are:

al_be_consts.version	Version number for this function	string
al_be_consts.prepdate	Date and time the constants structure was created	1x6 double array
al_be_consts.l10	$^{10}$ Be decay constant (yr $^{-1}$ )	double
al_be_consts.l26	$^{26}$ Al decay constant (yr $^{-1}$ )	double
al_be_consts.Lsp	Effective attenuation length for production by spallation (g $\cdot$ cm $^{-2}$ )	double
al_be_consts.Natoms10	number density of O atoms in quartz (atoms $\cdot$ g $^{-1}$ )	double
al_be_consts.Natoms26	number density of Si atoms in quartz (atoms $\cdot$ g $^{-1}$ )	double
al_be_consts.k_neg10	summary yield for $^{10}$ Be production by negative muon capture in quartz (atoms $\cdot$ (stopped $\mu_{-}$ ) $^{-1}$ )	double
al_be_consts.delk_neg10	uncertainty in above (atoms $\cdot$ (stopped $\mu_{-}$ ) $^{-1}$ )	double
al_be_consts.k_neg26	summary yield for $^{26}$ Al production by negative muon capture in quartz (atoms $\cdot$ (stopped $\mu_{-}$ ) $^{-1}$ )	double
al_be_consts.delk_neg26	uncertainty in above (atoms $\cdot$ (stopped $\mu_{-}$ ) $^{-1}$ )	double
al_be_consts.sigma190_10	measured cross-section at 190 GeV for $^{10}$ Be production from O by fast muon reactions (cm $^{-2}$ )	double
al_be_consts.delsigma190_10	uncertainty in above (cm $^{-2}$ )	double
al_be_consts.sigma190_26	measured cross-section at 190 GeV for $^{26}$ Al production from Si by fast muon reactions (cm $^{-2}$ )	double
al_be_consts.delsigma190_26	uncertainty in above (cm $^{-2}$ )	double

al_be_consts.Fsp10	fraction of reference $^{10}$ Be production rate attributable to spallation according to Stone (2000) (nondimensional)	double
al_be_consts.Fsp26	fraction of reference $^{26}$ Al production rate attributable to spallation according to Stone (2000) (nondimensional)	double
al_be_consts.P10_ref_St	reference $^{10}$ Be production rate due to spallation for Lal(1991)/Stone(2000) scaling scheme (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.delP10_ref_St	uncertainty in above (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.P26_ref_St	reference $^{26}$ Al production rate due to spallation for Lal(1991)/Stone(2000) scaling scheme (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.delP26_ref _St	uncertainty in above (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.P10_ref_De	reference $^{10}$ Be production rate due to spallation for Desilets et al. (2006) scaling scheme (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.delP10_ref_De	uncertainty in above (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.P26_ref_De	reference $^{26}$ Al production rate due to spallation for Desilets et al. (2006) scaling scheme (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.delP26_ref_De	uncertainty in above (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.P10_ref_Du	reference $^{10}$ Be production rate due to spallation for Dunai (2001) scaling scheme (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.delP10_ref_Du	uncertainty in above (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.P26_ref_Du	reference $^{26}$ Al production rate due to spallation for Dunai (2001) scaling scheme (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.delP26_ref_Du	uncertainty in above (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.P10_ref_Li	reference $^{10}$ Be production rate due to spallation for Lifton et al. (2005) scaling scheme (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.delP10_ref_Li	uncertainty in above (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.P26_ref_Li	reference $^{26}$ Al production rate due to spallation for Lifton et al. (2005) scaling scheme (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.delP26_ref_Li	uncertainty in above (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.P10_ref_Lm	reference $^{10}$ Be production rate due to spallation for time-dependent adaptation of the Lal(1991)/Stone(2000) scaling scheme (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.delP10_ref_Lm	uncertainty in above (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.P26_ref_Lm	reference $^{26}$ Al production rate due to spallation for time-dependent adaptation of the Lal(1991)/Stone(2000) scaling scheme (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double
al_be_consts.delP26_ref_Lm	uncertainty in above (atoms $\cdot$ g $^{-1} \cdot$ yr $^{-1}$ )	double

al_be_consts.M	Ratio of past to present magnetic field intensity. Values from 7500-11500 yr are taken from Yang et al.; values from 12000 to 800000 are taken from SINT800; final two values (corresponding to times of 801,000 yr and Inf) are the average of the SINT800 record. (nondimensional)	1x796 double array
al_be_consts.t_M	Time vector to match above. Times are [7500:1000:11500 12000:1000:801000 Inf] (yr)	1x796 double array
al_be_consts.TTRc	3-d data block containing global maps of cutoff rigidities derived from the Korte and Constable magnetic field model by trajectory tracing (Lifton and others, in prep). Used in the De and Li scaling schemes. Dimensions are latitude, longitude, time - see below for step size in these. (GV)	37x25x15 double array
al_be_consts.IHRc	3-d data block containing global maps of cutoff rigidities derived from the Korte and Constable magnetic field model by applying Equation 2 of Dunai (2001) to inclination and horizontal field intensity values. These data are also due to Nat Lifton. Used in the Du scaling scheme. Dimensions are latitude, longitude, time - see below for step size in these. (GV)	37x25x15 double array
al_be_consts.lat_Rc	Latitude vector for interpolating cutoff rigidity data blocks. 5°spacing. (DD)	1x37 double array
al_be_consts.lon_Rc	Longitude vector for interpolating cutoff rigidity data blocks. 15°spacing. (DD)	1x25 double array
al_be_consts.t_Rc	Time vector for interpolating cutoff rigidity data blocks. Times are [0:500:6500 6900]. (yr)	1x15 double array
al_be_consts.MM0_KCL	Ratio of past to present magnetic field intensity at past times obtained by approximating the Korte and Constable magnetic field reconstruction with a geocentric dipole. Used in the Lm scaling scheme. Times match the al_be_consts.t_Rc time vector. (nondimensional)	15x1 double array
al_be_consts.lat_pp_KCL	Latitude of north magnetic pole position at past times obtained by approximating the Korte and Constable magnetic field reconstruction with a geocentric dipole. Used in the Lm scaling scheme. Times match the al_be_consts.t_Rc time vector. (nondimensional)	15x1 double array
al_be_consts.lon_pp_KCL	Longitude of north magnetic pole position at past times obtained by approximating the Korte and Constable magnetic field reconstruction with a geocentric dipole. Used in the Lm scaling scheme. Times match the al_be_consts.t_Rc time vector. (nondimensional)	15x1 double array
al_be_consts.S	Solar variability parameter used in Li scaling scheme. Resampled from data given by Lifton et al. (2005) to [0:500:6500 6900 7500:1000:10500] to match time steps in other magnetic field records. (nondimensional)	1x19 double array
al_be_consts.SInf	Long-term average solar variability parameter for use in Li scaling scheme before 10500 yr BP. See Lifton et al. (2005) for details. (nondimensional)	double

make_al_be_consts_v2.m