Rahul Gaur

https://iopscience.iop.org/article/10.1088/0029-5515/56/9/092006/pdf section 3 actually explains how to get J contours.

> I'd like to figure out what we need to fix so that calling eq.compute with any grid "just works". That should also fix the issues with plotting since all...

```python theta_unwrapped = np.unwrap(theta_values) phi_unwrapped = np.unwrap(phi_values) iota = (theta_unwrapped[-1] - theta_unwrapped[0]) / ( phi_unwrapped[-1] - phi_unwrapped[0] ) ``` Once you have your theta, zeta along a fieldline, the above...

Good idea. But I would be surprised if these folks are coupling flux tube geometry calculation with an actual gyrokinetic solver for speed. It is important to note that even...