diff --git a/scripts_python/plotCumF.py b/scripts_python/plotCumF.py
index cc8a0a3..f04ccfb 100644
--- a/scripts_python/plotCumF.py
+++ b/scripts_python/plotCumF.py
@@ -9,22 +9,23 @@ from scipy.constants import e, k
 m_i = 1.9712e-25
 
 #  paths = ['../quasiNeutral_fullAblation/','../Poisson_fullAblation/']
-#  paths = ['../quasiNeutral_partialAblation/','../Poisson_partialAblation/']
-paths = ['../2024-09-27_17.41.21/']
-for path in paths:
+#  paths = ['../quasiNeutral_partialAblation/','../quasiNeutral_partialAblation_cs/','../Poisson_partialAblation/']
+paths = ['../2024-10-04_12.10.38/']
+labels = [path[3:-1] for path in paths]
+
+for path, label in zip(paths, labels):
     filesCum_i = sorted(glob.glob(path+'time_*_fCum_i.csv'))
     start = 0
     end   = len(filesCum_i)
     every = 100
 
-    for fileCum_i in filesCum_i[start:end:every]:
+    for fileCum_i in filesCum_i[start:end+1:every]:
         time, x, v, f_i = readF.read(fileCum_i)
 
         plt.plot(v**2*m_i*0.5/e, f_i[0]*e/m_i/v, label='{:.3f} ns'.format(time*1e9))
 
     time, x, v, f_i = readF.read(filesCum_i[-1])
-
-    plt.plot(v**2*m_i*0.5/e, f_i[0]*e/m_i/v, label='{:.3f} ns'.format(time*1e9), color='k')
+    plt.plot(v**2*m_i*0.5/e, f_i[0]*e/m_i/v, label='time_max', color='k')
 
     plt.yscale('log')
     plt.ylim([1e18,1e24])
@@ -35,7 +36,7 @@ for path in paths:
 
     plt.legend()
 
-    plt.title('r = {:.1f} mm, time_max={:.1f} ns '.format(x[0]*1e3, time*1e9) + path)
+    plt.title('r = {:.1f} mm, time_max={:.1f} ns, '.format(x[0]*1e3, time*1e9) + label)
 
     plt.show()
  
diff --git a/scripts_python/plotMom.py b/scripts_python/plotMom.py
index ce5a45c..f46e3dd 100644
--- a/scripts_python/plotMom.py
+++ b/scripts_python/plotMom.py
@@ -9,31 +9,33 @@ from scipy.constants import e, k
 #  paths = ['../quasiNeutral_fullAblation/','../2024-09-26_11.48.04/']
 #  paths = ['../2024-09-26_12.47.11/']
 #  paths = ['../quasiNeutral_partialAblation/','../2024-09-26_13.58.24/']
-#  path = '../quasiNeutral_fullAblation/'
-#  path = '../quasiNeutral_partialAblatio/'
-paths = ['../2024-09-27_17.41.21/']
+#  paths = ['../2024-10-02_14.30.44/']
+#  paths = ['../quasiNeutral_fullAblation/','../Poisson_fullAblation/']
+paths = ['../2024-10-03_19.28.22/']
+labels = [path[3:-1] for path in paths]
 
-for path in paths:
+for path, label in zip(paths, labels):
     filesPhi   = sorted(glob.glob(path+'time_*_phi.csv'))
     filesMom_i = sorted(glob.glob(path+'time_*_mom_i.csv'))
-    start = 0
+    start = 100
     end   = len(filesMom_i)
-    every = 50
+    every = 100
     fig, ax = plt.subplots(4, sharex='all')
     for fileMom_i, filePhi in zip(filesMom_i[start:end+1:every], filesPhi[start:end+1:every]):
         time, r, phi, E, n_e         = readPhi.read(filePhi)
         time, r, n_i, u_i, T_i, Zave = readMom.read(fileMom_i)
 
-        #  ax[0].plot(r, phi, label='t = {:.3f} ns'.format(time*1e9))
-        ax[0].plot(r, E, label='t = {:.3f} ns'.format(time*1e9))
+        ax[0].plot(r, phi, label='t = {:.1f} ns'.format(time*1e9))
+        #  ax[0].plot(r, E, label='t = {:.1f} ns'.format(time*1e9))
+        #  ax[0].plot(r, (Zave*n_i - n_e)/n_e, label='t = {:.1f} ns'.format(time*1e9))
         ax[1].set_yscale('log')
-        ax[1].set_ylim([1e14,2e25])
+        ax[1].set_ylim([1e20,2e25])
         ax[1].plot(r, Zave*n_i)
         ax[1].plot(r,      n_e, color='k', linestyle='dashed')
         ax[2].plot(r,      u_i)
         ax[3].plot(r,      T_i)
 
-    ax[0].set_title(path)
+    ax[0].set_title(label)
     ax[0].legend()
 
     plt.show()