From 2931eaff00c741c4e39e8a60e2d683375256e879 Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Fri, 15 Dec 2023 14:38:45 +0100
Subject: [PATCH] ExB shear flow is back
---
src/ExB_shear_flow_mod.F90 | 10 +++++-----
1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/src/ExB_shear_flow_mod.F90 b/src/ExB_shear_flow_mod.F90
index dadf1aa..8d76ef3 100644
--- a/src/ExB_shear_flow_mod.F90
+++ b/src/ExB_shear_flow_mod.F90
@@ -26,7 +26,7 @@ CONTAINS
USE grid, ONLY: Nx, local_nky, total_nky, local_nx, Ny, deltakx, deltaky,&
kx_max, kx_min !kyarray, kyarray_full
USE geometry, ONLY: Cyq0_x0, C_y
- ! USE basic, ONLY: dt
+ USE basic, ONLY: speak
USE model, ONLY: LINEARITY
IMPLICIT NONE
INTEGER :: iky
@@ -36,6 +36,7 @@ CONTAINS
! In GENE, there is a minus sign here...
gamma_E = ExBrate*C_y*abs(Cyq0_x0/C_y)
IF(abs(gamma_E) .GT. EPSILON(gamma_E)) THEN
+ CALL speak('-ExB background flow detected-')
ExB = .TRUE.
t0 = deltakx/deltaky/gamma_E
inv_t0 = 1._xp/t0
@@ -70,7 +71,7 @@ CONTAINS
ALLOCATE(shiftnow_ExB(local_nky))
shiftnow_ExB = .FALSE.
- ! Setup nonlinear factor
+ ! Setup nonlinear factor (McMillan 2019)
ALLOCATE( ExB_NL_factor(Nx,local_nky))
ALLOCATE(inv_ExB_NL_factor(Ny/2+1,local_nx))
ExB_NL_factor = 1._xp
@@ -128,7 +129,6 @@ CONTAINS
sky_ExB_full(iky) = sky_ExB_full(iky) - REAL(iky-1,xp)*deltaky*gamma_E*dt_sub
ENDDO
! Shift the arrays if the shear value sky is too high
- IF(LINEARITY .EQ. 'nonlinear') &
CALL Array_shift_ExB_shear_flow
! We update the operators and grids
@@ -139,8 +139,8 @@ CONTAINS
CALL evaluate_EM_op
CALL evaluate_magn_curv
! update the ExB nonlinear factor...
- IF(LINEARITY .EQ. 'nonlinear') &
- CALL update_nonlinear_ExB_factors(dt_sub)
+ ! IF(LINEARITY .EQ. 'nonlinear') &
+ ! CALL Update_nonlinear_ExB_factors(dt_sub)
ENDIF
ENDIF
END SUBROUTINE Update_ExB_shear_flow
--
GitLab