diff --git a/moduleOutput.f90 b/moduleOutput.f90
index 74f41f9..a463d77 100644
--- a/moduleOutput.f90
+++ b/moduleOutput.f90
@@ -222,6 +222,24 @@ module output
 
     end subroutine writeOutputTime
 
+    subroutine writeOutputZList(nz, Z_list)
+      integer,  intent(in):: nz
+      real(dp), intent(in):: Z_list(1:nz)
+      character(:), allocatable:: filename
+      integer:: i
+
+      filename = 'ZList.csv'
+      write (*, '(A, A)') 'Writing: ', filename
+      open(unit=dataPhi_id, file=pathOutput//filename)
+      write(dataPhi_id, '(A)') "Z_list"
+      do i = 1, nz
+        write(dataPhi_id, '('//formatFloat//')') Z_list(i)
+
+      end do
+      close(unit=dataPhi_id)
+
+    end subroutine writeOutputZList
+
     subroutine writeOutputRef()
       use referenceValues,    only: t_ref, L_ref, n_ref, u_ref, Temp_ref, phi_ref
       use constantParameters, only: eV_to_K
diff --git a/scripts_python/readZlist.py b/scripts_python/readZlist.py
new file mode 100644
index 0000000..68488f7
--- /dev/null
+++ b/scripts_python/readZlist.py
@@ -0,0 +1,10 @@
+import pandas
+
+def read(filename):
+    # Get time
+    df = pandas.read_csv(filename)
+    Zlist = df['Z_list'].to_numpy()
+
+    return Zlist
+
+