PrincetonUniversity / SPEC

The Stepped-Pressure Equilibrium Code, an advanced MRxMHD equilibrium solver.
https://princetonuniversity.github.io/SPEC/
GNU General Public License v3.0
25 stars 6 forks source link

Changeangle #54

Closed jloizu closed 6 years ago

jloizu commented 6 years ago

I have modified the routine global.h in order to make SPEC deal with cases in which an initial guess is given in the input file for the inner surfaces (Linitialize=0) and such that the poloidal angle corresponding to the provided geometric representation has opposite sign to the SPEC poloidal angle. Namely, the Fourier coefficients for R,Z are transformed in such a way that the poloidal angle has the "correct" sign. This was already being done in SPEC for the boundary but not for internal surfaces.

This modification has enabled to run SPEC for a CFQS equilibrium in which initial guess had this property.

Please check that the branch compiles and runs for some test cases of your choice.

Then we can merge this into master.