diff --git a/plasmaExpansion.f90 b/plasmaExpansion.f90
index 2a4a331..f8f9120 100644
--- a/plasmaExpansion.f90
+++ b/plasmaExpansion.f90
@@ -353,15 +353,15 @@ program plasmaExpansion
 
   ! Set input parameters (remember these have to be in non-dimensional units)
   Temp0    = 60.0_dp * eV_to_K / Temp_ref
-  TempF    = 10.0_dp * eV_to_K / Temp_ref
+  TempF    =  5.0_dp * eV_to_K / Temp_ref
   Zave0    = 12.0_dp
   ZaveF    =  3.0_dp
   n_ecr    = 1.0e19_dp * cm3_to_m3 / n_ref
   c_s      = sqrt(Zave0 * gam * Temp0)
-  u_bc0    = c_s!sqrt(Temp0)
+  u_bc0    = sqrt(Temp0)
   u_bcF    = sqrt(TempF)
   n_bc0    = n_ecr / Zave0
-  n_bcF    = 1.0e-1*n_bc0!n_ecr*1.0e-1 / T_to_Z(Temp0)
+  n_bcF    = 1.0e-1_dp*n_bc0!n_ecr*1.0e-1 / T_to_Z(Temp0)
 
   ! Set domain boundaries (non-dimensional units)
   r0 = 200.0e-6_dp / L_ref
@@ -404,7 +404,7 @@ program plasmaExpansion
   dt =      (tf - t0) / float(nt)
 
   t_bc0 = nint(100.0e-9_dp / t_ref / dt)
-  t_bcF = nint(105.0e-9_dp / t_ref / dt)
+  t_bcF = nint(200.0e-9_dp / t_ref / dt)
 
   everyOutput = nint(1.0e-9_dp/t_ref/dt)
   if (everyOutput  == 0) then