From b4ff41187958413c55c2e7f596191d06c1f26aef Mon Sep 17 00:00:00 2001
From: Antoine <antoine.hoffmann@epfl.ch>
Date: Wed, 20 Sep 2023 09:16:57 +0200
Subject: [PATCH] Add params to initialize single mode
---
src/initial_mod.F90 | 14 +++++++++-----
1 file changed, 9 insertions(+), 5 deletions(-)
diff --git a/src/initial_mod.F90 b/src/initial_mod.F90
index c4ab7b9a..befe8f65 100644
--- a/src/initial_mod.F90
+++ b/src/initial_mod.F90
@@ -18,6 +18,8 @@ MODULE initial
REAL(xp), PUBLIC, PROTECTED :: init_noiselvl = 1E-6_xp
! Initialization for random number generator
INTEGER, PUBLIC, PROTECTED :: iseed=42
+ ! Single mode initialization
+ INTEGER, PUBLIC, PROTECTED :: ikx_init,iky_init
PUBLIC :: initial_outputinputs, initial_readinputs, initialize
@@ -31,7 +33,8 @@ CONTAINS
IMPLICIT NONE
NAMELIST /INITIAL/ INIT_OPT,ACT_ON_MODES,&
- init_amp,init_background,init_noiselvl,iseed
+ init_amp,init_background,init_noiselvl,iseed,&
+ ikx_init,iky_init
READ(lu_in,initial)
END SUBROUTINE initial_readinputs
@@ -271,10 +274,11 @@ CONTAINS
IMPLICIT NONE
moments = 0._xp
IF (my_id .EQ. 0) THEN
- moments(:,:,:,2,2,:,:) = init_amp
- moments(:,:,:,3,3,:,:) = init_amp
- moments(:,:,:,4,4,:,:) = init_amp
- moments(:,:,:,5,5,:,:) = init_amp
+ moments(:,:,:,iky_init,ikx_init,:,:) = init_amp
+ ! moments(:,:,:,2,2,:,:) = init_amp
+ ! moments(:,:,:,4,3,:,:) = init_amp
+ ! moments(:,:,:,4,4,:,:) = init_amp
+ ! moments(:,:,:,5,5,:,:) = init_amp
ENDIF
END SUBROUTINE init_single_mode
!******************************************************************************!
--
GitLab