spacepy.empiricals.getLmax

spacepy.empiricals.getLmax(ticks, model='JKemp', dbase='QDhourly')[source]

calculate a simple empirical model for Lmax - last closed drift-shell

Uses the parametrized Lmax from: Koller and Morley (2010) ‘Magnetopause shadowing effects for radiation belt models during high-speed solar wind streams’ American Geophysical Union, Fall Meeting 2010, abstract #SM13A-1787

Parameters
ticksspacepy.time.Ticktock

Ticktock object of desired times

modelstring, optional

‘JKemp’ (default - empirical model of J. Koller)

Returns
outnp.ndarray

Lmax - L* of last closed drift shell

Examples

>>> from spacepy.empiricals import getLmax
>>> import spacepy.time as st
>>> import datetime
>>> ticks = st.tickrange(datetime.datetime(2000, 1, 1), datetime.datetime(2000, 1, 3), deltadays=1)
array([ 7.4928412,  8.3585632,  8.6463423])