It would be great to have the ability to evaluate the Tsyganenko models at multiple points using vectorized processing. For example, evaluating B for every point on a grid. Some scientific use cases:
Model comparison in the global perspective
Evaluating the models on a grid, and then plugging the result into other code that uses gridded output
This should be possible if operations are all done in a way that supports arrays in addition to individual points. Right now, I can iterate on every point in the grid, but it's very slow to do it this way because the Python interpreter is doing a lot of unnecessary work.
This is an excellent suggestion. I'll find some time this summer (but unfortunately no guarantee) to do a significant upgrade, including improving efficiency and adding functions as you and others have suggested.
It would be great to have the ability to evaluate the Tsyganenko models at multiple points using vectorized processing. For example, evaluating
B
for every point on a grid. Some scientific use cases:This should be possible if operations are all done in a way that supports arrays in addition to individual points. Right now, I can iterate on every point in the grid, but it's very slow to do it this way because the Python interpreter is doing a lot of unnecessary work.