Missing allocation in disloop code

5 years ago · 94e1c9fd7b
parent 071e8169b9
commit 94e1c9fd7b
1 changed files with 1 additions and 0 deletions
--- a/src/opt_disl.f90
+++ b/src/opt_disl.f90
@ -292,6 +292,7 @@ module opt_disl
        end select

        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