Secure the moments artificial dissipation w.r.t. the p paralleliz.

src/moments_eq_rhs_mod.F90

@@ -140,8 +140,9 @@ SUBROUTINE moments_eq_rhs_e
                 ! Nonlinear term
                 -hatB_NL(iz,eo) * Sepj(ip,ij,iky,ikx,iz)
 
-            IF(ip-4 .GT. 0) &
+              IF( (ip-4 .GT. 0) .AND. (num_procs_p .EQ. 1) ) &
               ! Numerical parallel velocity hyperdiffusion "+ dvpar4 g_a" see Pueschel 2010 (eq 33)
+              ! (not used often so not parallelized)
               moments_rhs_e(ip,ij,iky,ikx,iz,updatetlevel) = &
                 moments_rhs_e(ip,ij,iky,ikx,iz,updatetlevel) &
                   + mu_p * moments_e(ip-4,ij,iky,ikx,iz,updatetlevel)
@@ -290,8 +291,9 @@ SUBROUTINE moments_eq_rhs_i
                   ! Nonlinear term with a (gxx*gxy - gxy**2)^1/2 factor
                   -hatB_NL(iz,eo) * Sipj(ip,ij,iky,ikx,iz)
 
-                  IF(ip-4 .GT. 0) &
+                  IF( (ip-4 .GT. 0) .AND. (num_procs_p .EQ. 1) ) &
                     ! Numerical parallel velocity hyperdiffusion "+ dvpar4 g_a" see Pueschel 2010 (eq 33)
+                    ! (not used often so not parallelized)
                     moments_rhs_i(ip,ij,iky,ikx,iz,updatetlevel) = &
                       moments_rhs_i(ip,ij,iky,ikx,iz,updatetlevel) &
                         + mu_p * moments_i(ip-4,ij,iky,ikx,iz,updatetlevel)