From a68c2ecba600c3b20211db03264092de35f64f5a Mon Sep 17 00:00:00 2001
From: Lionel GUEZ <guez@lmd.ipsl.fr>
Date: Thu, 4 Jul 2024 16:19:09 +0200
Subject: [PATCH] Polish
---
Inst_eddies/get_1_outerm.f90 | 3 ++-
Inst_eddies/set_contours.f90 | 1 +
2 files changed, 3 insertions(+), 1 deletion(-)
diff --git a/Inst_eddies/get_1_outerm.f90 b/Inst_eddies/get_1_outerm.f90
index b60a9b00..ce73601d 100644
--- a/Inst_eddies/get_1_outerm.f90
+++ b/Inst_eddies/get_1_outerm.f90
@@ -29,7 +29,8 @@ contains
! abs(ediff1d(cont_list_proj(:n_cont)%ssh)) is in descending
! order. (See notes for a proof.) On return, n_cont == 1 means we
! have only found one good contour at innermost_level_2 and it has
- ! sufficient area.
+ ! sufficient area. On return, n_cont == 0 .eqv. out_cont ==
+ ! null_polyline().
! Method: dichotomy on level of ssh.
diff --git a/Inst_eddies/set_contours.f90 b/Inst_eddies/set_contours.f90
index 2d12b75e..d7477c0d 100644
--- a/Inst_eddies/set_contours.f90
+++ b/Inst_eddies/set_contours.f90
@@ -100,6 +100,7 @@ contains
! Done with outermost contour, now let us take care of the
! max-speed contour.
if (out_cont%closed) then
+ ! assert n_cont >= 1
! {begin} Restrict the field to the outermost contour:
llc = floor((corner_whole - corner) / step &
--
GitLab