From ed2abbd100dafc25871ccf0516e864baac15989e Mon Sep 17 00:00:00 2001
From: Lionel GUEZ <guez@lmd.ens.fr>
Date: Mon, 22 Nov 2021 16:33:21 +0100
Subject: [PATCH] Shortcut if `radius <= 2`
---
Inst_eddies/max_speed_contour_ssh.f90 | 7 ++++++-
1 file changed, 6 insertions(+), 1 deletion(-)
diff --git a/Inst_eddies/max_speed_contour_ssh.f90 b/Inst_eddies/max_speed_contour_ssh.f90
index 7cb7e84f..8b70fcef 100644
--- a/Inst_eddies/max_speed_contour_ssh.f90
+++ b/Inst_eddies/max_speed_contour_ssh.f90
@@ -50,7 +50,12 @@ contains
! outermost contour as max-speed contour.
max_speed_contour_ssh = missing_ssh
else
- l = maxloc(abs(sum(v_azim, dim = 1) / 4.), dim = 1)
+ if (radius >= 3) then
+ l = maxloc(abs(sum(v_azim, dim = 1) / 4.), dim = 1)
+ else
+ l = 1
+ end if
+
direction = maxloc(abs(v_azim(:, l)), dim = 1)
select case (direction)
--
GitLab