gs2 ignores n0 and rhostar_single and just uses default aky instead

I have discovered the following issue in 8.0.3-RC compiled with make -j USE_NEW_DIAG= gs2 on one of the YPI servers (pinch). I have not tested other systems or versions.

With the following minimal input file, input.in:

&kt_grids_knobs
  grid_option = "single"
/

kt_grids_single_parameters
  rhostar_single = 0.0025
  n0 = 100
/

Running gs2 input.in and then ncdump -v ky input.out.nc results in ky = 0.4. Meanwhile, ncdump -v drhodpsi input.out.nc results in drhodpsi = 2.96201569410306. According to the documentation (old wiki), aky = n0 * drhodpsi * rhostar_single. Using the values above, ky should be 0.740503923525765000 so the four numbers are not consistent. Coincidentally, the default value for ky is 0.4 suggesting that gs2 is ignoring n0 and rhostar_single and just using the default value of ky. Obviously, this is not good if I want gs2 to calculate ky from the n0 and rhostar_single that I specify.

Comments (4)