From ee6a8760a3efdd8806fba9ac8bd2d4573cef62be Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Fri, 19 Feb 2021 11:24:12 +0100
Subject: [PATCH] added parameters for kinetic hyperdiffusion
---
src/model_mod.F90 | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/src/model_mod.F90 b/src/model_mod.F90
index 264a640e..00fb21cb 100644
--- a/src/model_mod.F90
+++ b/src/model_mod.F90
@@ -9,7 +9,9 @@ MODULE model
INTEGER, PUBLIC, PROTECTED :: CLOS = 0 ! Truncation method
INTEGER, PUBLIC, PROTECTED :: KERN = 0 ! Kernel model
LOGICAL, PUBLIC, PROTECTED :: NON_LIN = .true. ! To turn on non linear bracket term
- REAL(dp), PUBLIC, PROTECTED :: mu = 0._dp ! Hyperdiffusivity coefficient (for num. stability)
+ REAL(dp), PUBLIC, PROTECTED :: mu = 0._dp ! spatial Hyperdiffusivity coefficient (for num. stability)
+ REAL(dp), PUBLIC, PROTECTED :: mu_p = 0._dp ! kinetic para hyperdiffusivity coefficient (for num. stability)
+ REAL(dp), PUBLIC, PROTECTED :: mu_j = 0._dp ! kinetic perp hyperdiffusivity coefficient (for num. stability)
REAL(dp), PUBLIC, PROTECTED :: nu = 1._dp ! Collision frequency
REAL(dp), PUBLIC, PROTECTED :: tau_e = 1._dp ! Temperature
REAL(dp), PUBLIC, PROTECTED :: tau_i = 1._dp !
@@ -45,7 +47,7 @@ CONTAINS
USE prec_const
IMPLICIT NONE
- NAMELIST /MODEL_PAR/ CO, CLOS, KERN, NON_LIN, mu, nu, tau_e, tau_i, sigma_e, sigma_i, &
+ NAMELIST /MODEL_PAR/ CO, CLOS, KERN, NON_LIN, mu, mu_p, mu_j, nu, tau_e, tau_i, sigma_e, sigma_i, &
q_e, q_i, eta_n, eta_T, eta_B, lambdaD
READ(lu_in,model_par)
--
GitLab