XGC-209 case

Error in the initialization of the case XGC-209 is discovered. The case is resubmitted with the following additional changes:

fm_nsmooth = 0   ! number of applications of smoothing in fmcfm_call
fm_nsmooth_er = 4 ! number of smoothing for radial electric field in the core region
fm_smooth_er_rhob = 0.95 ! rho, where transition between different level of smoothing is applied
fm_smooth_er_coef = 50. ! coefficient that controls sharpness of transition between different level of smoothness for Er
fm_nsmooth_chi = 2 ! number of applications of diffusivity smoothing in fmcfm_call
fm_wexb   = 1.0d-6

Note that the flow shear stabilization coefficient is reduced from 1e-5 to 1e-6.

Averaging over three neighboring points is used in this case.


New two short simulations to test smoothing are submitted

New short XGC-0 simulations are submitted on FRANKLIN in the debug mode. The case XGC-210 is based on the case XGC-209. The case XGC-211 uses modified source code where smoothing for the radial electric field is done for three neighboring grid points instead of one (as in the case XGC-210).

IF ( fm_nsmooth_er > 0 ) THEN
  DO i=1, fm_nsmooth_er
    CALL smoothb0( z_e_radial, diag_flow_npsi, 3 )
  ENDDO
ENDIF

The radial electric field profiles as well as flow shear rates are significantly smoother comparing to the case where is smoothing is done over one neighboring grid point.

Results that use the same smoothing in the plasma core, but use averaging over one neighboring points are given below for comparison.