From 313f33a2e4a1cc49f500658b74890a63651e4587 Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Mon, 13 Feb 2023 13:45:25 +0100
Subject: [PATCH] put a new BC for sheared and renamed the former periodic as
cyclic
---
src/geometry_mod.F90 | 10 +++++++---
1 file changed, 7 insertions(+), 3 deletions(-)
diff --git a/src/geometry_mod.F90 b/src/geometry_mod.F90
index 7105088d..41c42936 100644
--- a/src/geometry_mod.F90
+++ b/src/geometry_mod.F90
@@ -369,7 +369,9 @@ CONTAINS
SELECT CASE(parallel_bc)
CASE ('dirichlet')! connected to 0
ikx_zBC_L(iky,ikx) = -99
- CASE ('periodic') !reroute it by cycling through modes
+ CASE ('periodic')
+ ikx_zBC_L(iky,ikx) = ikx
+ CASE ('cyclic')! reroute it by cycling through modes
ikx_zBC_L(iky,ikx) = MODULO(ikx_zBC_L(iky,ikx)-1,Nkx)+1
END SELECT
ENDIF
@@ -395,8 +397,10 @@ CONTAINS
SELECT CASE(parallel_bc)
CASE ('dirichlet') ! connected to 0
ikx_zBC_R(iky,ikx) = -99
- CASE ('periodic') !reroute it by cycling through modes
- write(*,*) 'check',ikx,iky, kxarray(ikx) + shift, '>', kx_max
+ CASE ('periodic') ! connected to itself as for shearless
+ ikx_zBC_R(iky,ikx) = ikx
+ CASE ('cyclic')
+ ! write(*,*) 'check',ikx,iky, kxarray(ikx) + shift, '>', kx_max
ikx_zBC_R(iky,ikx) = MODULO(ikx_zBC_R(iky,ikx)-1,Nkx)+1
END SELECT
ENDIF
--
GitLab