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).