diff --git a/src/opt_disl.f90 b/src/opt_disl.f90
index 9450171..4281d12 100644
--- a/src/opt_disl.f90
+++ b/src/opt_disl.f90
@@ -285,18 +285,18 @@ module opt_disl
         if(loop_radius < 0.0_dp) then 
                 ALLOCATE(xLoop(4,3))
                 xLoop(:,:) = 0.d0
-                xLoop(1,a1) = centroid(1) - loop_radius
-                xLoop(1,a2) = centroid(2) - loop_radius
-                xLoop(1,a3) = centroid(3)
-                 xLoop(2,a1) = centroid(1) + loop_radius
-                 xLoop(2,a2) = centroid(2) - loop_radius
-                 xLoop(2,a3) = centroid(3)
-                xLoop(3,a1) = centroid(1) + loop_radius
-                xLoop(3,a2) = centroid(2) + loop_radius
-                xLoop(3,a3) = centroid(3)
-                 xLoop(4,a1) = centroid(1) - loop_radius
-                 xLoop(4,a2) = centroid(2) + loop_radius
-                 xLoop(4,a3) = centroid(3)
+                xLoop(1,a1) = centroid(a1) + loop_radius
+                xLoop(1,a2) = centroid(a2) + loop_radius
+                xLoop(1,a3) = centroid(a3)
+                 xLoop(2,a1) = centroid(a1) - loop_radius
+                 xLoop(2,a2) = centroid(a2) + loop_radius
+                 xLoop(2,a3) = centroid(a3)
+                xLoop(3,a1) = centroid(a1) - loop_radius
+                xLoop(3,a2) = centroid(a2) - loop_radius
+                xLoop(3,a3) = centroid(a3)
+                 xLoop(4,a1) = centroid(a1) + loop_radius
+                 xLoop(4,a2) = centroid(a2) - loop_radius
+                 xLoop(4,a3) = centroid(a3)
         else
             !Calculate loop perimeter
             perimeter = 2.0_dp*pi*loop_radius