From 6796aee960a28d65759ab352485fa95a0f7677b8 Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Wed, 21 Sep 2022 11:48:14 +0200
Subject: [PATCH] flag for EM effects and separation of species background
gradients
---
src/model_mod.F90 | 24 +++++++++++++++---------
1 file changed, 15 insertions(+), 9 deletions(-)
diff --git a/src/model_mod.F90 b/src/model_mod.F90
index a7209fbe..75a27687 100644
--- a/src/model_mod.F90
+++ b/src/model_mod.F90
@@ -24,16 +24,17 @@ MODULE model
REAL(dp), PUBLIC, PROTECTED :: sigma_i = 1._dp !
REAL(dp), PUBLIC, PROTECTED :: q_e = -1._dp ! Charge
REAL(dp), PUBLIC, PROTECTED :: q_i = 1._dp !
- REAL(dp), PUBLIC, PROTECTED :: k_N = 1._dp ! Ion density drive
- REAL(dp), PUBLIC, PROTECTED :: eta_N = 1._dp ! electron-ion density drive ratio (k_Ne/k_Ni)
- REAL(dp), PUBLIC, PROTECTED :: k_T = 0._dp ! Temperature drive
- REAL(dp), PUBLIC, PROTECTED :: eta_T = 0._dp ! electron-ion temperature drive ratio (k_Te/k_Ti)
+ REAL(dp), PUBLIC, PROTECTED :: k_Ni = 7._dp ! Ion density drive
+ REAL(dp), PUBLIC, PROTECTED :: k_Ne = 7._dp ! Ele ''
+ REAL(dp), PUBLIC, PROTECTED :: k_Ti = 2._dp ! Ion temperature drive
+ REAL(dp), PUBLIC, PROTECTED :: k_Te = 2._dp ! Ele ''
REAL(dp), PUBLIC, PROTECTED :: K_E = 0._dp ! Backg. electric field drive
REAL(dp), PUBLIC, PROTECTED :: GradB = 1._dp ! Magnetic gradient
REAL(dp), PUBLIC, PROTECTED :: CurvB = 1._dp ! Magnetic curvature
REAL(dp), PUBLIC, PROTECTED :: lambdaD = 1._dp ! Debye length
REAL(dp), PUBLIC, PROTECTED :: beta = 0._dp ! electron plasma Beta (8piNT_e/B0^2)
+ LOGICAL, PUBLIC, PROTECTED :: EM = .false. ! Electromagnetic effects flag
REAL(dp), PUBLIC, PROTECTED :: taue_qe ! factor of the magnetic moment coupling
REAL(dp), PUBLIC, PROTECTED :: taui_qi !
REAL(dp), PUBLIC, PROTECTED :: qi_taui !
@@ -64,7 +65,7 @@ CONTAINS
NAMELIST /MODEL_PAR/ CLOS, NL_CLOS, KERN, LINEARITY, KIN_E, &
mu_x, mu_y, N_HD, mu_z, mu_p, mu_j, nu,&
tau_e, tau_i, sigma_e, sigma_i, q_e, q_i,&
- k_N, eta_N, k_T, eta_T, K_E, GradB, CurvB, lambdaD, beta
+ k_Ne, k_Ni, k_Te, k_Ti, GradB, CurvB, lambdaD, beta
READ(lu_in,model_par)
@@ -73,6 +74,11 @@ CONTAINS
beta = 0._dp
ENDIF
+ IF(beta .GT. 0) THEN
+ EM = .TRUE.
+ IF(my_id.EQ.0) print*, 'Electromagnetic effects are included'
+ ENDIF
+
taue_qe = tau_e/q_e ! factor of the magnetic moment coupling
taui_qi = tau_i/q_i ! factor of the magnetic moment coupling
qe_taue = q_e/tau_e
@@ -134,10 +140,10 @@ CONTAINS
CALL attach(fidres, TRIM(str), "sigma_i", sigma_i)
CALL attach(fidres, TRIM(str), "q_e", q_e)
CALL attach(fidres, TRIM(str), "q_i", q_i)
- CALL attach(fidres, TRIM(str), "k_N", k_N)
- CALL attach(fidres, TRIM(str), "eta_N", eta_N)
- CALL attach(fidres, TRIM(str), "k_T", k_T)
- CALL attach(fidres, TRIM(str), "eta_T", eta_T)
+ CALL attach(fidres, TRIM(str), "k_Ne", k_Ne)
+ CALL attach(fidres, TRIM(str), "k_Ni", k_Ni)
+ CALL attach(fidres, TRIM(str), "k_Te", k_Te)
+ CALL attach(fidres, TRIM(str), "k_Ti", k_Ti)
CALL attach(fidres, TRIM(str), "K_E", K_E)
CALL attach(fidres, TRIM(str), "lambdaD", lambdaD)
CALL attach(fidres, TRIM(str), "beta", beta)
--
GitLab