Geopack.trace not providing vector of output values

tsssss / geopack

Python version of geopack and Tsyganenko models

MIT License

30 stars 12 forks source link

Geopack.trace not providing vector of output values #1

Closed w2ruf closed 5 years ago

w2ruf commented 5 years ago

Hi,

I am trying to do a field line trace and get the points for the entire field line. However, I am only getting a single set of coordinates even though the documentation for the function says it should also output a vector of coordinates. Could you please help?

Example Code: geopack.trace(1,0,.5, 1,10,1,2,'t89','igrf')

Output: (0.9921003911113675, 0.00048052726154043516, 0.1254269006978363)

Thank you!

tsssss commented 5 years ago

Hello, I just made trace to return the field lines as well. I haven't had the setup at hand to test. So could you please verify if the change works? BTW: Thanks for catching this! -Sheng

w2naf commented 5 years ago

Hi Sheng,

Thanks for your help on this! There was another small bug, but I fixed it in https://github.com/tsssss/geopack/pull/2, added a keyword for maxloop, and created a notebook demonstrating field line tracing.

Best, Nathaniel