diff --git a/figures/MAM2EBRAINS/M2E_compute_corrcoeff.py b/figures/MAM2EBRAINS/M2E_compute_corrcoeff.py
index 37826a7a24a3c2756a84a016681a37fd36648bd3..1b486296c04d90927b025043cbc08c16185da0c2 100644
--- a/figures/MAM2EBRAINS/M2E_compute_corrcoeff.py
+++ b/figures/MAM2EBRAINS/M2E_compute_corrcoeff.py
@@ -32,11 +32,8 @@ def compute_corrcoeff(M, data_path, label):
     resolution = 1.
 
     cc_dict = {}
-
     for area in areas_simulated:
         cc_dict[area] = {}
-        LvR_list = []
-        N = []
         for pop in M.structure[area]:
             fp = '-'.join((label,
                            'spikes',  # assumes that the default label for spike files was used
diff --git a/figures/MAM2EBRAINS/M2E_compute_pop_LvR.py b/figures/MAM2EBRAINS/M2E_compute_pop_LvR.py
index bd6ca730ff581c80b626245eeaae1059c59eadb1..0cb91a6e4e5fcc16416f08546415d851985060ed 100644
--- a/figures/MAM2EBRAINS/M2E_compute_pop_LvR.py
+++ b/figures/MAM2EBRAINS/M2E_compute_pop_LvR.py
@@ -32,8 +32,6 @@ def compute_pop_LvR(M, data_path, label):
 
     for area in areas_simulated:
         pop_LvR_dict[area] = {}
-        LvR_list = []
-        N = []
         for pop in M.structure[area]:
             fp = '-'.join((label,
                            'spikes',  # Assumes that the default label for spike files was used
diff --git a/figures/MAM2EBRAINS/M2E_compute_rate_time_series.py b/figures/MAM2EBRAINS/M2E_compute_rate_time_series.py
index ff594137c775ea7f036891bf4c9bbf56e6930bad..ce83fd1b9dc46edfe01d732c88b3f93d10dae32d 100644
--- a/figures/MAM2EBRAINS/M2E_compute_rate_time_series.py
+++ b/figures/MAM2EBRAINS/M2E_compute_rate_time_series.py
@@ -52,7 +52,7 @@ def compute_rate_time_series(M, data_path, label, area, method):
     N_list = []
     for pop in M.structure[area]:
         fp = '-'.join((label,
-                       'spikes',  # assumes that the default label for spike files was used
+                       'spikes',  # Assumes that the default label for spike files was used
                        area,
                        pop))
         fn = '{}/{}.npy'.format(load_path, fp)
diff --git a/figures/MAM2EBRAINS/M2E_firing_rate.py b/figures/MAM2EBRAINS/M2E_firing_rate.py
index a59a332d131adf772e95a3bff02f010282e61d80..c89214eab201aaccd8a6e35641db117de711d2be 100644
--- a/figures/MAM2EBRAINS/M2E_firing_rate.py
+++ b/figures/MAM2EBRAINS/M2E_firing_rate.py
@@ -73,7 +73,6 @@ def plot_firing_rate_over_areas(M, data_path):
 
     t_min = 500.
     t_max = M.simulation.params['t_sim']
-    time = np.arange(t_min, t_max)
     matrix = []
     for area in area_list:
         binned_spikes = rate_time_series[area][500:600]
diff --git a/figures/MAM2EBRAINS/M2E_visualize_fc.py b/figures/MAM2EBRAINS/M2E_visualize_fc.py
index ca9b60c43626b5185554d28c03ab873f2e159b0c..749144ed1a5b08610f9c30d27dd481017b9a26fa 100644
--- a/figures/MAM2EBRAINS/M2E_visualize_fc.py
+++ b/figures/MAM2EBRAINS/M2E_visualize_fc.py
@@ -90,7 +90,6 @@ def visualize_fc(M, data_path):
     Returns:
         None
     """
-    area_list = M.simulation.params["areas_simulated"]
     label = M.simulation.label
     
     # Compute functional connectivity