diff --git a/swash/config.ini b/swash/config.ini
index d9adee9..036fd96 100644
--- a/swash/config.ini
+++ b/swash/config.ini
@@ -9,7 +9,7 @@ mpi=8
 
 [post]
 inp=inp_post/spec_2lay
-compare=inp_post/spec_2lay
-out=out_post/spec_2lay
+compare=inp_post/spec_nolay
+out=out_post/spec_2layvsnolay
 x0=-1250
 t0=180
diff --git a/swash/processing/post.py b/swash/processing/post.py
index 8710133..b8e6d17 100644
--- a/swash/processing/post.py
+++ b/swash/processing/post.py
@@ -59,10 +59,10 @@ R = np.sqrt(
     (np.abs(phi_eta[1]) + np.abs(phi_u[1]) - 2 * np.abs(phi_eta_u[1]))
     / (np.abs(phi_eta[1]) + np.abs(phi_u[1]) + 2 * np.abs(phi_eta_u[1]))
 )
-# R = np.sqrt(
-#    (1 + G**2 - 2 * G * np.cos(th_eta_u))
-#    / (1 + G**2 + 2 * G * np.cos(th_eta_u))
-# )
+#R = np.sqrt(
+#   (1 + G**2 - 2 * G * np.cos(th_eta_u))
+#   / (1 + G**2 + 2 * G * np.cos(th_eta_u))
+#)
 
 if config.has_option("post", "compare"):
     inp_comp = pathlib.Path(config.get("post", "compare"))
@@ -88,14 +88,14 @@ if config.has_option("post", "compare"):
     G_ = H_ / U_
     th_eta_u_ = np.angle(phi_eta_u_[1])
 
-    #R_ = np.sqrt(
-    #    (np.abs(phi_eta_[1]) + np.abs(phi_u_[1]) - 2 * np.abs(phi_eta_u_[1]))
-    #    / (np.abs(phi_eta_[1]) + np.abs(phi_u_[1]) + 2 * np.abs(phi_eta_u_[1]))
-    #)
     R_ = np.sqrt(
-       (1 + G_**2 - 2 * G_ * np.cos(th_eta_u_))
-       / (1 + G_**2 + 2 * G_ * np.cos(th_eta_u_))
+        (np.abs(phi_eta_[1]) + np.abs(phi_u_[1]) - 2 * np.abs(phi_eta_u_[1]))
+        / (np.abs(phi_eta_[1]) + np.abs(phi_u_[1]) + 2 * np.abs(phi_eta_u_[1]))
     )
+    #R_ = np.sqrt(
+    #   (1 + G_**2 - 2 * G_ * np.cos(th_eta_u_))
+    #   / (1 + G_**2 + 2 * G_ * np.cos(th_eta_u_))
+    #)
 
 
 # Plotting