spacepy.radbelt.get_modelop_L
spacepy.radbelt.get_local_accel
calculate the diffusion in L at constant mu,K coordinates time units