Fixed bug in fragment position initialization that accidentally rando…

…mized the 2nd largest fragment instead of putting it at the opposite end
MintonGroup · May 28, 2021 · f15dc6e · f15dc6e
1 parent 7150acd
commit f15dc6e
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90
@@ -453,7 +453,7 @@ subroutine set_fragment_position_vectors()
       x_frag(:, 1) = -y_col_unit(:) * r_max  
       x_frag(:, 2) =  y_col_unit(:) * r_max  
 
-      call random_number(x_frag(:,2:nfrag))
+      call random_number(x_frag(:,3:nfrag))
       loverlap(:) = .true.
       do while (any(loverlap(3:nfrag)))
          do i = 3, nfrag