geopack
geopack copied to clipboard
Feature Request: Evaluate models at multiple points
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.