Demonstrated correct use of ExtrapolateMomentsToMesh to compute plasma...

&params_nml

! Simulation name:

! ===============

params%fileDescriptor = '2021_03_09b',

params%fileDescriptor = '2021_03_10a',

! Magnetic field input data:

! =========================

! Define local variables:

INTEGER(i4) :: i, ix, frame

REAL(r8) :: vpar, Ma, Ep, alpha, xip, vper, v, a

REAL(r8), DIMENSION(mesh%NZmesh + 4) :: n

! Initialize mesh quantities:

mesh%n = 0.

mesh%nU = 0.

mesh%unU = 0.

mesh%P11 = 0.

mesh%P22 = 0.

mesh%nuE = 0.

n = 0.

!mesh%nU = 0.

!mesh%unU = 0.

!mesh%P11 = 0.

!mesh%P22 = 0.

!mesh%nuE = 0.

! Species "a" mass:

Ma = params%Ma

alpha = plasma%alpha

! Calculate moments and extrapolate to mesh points

!$OMP PARALLEL DO private(ix,a,Ep,v,xip,vpar,vper)

! mesh quantities need to be private to prevent race conditions 

! if they are declared shared

!$OMP PARALLEL DO PRIVATE(ix,a,Ep,v,xip,vpar,vper) REDUCTION(+:n)

DO i = 1,params%NC

	IF (plasma%f1(i) .EQ. 0 .AND. plasma%f2(i) .EQ. 0) THEN

		ix = plasma%m(i) + 2

		! Plasma density: [NSP m^-3]

		mesh%n(ix-1) = mesh%n(ix-1) + plasma%wL(i)*a;

		mesh%n(ix)   = mesh%n(ix)   + plasma%wC(i)*a;

		mesh%n(ix+1) = mesh%n(ix+1) + plasma%wR(i)*a;

		n(ix-1) = n(ix-1) + plasma%wL(i)*a;

		n(ix)   = n(ix)   + plasma%wC(i)*a;

		n(ix+1) = n(ix+1) + plasma%wR(i)*a;

	END IF

END DO

!$OMP END PARALLEL DO

! Apply magnetic compression:

mesh%n = mesh%n/params%Area0

mesh%n = n/params%Area0

! Apply scaling factor:

mesh%n = alpha*mesh%n/mesh%dzm

# Set number of threads:

# ===================================================

export OMP_NUM_THREADS=96

export OMP_NUM_THREADS=48

# Set processor binding for openMP:

# ===================================================