From d543ce4316d37b45a6c8c269b3d6f6eff44dce32 Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Wed, 21 Sep 2022 11:36:12 +0200
Subject: [PATCH] init ppj is now deltak scaled
---
src/inital.F90 | 10 +++++-----
1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/src/inital.F90 b/src/inital.F90
index d1958e72..dbdc5007 100644
--- a/src/inital.F90
+++ b/src/inital.F90
@@ -221,7 +221,7 @@ SUBROUTINE init_gyrodens
END DO
IF ( contains_ky0 ) THEN
DO ikx=2,Nkx/2 !symmetry at ky = 0 for all z
- moments_e(ip,ij,iky_0,ikx,:,:) = moments_e( ip,ij,iky_0,Nkx+2-ikx,:, :)
+ moments_e(ip,ij,iky_0,ikx,:,:) = moments_e(ip,ij,iky_0,Nkx+2-ikx,:,:)
END DO
ENDIF
END DO
@@ -422,9 +422,9 @@ SUBROUTINE init_ppj
ENDIF
ELSE
IF(ky .GT. 0) THEN
- moments_e(ip,ij,iky,ikx,iz,:) = (kx_min/(ABS(kx)+kx_min))*(ky_min/(ABS(ky)+ky_min))
+ moments_e(ip,ij,iky,ikx,iz,:) = (deltakx/(ABS(kx)+deltakx))*(ky_min/(ABS(ky)+ky_min))
ELSE
- moments_e(ip,ij,iky,ikx,iz,:) = 0.5_dp*amp*(kx_min/(ABS(kx)+kx_min))
+ moments_e(ip,ij,iky,ikx,iz,:) = 0.5_dp*amp*(deltakx/(ABS(kx)+deltakx))
ENDIF
ENDIF
CALL RANDOM_NUMBER(noise)
@@ -465,9 +465,9 @@ SUBROUTINE init_ppj
ENDIF
ELSE
IF(ky .GT. 0) THEN
- moments_i(ip,ij,iky,ikx,iz,:) = (kx_min/(ABS(kx)+kx_min))*(ky_min/(ABS(ky)+ky_min))
+ moments_i(ip,ij,iky,ikx,iz,:) = (deltakx/(ABS(kx)+deltakx))*(ky_min/(ABS(ky)+ky_min))
ELSE
- moments_i(ip,ij,iky,ikx,iz,:) = 0.5_dp*amp*(kx_min/(ABS(kx)+kx_min))
+ moments_i(ip,ij,iky,ikx,iz,:) = 0.5_dp*amp*(deltakx/(ABS(kx)+deltakx))
ENDIF
ENDIF
CALL RANDOM_NUMBER(noise)
--
GitLab