The ExB flow shear rate has been implemented as simple derivative of the radial electric field with respect to minor radius:
New implementation is according to Ref. [T.S. Hahm and K.H. Burrell Phys. Plasmas 2, 1648 (1995)]:
where 
The FMCFM interface is modified and particle diffusivity is updated for this case with tbl_d_coeff_in;:
fm_ion_diff_norm%particle(i) = tbl_d_coeff_in * fm_norm_diff
The particle diffusivity in the core region is selected to be 
The fm_wexb coefficient has been increased from 1e-6 to 1e-4. This factor determines the effect of ExB flow shear suppression on the anomalous transport:
There is no dramatical particle losses as in the case of MMM95 model implemented using the FMCFM interface.
The time step in the case XGC-213 has been increased from 2e-3 to 2e-2. The plasma profiles looks more noisy.
Summary: I should try to increase time step by a factor of 2 rather than 10.
Simulation results for the case XGC-209 that uses smoothing procedure introduced on January 12. In addition, effect of fm_wexb parameter on the particle transport is studied. The results show noticeable smoothing for the radial electric field profiles in the plasma core and very little change in the plasma density dynamics.
Smoothing procedure adopted from the ASTRA code on Friday (see previous post) is tested for set of 200 random numbers. The plots are generated for the smoothness parameter 
In order to introduce more smoothing in the plasma core which is needed to compensate additional numerical noise in the core introduced by volume effect, 
where 
where R is a random number in the range from 0 to 1, 
The level of smoothness in the plasma core is controlled by 
The idea for smoothing of plasma profiles is based on routine implemented in the ASTRA code. The resulting function is based on finding of the minimum of the functional
with respect to U. The functional is design for the cylindrical case and the resulting function has zero derivative for x=0 and preserve boundary condition at x=1. The description of the method can be found at:
New simulation for the whole plasma profile with XGC-0 standard anomalous transport model is submitted on FRANKLIN. The FMCFM interface is not used. All other settings are the same.
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.
