Fixed bad spill operation on the symba_kinship type

MintonGroup · Oct 13, 2021 · aa37c3d · aa37c3d
1 parent c3ef158
commit aa37c3d
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/symba/symba_util.f90 b/src/symba/symba_util.f90
@@ -1102,6 +1102,7 @@ module subroutine symba_util_spill_arr_kin(keeps, discards, lspill_list, ldestru
       logical,                                        intent(in)    :: ldestructive !! Logical flag indicating whether or not this operation should alter the keeps array or not
       ! Internals
       integer(I4B) :: nspill, nkeep, nlist
+      type(symba_kinship), dimension(:), allocatable :: tmp
 
       nkeep = count(.not.lspill_list(:))
       nspill = count(lspill_list(:))
@@ -1118,7 +1119,9 @@ module subroutine symba_util_spill_arr_kin(keeps, discards, lspill_list, ldestru
       discards(:) = pack(keeps(1:nlist), lspill_list(1:nlist))
       if (ldestructive) then
          if (nkeep > 0) then
-            keeps(:) = pack(keeps(1:nlist), .not. lspill_list(1:nlist))
+            allocate(tmp(nkeep))
+            tmp(:) = pack(keeps(1:nlist), .not. lspill_list(1:nlist))
+            call move_alloc(tmp, keeps)
          else
             deallocate(keeps)
          end if