Restore definition of ne for fluxlim case

Fixes FusionDLS#27
PlasmaFAIR · Jul 11, 2024 · 4ad12e1 · 4ad12e1
1 parent 2729397
commit 4ad12e1
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/fusiondls/Iterate.py b/src/fusiondls/Iterate.py
@@ -73,6 +73,8 @@ def LengFunc(s, y, si, st):
     # Flux limiter
     dtds = 0
     if radios["fluxlim"]:
+        # set density using constant pressure assumption (missing factor of 2 at target due to lack of Bohm condition)
+        ne = nu * Tu / T
         dtds = (
             qoverB
             * fieldValue