From 17021f819412f48432a9fbf8f88a59957c2f80d8 Mon Sep 17 00:00:00 2001
From: JGonzalez <jorge.gonzalez@upm.es>
Date: Wed, 20 Nov 2024 13:58:16 +0100
Subject: [PATCH] Adjust length and simulation time

---
 plasmaExpansion.f90 | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/plasmaExpansion.f90 b/plasmaExpansion.f90
index cd103d1..20b456c 100644
--- a/plasmaExpansion.f90
+++ b/plasmaExpansion.f90
@@ -99,7 +99,7 @@ program plasmaExpansion
 
   ! Set domain boundaries (non-dimensional units)
   r0 =  10.0e-6_dp / L_ref
-  rf =   1.0e-3_dp / L_ref
+  rf =   2.0e-3_dp / L_ref
   dr = 5.0e2_dp
   nr = nint((rf - r0) / dr) + 1
   dr =      (rf - r0) / float(nr-1)
@@ -109,7 +109,7 @@ program plasmaExpansion
   end do
 
   ! Set position to calculate cumulative sum of f (non-dimensional units)
-  rCum =  5.0e-4 / L_ref
+  rCum =  1.0e-3 / L_ref
 
   ! Index for cumulative sum
   rCum_index = minloc(abs(r - rCum), 1)
@@ -131,7 +131,7 @@ program plasmaExpansion
   end if
 
   t0 = 0.0_dp
-  tf = 5.0e-7_dp / t_ref
+  tf = 2.0e-7_dp / t_ref
   ! tf = 1.0e1_dp * (rf - r0) / c_s
   dt = 1.0e-2_dp*dr/c_s
   nt = nint((tf - t0) / dt)