From 4ad67ff778fab3cb60ce46204aa4974dcbdd61f3 Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Wed, 26 May 2021 13:54:59 +0200
Subject: [PATCH] Correction of the e-e collision frequency in dougherty
---
src/collision_mod.F90 | 4 ++--
src/model_mod.F90 | 8 ++++++++
2 files changed, 10 insertions(+), 2 deletions(-)
diff --git a/src/collision_mod.F90 b/src/collision_mod.F90
index 2bdeb5ce..85db678e 100644
--- a/src/collision_mod.F90
+++ b/src/collision_mod.F90
@@ -111,8 +111,8 @@ CONTAINS
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ENDIF
- ! Multiply by electron-ion collision coefficient
- TColl_ = nu_e * TColl_
+ ! Multiply by electron-electron collision coefficient
+ TColl_ = nu_ee * TColl_
END SUBROUTINE DoughertyGK_e
diff --git a/src/model_mod.F90 b/src/model_mod.F90
index 3b77badc..4c9e298c 100644
--- a/src/model_mod.F90
+++ b/src/model_mod.F90
@@ -76,6 +76,14 @@ CONTAINS
nu_i = nu ! ion-ion collision frequ.
nu_ee = nu_e ! e-e coll. frequ.
nu_ie = nu_i ! i-e coll. frequ.
+
+
+ ! Old normalization (MOLI Jorge/Frei)
+ ! nu_e = 0.532_dp*nu ! electron-ion collision frequency (where already multiplied by 0.532)
+ ! nu_i = 0.532_dp*nu*sigma_e*tau_e**(-3._dp/2._dp)/SQRT2 ! ion-ion collision frequ.
+ ! nu_ee = nu_e/SQRT2 ! e-e coll. frequ.
+ ! nu_ie = 0.532_dp*nu*sigma_e**2 ! i-e coll. frequ.
+
END SUBROUTINE model_readinputs
--
GitLab