diff --git a/src/dcore/impurity_solvers/hphi_spectrum.py b/src/dcore/impurity_solvers/hphi_spectrum.py
index 38c65637..86b0d9cc 100644
--- a/src/dcore/impurity_solvers/hphi_spectrum.py
+++ b/src/dcore/impurity_solvers/hphi_spectrum.py
@@ -124,6 +124,14 @@ def calc_one_body_green(one_body_green_core):
             G[i, j] = G_ij[i, j]
             G[j, i] = G_ji[i, j]
 
+    # for i in range(G_shape[0]):
+    #     for j in range(i+1, G_shape[1]):  # i<j
+    #         g_iijj = G_diag[i] - G_diag[j]
+    #         A2 = A[i, j] - g_iijj
+    #         B2 = B[i, j] - g_iijj
+    #         G[i, j] = (A2 - 1J * B2) / 2.
+    #         G[j, i] = (A2 + 1J * B2) / 2.
+
     return G.reshape(n_site, n_sigma, n_site, n_sigma, n_T, n_omega)
 
 
diff --git a/src/dcore/impurity_solvers/pomerol.py b/src/dcore/impurity_solvers/pomerol.py
index eb27172c..f9cc3fc2 100644
--- a/src/dcore/impurity_solvers/pomerol.py
+++ b/src/dcore/impurity_solvers/pomerol.py
@@ -83,14 +83,16 @@ def assign_from_numpy_array(g_block, data, block_names):
     for b in range(data.shape[0]):
         for o1 in range(data.shape[1]):
             for o2 in range(data.shape[2]):
-                ax.plot(data[b, o1, o2, :30].real, marker='.')
+                ax.plot(data[b, o1, o2, :30].real, marker='.', label=f"{b}, ({o1},{o2})")
+    ax.legend()
     fig.savefig("test_gf_re.pdf")
 
     fig, ax = plt.subplots()
     for b in range(data.shape[0]):
         for o1 in range(data.shape[1]):
             for o2 in range(data.shape[2]):
-                ax.plot(data[b, o1, o2, :30].imag, marker='.')
+                ax.plot(data[b, o1, o2, :30].imag, marker='.', label=f"{b}, ({o1},{o2})")
+    ax.legend()
     fig.savefig("test_gf_im.pdf")