If 'include_neoclassical_terms' = True in the 'neoclassical_input' namelist, the code will abort, since I found bugs in both the Miller and VMEC implementations of gds23, gds24, gds25 and gds26.
Specifically gds25 and gds26 were badly defined within VMEC and are never used within stella. In Miller, gds25 and gds26 are not defined, and gds23 and gds24 were lost during the geo_spline() routine, where they were overwritten by gds21 and gds22 instead.
Based on the amount of obvious mistakes and loose ends, the inclusion of neoclassical terms now has an abort statement, since it needs to be properly implemented first.
If 'include_neoclassical_terms' = True in the 'neoclassical_input' namelist, the code will abort, since I found bugs in both the Miller and VMEC implementations of gds23, gds24, gds25 and gds26.
Specifically gds25 and gds26 were badly defined within VMEC and are never used within stella. In Miller, gds25 and gds26 are not defined, and gds23 and gds24 were lost during the geo_spline() routine, where they were overwritten by gds21 and gds22 instead.
Based on the amount of obvious mistakes and loose ends, the inclusion of neoclassical terms now has an abort statement, since it needs to be properly implemented first.