From a58f14c587be4a93c54897128fe156c9e3205c75 Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Wed, 4 Aug 2021 14:01:54 +0200
Subject: [PATCH] Now working for 1D linear analysis
---
src/geometry_mod.F90 | 10 +++++++++-
1 file changed, 9 insertions(+), 1 deletion(-)
diff --git a/src/geometry_mod.F90 b/src/geometry_mod.F90
index 02a9edb5..4703bf1c 100644
--- a/src/geometry_mod.F90
+++ b/src/geometry_mod.F90
@@ -63,7 +63,15 @@ contains
Ckxky(ikx,iky,iz) = -SIN(z)*kx -(COS(z)+shear*z*SIN(z))*ky
ENDDO
ENDDO
-
+ IF (Nky .EQ. 1) THEN ! linear 1D run we switch kx and ky for parallel opt
+ DO iky = ikys, ikye
+ ky = kyarray(iky)
+ DO ikx= ikxs, ikxe
+ kx = kxarray(ikx)
+ Ckxky(ikx,iky,iz) = -SIN(z)*ky -(COS(z)+shear*z*SIN(z))*kx
+ ENDDO
+ ENDDO
+ ENDIF
ENDDO zloop
END SUBROUTINE eval_s_alpha_geometry
--
GitLab