diff --git a/data/collisions/IO_e-Xe.dat b/data/collisions/IO_e-Xe.dat
new file mode 100644
index 0000000..3067578
--- /dev/null
+++ b/data/collisions/IO_e-Xe.dat
@@ -0,0 +1,52 @@
+# EL cross sections extracted from PROGRAM MAGBOLTZ, VERSION 7.1 JUNE 2004 www.lxcat.net/Biagi-v7.1
+# Relative energy (eV) cross section (m^2)
+1.21E+01 0
+1.41E+01 3.923E-21
+1.64E+01 1.194E-20
+1.91E+01 2.1E-20
+2.22E+01 2.946E-20
+2.58E+01 3.65E-20
+3.00E+01 4.185E-20
+3.49E+01 4.552E-20
+4.06E+01 4.766E-20
+4.72E+01 4.85E-20
+5.49E+01 4.828E-20
+6.39E+01 5.031E-20
+7.43E+01 5.1E-20
+8.64E+01 5.1E-20
+1.01E+02 5.032E-20
+1.17E+02 4.906E-20
+1.36E+02 4.732E-20
+1.58E+02 4.521E-20
+1.84E+02 4.283E-20
+2.14E+02 4.029E-20
+2.49E+02 3.764E-20
+2.90E+02 3.497E-20
+3.37E+02 3.233E-20
+3.92E+02 2.975E-20
+4.56E+02 2.726E-20
+5.31E+02 2.489E-20
+6.17E+02 2.266E-20
+7.18E+02 2.056E-20
+8.35E+02 1.861E-20
+9.72E+02 1.68E-20
+1.13E+03 1.514E-20
+1.32E+03 1.361E-20
+1.53E+03 1.221E-20
+1.78E+03 1.094E-20
+2.07E+03 9.781E-21
+2.41E+03 8.735E-21
+2.80E+03 7.789E-21
+3.26E+03 6.938E-21
+3.79E+03 6.171E-21
+4.41E+03 5.484E-21
+5.13E+03 4.868E-21
+5.97E+03 4.316E-21
+6.94E+03 3.824E-21
+8.07E+03 3.385E-21
+9.39E+03 2.994E-21
+1.09E+04 2.646E-21
+1.27E+04 2.336E-21
+1.48E+04 2.062E-21
+1.72E+04 1.818E-21
+2.00E+04 1.602E-21
diff --git a/src/modules/init/moduleInput.f90 b/src/modules/init/moduleInput.f90
index 5090942..18c6e65 100644
--- a/src/modules/init/moduleInput.f90
+++ b/src/modules/init/moduleInput.f90
@@ -1331,6 +1331,7 @@ MODULE moduleInput
       USE moduleCaseParam, ONLY: tauMin
       USE moduleMesh, ONLY: mesh
       USE moduleSpecies, ONLY: nSpecies
+      USE moduleRefParam, ONLY: ti_ref
       IMPLICIT NONE
 
       TYPE(json_file), INTENT(inout):: config
@@ -1344,8 +1345,11 @@ MODULE moduleInput
         CALL config%get('average.startTime', tStart, found)
 
         IF (found) THEN
-          tAverageStart = INT(tStart / tauMin)
+          tAverageStart = INT(tStart / ti_ref / tauMin)
 
+        ELSE
+          tAverageStart = 0
+          
         END IF
 
         ALLOCATE(averageScheme(1:mesh%numNodes))