From e382c96233f1d5cf3bfb4548ffdc46dc6ad47131 Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Fri, 8 Dec 2023 15:59:00 +0100
Subject: [PATCH] add an ion charge for adiabatic ion runs
---
src/model_mod.F90 | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/model_mod.F90 b/src/model_mod.F90
index 1eff2bcf..7f23ce85 100644
--- a/src/model_mod.F90
+++ b/src/model_mod.F90
@@ -20,6 +20,7 @@ MODULE model
LOGICAL, PUBLIC :: ADIAB_E = .false. ! adiabatic electron model
LOGICAL, PUBLIC :: ADIAB_I = .false. ! adiabatic ion model
REAL(xp), PUBLIC, PROTECTED :: tau_i = 1.0 ! electron-ion temperature ratio for ion adiabatic model
+ REAL(xp), PUBLIC, PROTECTED :: q_i = 1.0 ! ion charge for ion adiabatic model
REAL(xp), PUBLIC, PROTECTED :: nu = 0._xp ! collision frequency parameter
REAL(xp), PUBLIC, PROTECTED :: k_gB = 1._xp ! Magnetic gradient strength (L_ref/L_gB)
REAL(xp), PUBLIC, PROTECTED :: k_cB = 1._xp ! Magnetic curvature strength (L_ref/L_cB)
@@ -51,7 +52,7 @@ CONTAINS
IMPLICIT NONE
NAMELIST /MODEL/ KERN, LINEARITY, RM_LD_T_EQ, FORCE_SYMMETRY, MHD_PD, &
- Na, ADIAB_E, ADIAB_I, tau_i, &
+ Na, ADIAB_E, ADIAB_I, q_i, tau_i, &
mu_x, mu_y, N_HD, HDz_h, mu_z, mu_p, mu_j, HYP_V, &
nu, k_gB, k_cB, lambdaD, beta, ExBrate, ExB_NL_CORRECTION,&
ikxZF, ZFamp
--
GitLab