From 43fea64986fc4f65ee61baab57f060df9ceedc9e Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Fri, 16 Feb 2024 16:08:12 +0100
Subject: [PATCH] species-wise hyperdiff
---
src/species_mod.F90 | 8 +++++++-
1 file changed, 7 insertions(+), 1 deletion(-)
diff --git a/src/species_mod.F90 b/src/species_mod.F90
index 2d86e06..43c8557 100644
--- a/src/species_mod.F90
+++ b/src/species_mod.F90
@@ -10,6 +10,7 @@ MODULE species
REAL(xp) :: q_ ! Charge
REAL(xp) :: k_N_ ! density drive (L_ref/L_Ni)
REAL(xp) :: k_T_ ! temperature drive (L_ref/L_Ti)
+ REAL(xp) :: mu_ ! species-wise hyperdiffusion (not tested)
!! Arrays to store all species features
CHARACTER(len=32),&
ALLOCATABLE, DIMENSION(:), PUBLIC, PROTECTED :: name ! name of the species
@@ -19,6 +20,7 @@ MODULE species
REAL(xp), ALLOCATABLE, DIMENSION(:), PUBLIC, PROTECTED :: k_N ! density drive (L_ref/L_Ni)
REAL(xp), ALLOCATABLE, DIMENSION(:), PUBLIC, PROTECTED :: k_T ! temperature drive (L_ref/L_Ti)
REAL(xp), ALLOCATABLE, DIMENSION(:,:),PUBLIC, PROTECTED :: nu_ab ! Collision frequency tensor
+ REAL(xp), ALLOCATABLE, DIMENSION(:), PUBLIC, PROTECTED :: mu ! Hyperdiffusion
!! Auxiliary variables to store precomputation
REAL(xp), ALLOCATABLE, DIMENSION(:),PUBLIC, PROTECTED :: tau_q ! factor of the magnetic moment coupling
REAL(xp), ALLOCATABLE, DIMENSION(:),PUBLIC, PROTECTED :: q_tau ! charge/temp ratio
@@ -43,7 +45,7 @@ CONTAINS
INTEGER :: ia,ib
! expected namelist in the input file
NAMELIST /SPECIES/ &
- name_, tau_, sigma_, q_, k_N_, k_T_
+ name_, tau_, sigma_, q_, k_N_, k_T_, mu_
! allocate the arrays of species parameters
CALL species_allocate
! loop over the species namelists in the input file
@@ -56,6 +58,7 @@ CONTAINS
q_ = 1._xp
k_N_ = 2.22_xp
k_T_ = 6.96_xp
+ mu_ = 0._xp
! read input
READ(lu_in,species)
! place values found in the arrays
@@ -66,6 +69,7 @@ CONTAINS
k_N(ia) = k_N_
k_T(ia) = k_T_
tau_q(ia) = tau_/q_
+ mu(ia) = mu_
! precompute factors
q_tau(ia) = q_/tau_
sqrtTau_q(ia) = sqrt(tau_)/q_
@@ -130,6 +134,7 @@ CONTAINS
CALL attach(fid, TRIM(str), "q", q(ia))
CALL attach(fid, TRIM(str), "k_N", k_N(ia))
CALL attach(fid, TRIM(str), "k_T", k_T(ia))
+ CALL attach(fid, TRIM(str), "mu", mu(ia))
ENDDO
END SUBROUTINE species_outputinputs
@@ -144,6 +149,7 @@ CONTAINS
ALLOCATE( q(Na))
ALLOCATE( k_N(Na))
ALLOCATE( k_T(Na))
+ ALLOCATE( mu(Na))
ALLOCATE( tau_q(Na))
ALLOCATE( q_tau(Na))
ALLOCATE( sqrtTau_q(Na))
--
GitLab