diff --git a/plasmaExpansion.f90 b/plasmaExpansion.f90
index 1877a88..513d2d1 100644
--- a/plasmaExpansion.f90
+++ b/plasmaExpansion.f90
@@ -191,7 +191,7 @@ module output
       if (.not. res) then
         write (*, '(A, A)') 'Writing: ', filename
         open(unit=dataBC_id, file=pathOutput // filename, action='write', position='append')
-        write(dataBC_id, '(A,4(' // formatSep // ',A))') 't (s)', 'n (m^-3)', 'u (m s^-2)', 'T (eV)', 'Z'
+        write(dataBC_id, '(A,4(' // formatSep // ',A))') 't (s)', 'n (m^-3)', 'u (m s^-1)', 'T (eV)', 'Z'
         close(dataBC_id)
 
       end if