sdelarquier / tsyganenko Goto Github PK

I ran the following code:

import tsyganenko
import numpy as np

lats = np.concatenate((np.arange(-90,95,5), np.arange(90,-95,-5))
lons = no.concatenate((np.zeros(len(lats) / 2.0), 180.0 * np.ones(len(lats) / 2.0)))
rhos = 6371.0 * np.ones(len(lats))

trace = tsyganenko.tsygTrace(lats, lons, rhos)

trace.plot()

The image below is the result, and appears to yield a magnetopause on the dayside of Earth at approximately 2 Earth radii. This is quite obviously far too close – there must be a bug somewhere in this T96 implementation which is causing incorrect field line traces. Is it possible to diagnose and fix this issue?