From db90bb88d6eaadfef86e5ba4f9bd8d51ff35a1c4 Mon Sep 17 00:00:00 2001
From: Lionel GUEZ <guez@lmd.ipsl.fr>
Date: Fri, 5 Jul 2024 19:53:11 +0200
Subject: [PATCH] Do not use `convert_to_ind`
Not worth it for a single point.
---
Overlap/read_snapshot.f90 | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/Overlap/read_snapshot.f90 b/Overlap/read_snapshot.f90
index 5f23b3cc..11673e27 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:
--
GitLab