From 2b6ea0ff8494693d9f5f7e5bfb3e5c3ec80472fa Mon Sep 17 00:00:00 2001
From: JGonzalez <jorge.gonzalez@upm.es>
Date: Tue, 3 Feb 2026 11:48:43 +0100
Subject: [PATCH] First attempt at working with a distribution function with u

---
 src/modules/mesh/moduleMeshBoundary.f90 | 14 ++++++++++----
 1 file changed, 10 insertions(+), 4 deletions(-)

diff --git a/src/modules/mesh/moduleMeshBoundary.f90 b/src/modules/mesh/moduleMeshBoundary.f90
index deb4459..eaa2d52 100644
--- a/src/modules/mesh/moduleMeshBoundary.f90
+++ b/src/modules/mesh/moduleMeshBoundary.f90
@@ -220,6 +220,7 @@ MODULE moduleMeshBoundary
 
     subroutine outflowAdaptive(edge, part)
       use moduleRandom
+      use moduleRefParam, only: v_ref
       implicit none
 
       class(meshEdge), intent(inout):: edge
@@ -228,14 +229,19 @@ MODULE moduleMeshBoundary
       select type(bound => edge%boundary%bTypes(part%species%n)%obj)
       type is(boundaryOutflowAdaptive)
 
-        if (random() < 0.844d0) then
-          call reflection(edge, part)
-
-        else
+        ! if (random() < 0.844d0) then
+        call reflection(edge, part)
+        part%v = part%v + 40e3/v_ref
+        if (dot_product(part%v, edge%normal) <= 0.d0) then
           call transparent(edge, part)
 
         end if
 
+        ! else
+          ! call transparent(edge, part)
+
+        ! end if
+
       end select
 
     end subroutine outflowAdaptive