Skip to content
Snippets Groups Projects
Commit 1ac9de43 authored by Antoine Cyril David Hoffmann's avatar Antoine Cyril David Hoffmann
Browse files

closure and kernel model flags

parent b82803a0
No related branches found
No related tags found
No related merge requests found
...@@ -6,7 +6,8 @@ MODULE model ...@@ -6,7 +6,8 @@ MODULE model
PRIVATE PRIVATE
INTEGER, PUBLIC, PROTECTED :: CO = 0 ! Collision Operator INTEGER, PUBLIC, PROTECTED :: CO = 0 ! Collision Operator
LOGICAL, PUBLIC, PROTECTED :: DK = 0 ! Drift kinetic 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 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 ! Hyperdiffusivity coefficient (for num. stability)
REAL(dp), PUBLIC, PROTECTED :: nu = 1._dp ! Collision frequency REAL(dp), PUBLIC, PROTECTED :: nu = 1._dp ! Collision frequency
...@@ -44,7 +45,7 @@ CONTAINS ...@@ -44,7 +45,7 @@ CONTAINS
USE prec_const USE prec_const
IMPLICIT NONE IMPLICIT NONE
NAMELIST /MODEL_PAR/ CO, DK, NON_LIN, mu, nu, tau_e, tau_i, sigma_e, sigma_i, & NAMELIST /MODEL_PAR/ CO, CLOS, KERN, NON_LIN, mu, nu, tau_e, tau_i, sigma_e, sigma_i, &
q_e, q_i, eta_n, eta_T, eta_B, lambdaD q_e, q_i, eta_n, eta_T, eta_B, lambdaD
READ(lu_in,model_par) READ(lu_in,model_par)
...@@ -88,6 +89,8 @@ CONTAINS ...@@ -88,6 +89,8 @@ CONTAINS
INTEGER, INTENT(in) :: fidres INTEGER, INTENT(in) :: fidres
CHARACTER(len=256), INTENT(in) :: str CHARACTER(len=256), INTENT(in) :: str
CALL attach(fidres, TRIM(str), "CO", CO) CALL attach(fidres, TRIM(str), "CO", CO)
CALL attach(fidres, TRIM(str), "CLOS", CLOS)
CALL attach(fidres, TRIM(str), "KERN", KERN)
CALL attach(fidres, TRIM(str), "NON_LIN", NON_LIN) CALL attach(fidres, TRIM(str), "NON_LIN", NON_LIN)
CALL attach(fidres, TRIM(str), "nu", nu) CALL attach(fidres, TRIM(str), "nu", nu)
CALL attach(fidres, TRIM(str), "tau_e", tau_e) CALL attach(fidres, TRIM(str), "tau_e", tau_e)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment