Incosistency in toroidal flux calculated for diferent cocos/es
As reported by @malinak
loaded equilibrium with cocos = 3, 13 shows the different magnitude of flux function.
shot_no = 3694
time_step = 2.5
eq = compass.cudb(shot_no, time_step)
eq.to_geqdsk(new_fname, 129, 257)
geqdsk = omfit_eqdsk.OMFITgeqdsk(new_fname)
Goal of this issues:
- check the correctness of toroidal flux function calculation in plaque and during g-file export
- improve documentation to clarify function calculation
- add
poloidal_fluxfunction