Redefine `k_test_[12]` as date indices

Redefine `k_test_[12]` as date indices (0-based) instead of dates.

Overlap/Tests/test_overlap.f90

@@ -22,7 +22,7 @@ program test_overlap
 
   character(len = :), allocatable:: shpc_dir
   integer d_init
-  integer:: k_test_1 = 20454, k_test_2 = 20455
+  integer:: k_test_1 = 0, k_test_2 = 1
   integer unit, i, copy, rank, n_proc
   integer, allocatable:: ishape_last(:)
   type(snapshot), allocatable:: flow(:) ! (max_delta + 1)
@@ -81,10 +81,10 @@ program test_overlap
   call read_column(ishape_last, file = trim(shpc_dir) // "/ishape_last.txt")
   call shpc_open(hshp, trim(shpc_dir), rank = 0)
   call dbf_read_attribute_03(d_init, hshp%extremum, hshp%extr_date, ishape = 0)
-  call read_snapshot(flow(1), hshp, nlon, nlat, d_init, k_test_1 - d_init, &
+  call read_snapshot(flow(1), hshp, nlon, nlat, d_init, k_test_1, corner, &
+       step, copy, ishape_last)
+  call read_snapshot(flow(max_delta + 1), hshp, nlon, nlat, d_init, k_test_2, &
        corner, step, copy, ishape_last)
-  call read_snapshot(flow(max_delta + 1), hshp, nlon, nlat, d_init, &
-       k_test_2 - d_init, corner, step, copy, ishape_last)
   CALL shpc_close(hshp)
   print *, "Enter flow(1)%list_vis%delta_out (array with ", &
        flow(1)%number_vis_extr, " values):"
@@ -103,8 +103,8 @@ program test_overlap
 
   call shpc_create(hshp, shpc_dir = "SHPC", cyclone = hshp%cyclone)
   call init_interpolated_eddy
-  call overlap(flow, nlon, nlat, periodic, dist_lim, hshp, d = k_test_2, &
-       delta = max_delta, j = max_delta + 1)
+  call overlap(flow, nlon, nlat, periodic, dist_lim, hshp, &
+       d = d_init + k_test_2, delta = max_delta, j = max_delta + 1)
   close(unit_edge)
   print *, 'Created file "edgelist.csv".'
   CALL shpc_close(hshp)