diff --git a/Overlap/read_snapshot.f90 b/Overlap/read_snapshot.f90 index 5f23b3cc7e31227c6f427af7faaa37eacfc040e5..11673e278a03c51ac0e2ba7412538b780ec941c7 100644 --- a/Overlap/read_snapshot.f90 +++ b/Overlap/read_snapshot.f90 @@ -7,7 +7,6 @@ contains subroutine read_snapshot(s, hshpc, ssm, k, copy) ! Libraries: - use contour_531, only: convert_to_ind use jumble, only: assert use numer_rec_95, only: hunt @@ -69,7 +68,7 @@ contains if (s%number_extr /= 0) then if (s%list(1)%extr%coord_proj(1) == huge(0)) & forall (eddy_i = 1:s%number_extr) s%list(eddy_i)%extr%coord_proj & - = nint(convert_to_ind(s%list(eddy_i)%extr%coord, corner, step)) + = nint((s%list(eddy_i)%extr%coord - corner) / step + 1.) end if ! Define s%extr_map from s%list%extr%coord_proj: