Fix bug in neighbors

3 years ago · 044bfe7e7f
parent d5d20e468c
commit 044bfe7e7f
1 changed files with 2 additions and 0 deletions
--- a/src/neighbors.f90
+++ b/src/neighbors.f90
@ -94,6 +94,8 @@ module neighbors
        end if

        allocate(nei_list(100,n),nei_num(n))
+        nei_list=0
+        nei_num=0

        !Now first pass the position list and and point num to the cell algorithm
        call build_cell_list(n, r_list, rc_off, cell_num, num_in_cell, cell_list, which_cell)