From 6939fa973beec0c38b85f6c1ef3bdb6a3e172a6c Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Tue, 15 Mar 2022 08:27:00 +0100
Subject: [PATCH] back to 4th HD
---
src/moments_eq_rhs.F90 | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/moments_eq_rhs.F90 b/src/moments_eq_rhs.F90
index 4e5ce055..a7d1a71b 100644
--- a/src/moments_eq_rhs.F90
+++ b/src/moments_eq_rhs.F90
@@ -121,7 +121,7 @@ SUBROUTINE moments_eq_rhs_e
! Electrostatic background gradients
- i_ky * K_E * moments_e(ip,ij,ikx,iky,iz,updatetlevel) &
! Numerical hyperdiffusion (totally artificial, for stability purpose)
- - (mu_x*kx**2 + mu_y*ky**2)*moments_e(ip,ij,ikx,iky,iz,updatetlevel) &
+ - (mu_x*kx**4 + mu_y*ky**4)*moments_e(ip,ij,ikx,iky,iz,updatetlevel) &
! Collision term
+ TColl_e(ip,ij,ikx,iky,iz)
@@ -258,7 +258,7 @@ SUBROUTINE moments_eq_rhs_i
! Electrostatic background gradients
- i_ky * K_E * moments_i(ip,ij,ikx,iky,iz,updatetlevel) &
! Numerical hyperdiffusion (totally artificial, for stability purpose)
- - (mu_x*kx**2 + mu_y*ky**2)*moments_i(ip,ij,ikx,iky,iz,updatetlevel) &
+ - (mu_x*kx**4 + mu_y*ky**4)*moments_i(ip,ij,ikx,iky,iz,updatetlevel) &
! Collision term
+ TColl_i(ip,ij,ikx,iky,iz)
--
GitLab