Compute cutoff rigidities time series for the calculation of scaling factors according to Lifton et al. (2014) 'LSD' scaling scheme

Adapted from the matlab code in the supplementary material of Lifton et al. (2014)

lsdrc(lat, lon, time)

Arguments

lat

site latitude (deg, -90->90 southern hemisphere negative)

lon

site longitude (deg, -180->180 western hemisphere negative)

time

time before present (a, with respect to 2010)