diff --git a/bc_50ns_T10Z6.csv b/bc_50ns_T10Z6.csv
new file mode 100644
index 0000000..98fb235
--- /dev/null
+++ b/bc_50ns_T10Z6.csv
@@ -0,0 +1,5 @@
+t (s),n (m^-3),u (m s^-1),T (eV),Z
+0.000000E-000,1.111111E+025,0.000000E+000,1.000000E+001,6.000000E+000
+5.000000E-008,1.111111E+025,0.000000E+000,1.000000E+001,6.000000E+000
+5.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,6.000000E+000
+2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,6.000000E+000
diff --git a/bc_50ns_T30Z11.csv b/bc_50ns_T30Z11.csv
new file mode 100644
index 0000000..856d5d6
--- /dev/null
+++ b/bc_50ns_T30Z11.csv
@@ -0,0 +1,5 @@
+t (s),n (m^-3),u (m s^-1),T (eV),Z
+0.000000E-000,1.111111E+025,0.000000E+000,3.000000E+001,1.100000E+001
+5.000000E-008,1.111111E+025,0.000000E+000,3.000000E+001,1.100000E+001
+5.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,1.100000E+001
+2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,1.100000E+001
diff --git a/bc_50ns_T60Z16.csv b/bc_50ns_T60Z16.csv
new file mode 100644
index 0000000..a3eff24
--- /dev/null
+++ b/bc_50ns_T60Z16.csv
@@ -0,0 +1,5 @@
+t (s),n (m^-3),u (m s^-1),T (eV),Z
+0.000000E-000,1.111111E+025,0.000000E+000,6.000000E+001,1.600000E+001
+5.000000E-008,1.111111E+025,0.000000E+000,6.000000E+001,1.600000E+001
+5.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,1.600000E+001
+2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,1.600000E+001
diff --git a/bc_50ns_T6Z4.csv b/bc_50ns_T6Z4.csv
new file mode 100644
index 0000000..4d623c7
--- /dev/null
+++ b/bc_50ns_T6Z4.csv
@@ -0,0 +1,5 @@
+t (s),n (m^-3),u (m s^-1),T (eV),Z
+0.000000E-000,1.111111E+025,0.000000E+000,6.000000E+000,4.000000E+000
+5.000000E-008,1.111111E+025,0.000000E+000,6.000000E+000,4.000000E+000
+5.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,4.000000E+000
+2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,4.000000E+000
diff --git a/bc_5ns.csv b/bc_5ns.csv
deleted file mode 100644
index 306b9b5..0000000
--- a/bc_5ns.csv
+++ /dev/null
@@ -1,5 +0,0 @@
-t (s),n (m^-3),u (m s^-1),T (eV),Z
-0.000000E-000,1.111111E+025,0.000000E+000,3.000000E+001,9.000000E+000
-5.000000E-008,1.111111E+025,0.000000E+000,3.000000E+001,9.000000E+000
-5.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,9.000000E+000
-2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,9.000000E+000
diff --git a/bc_fa_T30Z11.csv b/bc_fa_T30Z11.csv
new file mode 100644
index 0000000..2c0c946
--- /dev/null
+++ b/bc_fa_T30Z11.csv
@@ -0,0 +1,3 @@
+t (s),n (m^-3),u (m s^-1),T (eV),Z
+0.000000E-000,1.111111E+025,0.000000E+000,3.000000E+001,1.100000E+001
+1.000000E-006,1.111111E+025,0.000000E+000,3.000000E+001,1.100000E+001
diff --git a/bc_fullAblation.csv b/bc_fullAblation.csv
deleted file mode 100644
index 487aee3..0000000
--- a/bc_fullAblation.csv
+++ /dev/null
@@ -1,3 +0,0 @@
-t (s),n (m^-3),u (m s^-1),T (eV),Z
-0.000000E-000,1.111111E+025,0.000000E+000,3.000000E+001,9.000000E+000
-1.000000E-006,1.111111E+025,0.000000E+000,3.000000E+001,9.000000E+000
diff --git a/plasmaExpansion.f90 b/plasmaExpansion.f90
index fd24795..f715510 100644
--- a/plasmaExpansion.f90
+++ b/plasmaExpansion.f90
@@ -93,7 +93,7 @@ program plasmaExpansion
   phi_ref  = kb * Temp_ref / qe
 
   ! Set input parameters (remember these have to be in non-dimensional units)
-  c_s      = sqrt(9.0_dp * gamma_i * 1.0_dp)
+  c_s      = sqrt(11.0_dp * gamma_i * 1.0_dp)
   bc_file = 'bc.csv'
   call boundaryConditions%init(bc_file)