From 6a0d441bd072658e359b1c9c0048bfd3a79c737d Mon Sep 17 00:00:00 2001 From: Alex Date: Fri, 21 Feb 2020 17:18:39 -0500 Subject: [PATCH] Fix to disloop code --- src/opt_disl.f90 | 24 ++++++++++++------------ 1 file changed, 12 insertions(+), 12 deletions(-) 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