From a51d92a0d1c3830e70b7dc8a6d056270ca6492ee Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Fri, 24 Jun 2022 09:40:29 +0200
Subject: [PATCH] shear setup
---
src/auxval.F90 | 2 +-
src/grid_mod.F90 | 8 ++++++--
2 files changed, 7 insertions(+), 3 deletions(-)
diff --git a/src/auxval.F90 b/src/auxval.F90
index 8ade4159..5db2df05 100644
--- a/src/auxval.F90
+++ b/src/auxval.F90
@@ -27,7 +27,7 @@ subroutine auxval
CALL set_jgrid ! perp kin
- CALL set_kxgrid ! radial modes (MPI distributed by FFTW)
+ CALL set_kxgrid(shear) ! radial modes (MPI distributed by FFTW)
CALL set_kygrid ! azymuthal modes
diff --git a/src/grid_mod.F90 b/src/grid_mod.F90
index f58cbb86..2645c259 100644
--- a/src/grid_mod.F90
+++ b/src/grid_mod.F90
@@ -341,12 +341,16 @@ CONTAINS
END DO
END SUBROUTINE set_kygrid
- SUBROUTINE set_kxgrid
+ SUBROUTINE set_kxgrid(shear)
USE prec_const
USE model, ONLY: LINEARITY
IMPLICIT NONE
+ REAL(dp), INTENT(IN) :: shear
INTEGER :: i_, counter
-
+ IF(shear .GT. 0._dp) THEN
+ IF(my_id.EQ.0) write(*,*) 'Magnetic shear detected: set up sheared kx grid..'
+ Lx = Ly/(2._dp*pi*shear)
+ ENDIF
Nkx = Nx;
! Local data
! Start and END indices of grid
--
GitLab