Cosmogenic nuclide calculation code

This appendix contains the functions used for cosmogenic nuclide exposure age calculations, erosion
simulations, and interquartile range calculations. The functions are identical to the functions in
Appendix 5 in SKB report TR-19-07 except that uplmode in expage_sealevel.m and glacialE_sealevel.m
has been changed to 'hunneberg'.

Files included:

 expage_sealevel.m     Function for calculating simple exposure age and deviation from expected
                       simple exposure age taking shielding during glacial isostatic uplift through
                       the water column into account.
 
 glacialE_sealevel.m   Function for simulation of glacial erosion taking shielding during glacial
                       isostatic uplift through the water column into account (based on glacialE.m
                       from expage-201902).
 
 depthcalc.m           Function for calculating depth profile exposure ages (derived from
                       expage-201902).
                       
Associated files used by expage_sealevel.m, glacialE_sealevel.m, and depthcalc.m:
 antatm.m
 ERA40.mat
 ERA40atm.m
 glacialE_LR04.txt
 interpolate.m
 LSD_fix.m
 LSDspal.m
 make_consts_expage.m
 Muons.m
 Neutrons.m
 PMag_Sep12.mat
 P_mu_expage.m
 Protons.m
 rawattenuationlength.m
 Reference.mat
 SKB_RSL.txt
 stone2000.m
 XSectsReedyAll.mat
 
 range_percentile.m    Function for calculating summed percentiles from a number of minimum-maximum
                       ranges, assuming uniform probability distribution between the individual
                       minimum and maximum values and equal weight of each individual minimum-
                       maximum range.
                       
All functions have been used in Octave 5.1. Most of the functions are expected to work with MATLAB
but some difference has been noted for the function glacialE_sealevel.m which does not work properly
when used in MATLAB.
