More options for Bharm calculation

[zhucaoxiang]

@logan-nc @ysm1887

The new feature of bharm_jsurf is now in the develop branch and ready for test.

To use it, just specify the bharm_jsurf = 0 (or 1, 2) in the input namelist. Here are the detalis of each option

Additional normalizations are added to be aligned with GPEC results.

• bharm_jsurf=0, GPEC = 2*FOCUS
• bharm_jsurf =1, GPEC = 2*(4*pi^2)*FOCUS
• bharm_jsurf=2, GPEC = 2*(2*pi)*FOCUS

Here is a simple test against GPEC using KSTAR RMP coils.

Zoomed in

Most of the modes are consistent. A discrepancy was observed in phasing angles for high m modes. One guess is that FOCUS is using a Fourier representation for coils, while GPEC might not (need confirmation).

[logan-nc]

Awesome! Thanks @zhucaoxiang!
@ysm1887 please let us know once you have put the necessary 2's and pi's into the OMFIT interface between GPEC and FOCUS. Once that is done, I'll close this issue.

[zhucaoxiang]

@logan-nc @ysm1887 I have put the additional normalizations (2’s, pi’s) into the FOCUS source code. The figure I showed above is the raw output from FOCUS.

[logan-nc]

Ok @zhucaoxiang, so we can use the GPEC "dominant mode" that has bharm_jsurf=2 as input to FOCUS without any additional normalization?

[logan-nc]

@zhucaoxiang the namelist I start from seems to contradict the documentation in this thread:

FOCUS/examples/d3d_RMP/d3d.input

Line 31 in 0505fd2

bharm_jsurf = 0 ! 0: no weightes; 1: weighted with area square; 2: weighted with area

Should we change all the namelists to read

! 0: no weightes;  1: weighted with area; 2: weighted with sqrt area

?