From 69215ef66d52557f84232c8ffedd906d0cb2bc90 Mon Sep 17 00:00:00 2001
From: Jorge Gonzalez - Galactica <jorge.gonzalez@upm.es>
Date: Tue, 22 Jul 2025 19:52:39 +0200
Subject: [PATCH] Change in calculating ionization

I don't know why I normalizing density n_0 by Vol_ref and not n_ref
---
 src/modules/moduleBoundary.f90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/modules/moduleBoundary.f90 b/src/modules/moduleBoundary.f90
index 0b76105..b5bda2f 100644
--- a/src/modules/moduleBoundary.f90
+++ b/src/modules/moduleBoundary.f90
@@ -126,7 +126,7 @@ MODULE moduleBoundary
       SELECT TYPE(boundary)
       TYPE IS(boundaryIonization)
         boundary%m0 = m0 / m_ref
-        boundary%n0 = n0 * Vol_ref
+        boundary%n0 = n0 / n_ref
         boundary%v0 = v0 / v_ref
         boundary%vTh = DSQRT(kb*T0/m0)/v_ref
         boundary%species => species(ion)%obj