hiddenSymmetries / booz_xform

Calculates Boozer coordinates for toroidal MHD equilibria, including stellarators and tokamaks.
https://hiddensymmetries.github.io/booz_xform/
BSD 2-Clause "Simplified" License
12 stars 3 forks source link

Phip fix #17

Closed ejpaul closed 4 months ago

ejpaul commented 4 months ago

This branch fixes a bug that was caused by #11, in which the poloidal flux was saved as the variable phip_b in the netcdf file. However, the toroidal flux derivative was saved with the same name in the STELLOPT version. The variable phip_b is now saved as the toroidal flux derivative, while the poloidal flux is saved as chi_b if read_wout is called with flux = True. The netcdf read error is handled in the case that chi_b is not written to the file. The pressure, pres, which is required for AE calculations, is also read from the wout file and saved in the netcdf file if flux = True.