From a2aae1bcd96a6d9d3f4cb83523b6889b1726d344 Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Fri, 17 Jun 2022 14:39:58 +0200
Subject: [PATCH] ppj init tests
---
src/inital.F90 | 24 ++++++++++++------------
1 file changed, 12 insertions(+), 12 deletions(-)
diff --git a/src/inital.F90 b/src/inital.F90
index 2bdcfe0a..c0261105 100644
--- a/src/inital.F90
+++ b/src/inital.F90
@@ -393,11 +393,11 @@ SUBROUTINE init_ppj
IMPLICIT NONE
REAL(dp) :: noise
- REAL(dp) :: kx, ky, sigma_z, amplitude, ky_shift, z
+ REAL(dp) :: kx, ky, sigma_z, amp, ky_shift, z
INTEGER, DIMENSION(12) :: iseedarr
sigma_z = pi/4.0
- amplitude = 0.1
+ amp = 1e4
!**** Broad noise initialization *******************************************
! Electrons
@@ -419,15 +419,15 @@ 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,:) = 0._dp!(kx_min/(ABS(kx)+kx_min))*(ky_min/(ABS(ky)+ky_min))
ELSE
- moments_e(ip,ij,iky,ikx,iz,:) = 0.5_dp*(kx_min/(ABS(kx)+kx_min))
+ moments_e(ip,ij,iky,ikx,iz,:) = 0.5_dp*amp*(kx_min/(ABS(kx)+kx_min))
ENDIF
ENDIF
- ! z-dep
+ CALL RANDOM_NUMBER(noise)
+ ! z-dep and noise
moments_e(ip,ij,iky,ikx,iz,:) = moments_e(ip,ij,iky,ikx,iz,:) * &
- ! (1 + exp(-(z/sigma_z)**2/2.0)*sqrt(2.0*sqrt(pi)/sigma_z))
- (Jacobian(iz,0)*iInt_Jacobian)**2
+ (Jacobian(iz,0)*iInt_Jacobian)**2 !+ (init_background + init_noiselvl*(noise-0.5_dp))
END DO
END DO
END DO
@@ -462,15 +462,15 @@ 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,:) = 0._dp!(kx_min/(ABS(kx)+kx_min))*(ky_min/(ABS(ky)+ky_min))
ELSE
- moments_i(ip,ij,iky,ikx,iz,:) = 0.5_dp*(kx_min/(ABS(kx)+kx_min))
+ moments_i(ip,ij,iky,ikx,iz,:) = 0.5_dp*amp*(kx_min/(ABS(kx)+kx_min))
ENDIF
ENDIF
- ! z-dep
+ CALL RANDOM_NUMBER(noise)
+ ! z-dep and noise
moments_i(ip,ij,iky,ikx,iz,:) = moments_i(ip,ij,iky,ikx,iz,:) * &
- ! (1 + exp(-(z/sigma_z)**2/2.0)*sqrt(2.0*sqrt(pi)/sigma_z))
- (Jacobian(iz,0)*iInt_Jacobian)**2
+ (Jacobian(iz,0)*iInt_Jacobian)**2 !+ (init_background + init_noiselvl*(noise-0.5_dp))
END DO
END DO
END DO
--
GitLab