tsssss/geopack

Geopack.trace not providing vector of output values

Closed this issue · 2 comments

w2ruf commented

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!

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

Hi Sheng,

Thanks for your help on this! There was another small bug, but I fixed it in #2, added a keyword for maxloop, and created a notebook demonstrating field line tracing.

Best,
Nathaniel