Dario Panici

- Should add documentation to the ensure positive jacobian saying that it will assume your inputs for current/Psi are what you physically want to preserve

> * check for stellarator if not negating current keeps iota profile magnitude the same This seems to work fine, at least going from VMEC input to DESC input, I...

However, the magnetic field poloidal direction seems to be flipped > * Should add documentation to the ensure nested jacobian saying that it will assume your inputs for current/Psi are...

you can just do eq.compute(["nu_mn","|B|_mn"])

```python from desc.examples import get for name in ["precise_QA","precise_QH","W7-X","ESTELL","NCSX","DSHAPE_current"]: eq = get(name) from desc.grid import * grid_conc = ConcentricGrid(L=eq.L_grid, M=eq.M_grid, N=eq.N_grid, NFP=eq.NFP, axis=False,sym=True) rho=grid_conc.nodes[grid_conc.unique_rho_idx,0] grid = LinearGrid(rho =rho, M=eq.M_grid, N=eq.N_grid,...

Try increasing innermost flux surface in ConcentricGrid poits to 3 instead of 1 and compare above plots again

For the mgrid, I did a cursory check to make sure that the dimension names and variable names that we use [match the mgrid example shared by Sam](https://github.com/PlasmaControl/DESC/issues/1550#issuecomment-2626532474), and it...

The below test works if `LambdaGauge` is disabled (on #1036 ). The equilibrium is pretty high aspect ratio though, so I should test this with a higher aspect ratio and...

Have `LambdaGauge` fix to whatever the current average is for each surface

Related to #1138, as I think right now NaN exist in plotted quantities (for example on axis) but the blanket suppression would prevent the warning from coming through