From 945aa653658fbe577d50c9b8058079400780ffbd Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Thu, 29 Sep 2022 17:59:08 +0200
Subject: [PATCH] added a collision cutoff frequency
---
src/collision_mod.F90 | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/src/collision_mod.F90 b/src/collision_mod.F90
index a9d496d7..b6291344 100644
--- a/src/collision_mod.F90
+++ b/src/collision_mod.F90
@@ -17,7 +17,7 @@ CHARACTER(len=2),PUBLIC,PROTECTED :: collision_model
LOGICAL, PUBLIC, PROTECTED :: gyrokin_CO =.true. ! activates GK effects on CO
LOGICAL, PUBLIC, PROTECTED :: interspecies =.true. ! activates interpecies collision
CHARACTER(len=128), PUBLIC :: mat_file ! COSOlver matrix file names
-
+REAL(dp), PUBLIC, PROTECTED :: collision_kcut = 100.0
LOGICAL, PUBLIC, PROTECTED :: cosolver_coll ! check if cosolver matrices are used
PUBLIC :: collision_readinputs, coll_outputinputs
@@ -33,7 +33,7 @@ CONTAINS
SUBROUTINE collision_readinputs
! Read the input parameters
IMPLICIT NONE
- NAMELIST /COLLISION_PAR/ collision_model, gyrokin_CO, interspecies, mat_file
+ NAMELIST /COLLISION_PAR/ collision_model, gyrokin_CO, interspecies, mat_file, collision_kcut
READ(lu_in,collision_par)
SELECT CASE(collision_model)
CASE ('LB') ! Lenhard bernstein
@@ -55,6 +55,7 @@ CONTAINS
CASE DEFAULT
ERROR STOP 'Error stop: collision model not recognized!!'
END SELECT
+ print*, collision_kcut
END SUBROUTINE collision_readinputs
@@ -840,7 +841,7 @@ CONTAINS
DO iky = ikys,ikye
DO iz = izs,ize
! Check for nonlinear case if we are in the anti aliased domain or the filtered one
- kperp_sim = kparray(iky,ikx,iz,0) ! current simulation kperp
+ kperp_sim = MIN(kparray(iky,ikx,iz,0),collision_kcut) ! current simulation kperp
! Find the interval in kp grid mat where kperp_sim is contained
! Loop over the whole kp mat grid to find the smallest kperp that is
! larger than the current kperp_sim (brute force...)
--
GitLab