From 62a9f4915c8359c2ec1c942bd566cf5770d999fc Mon Sep 17 00:00:00 2001
From: Lionel GUEZ <guez@lmd.ens.fr>
Date: Thu, 4 Apr 2019 16:32:48 +0200
Subject: [PATCH] Correct date subscript in Tests/Input/Region_4.

---
 Documentation_texfol/documentation.tex        |  29 +++++++++++-------
 .../Region_4/extraction_eddies_stdout.txt     |  16 ++++++++++
 Tests/Input/Region_4/extremum.dbf             | Bin 452 -> 458 bytes
 Tests/Input/Region_4/max_speed_contour.dbf    | Bin 326 -> 332 bytes
 Tests/Input/Region_4/outermost_contour.dbf    | Bin 368 -> 374 bytes
 Tests/test_successive_overlap.f90             |   6 ++--
 derived_types.f90                             |   6 ++--
 7 files changed, 40 insertions(+), 17 deletions(-)
 create mode 100644 Tests/Input/Region_4/extraction_eddies_stdout.txt

diff --git a/Documentation_texfol/documentation.tex b/Documentation_texfol/documentation.tex
index de094c8c..f1a35156 100644
--- a/Documentation_texfol/documentation.tex
+++ b/Documentation_texfol/documentation.tex
@@ -59,7 +59,9 @@ de la vitesse.
 But du programme : produire une base de données de tourbillons. Le
 programme ne traite qu'une date.
 
-Cas pathologiques. \verb+out_cont+ est vide si et seulement s'il n'y a
+\subsection{Cas pathologiques}
+
+\verb+out_cont+ est vide si et seulement s'il n'y a
 pas de bon contour à \verb+innermost_level+. Le champ speed de
 extremum.dbf vaut \verb+missing_speed+ si et seulement si radius4 == 0
 ou \verb+speed_outerm+ dans \verb+set_max_speed+ est NaN. Si radius4
@@ -76,8 +78,9 @@ speed dans extremum.dbf vaut \verb+missing_speed+).
 
 Cas où les composantes \verb+out_cont+ ou \verb+speed_cont+ d'un
 tourbillon sont vides, c'est-à-dire égales à
-\verb+null_ssh_contour()+. On peut distinguer les cas du tableau
-(\ref{tab:null_ssh_contour}).
+\verb+null_ssh_contour()+. Si \verb+out_cont+ est vide alors
+\verb+speed_cont+ est vide aussi. On peut distinguer les cas du
+tableau (\ref{tab:null_ssh_contour}).
 \begin{table}[htbp]
   \centering
   \begin{tabular}{lllll|l}
@@ -91,13 +94,12 @@ tourbillon sont vides, c'est-à-dire égales à
     $\ne$ null & T & $\ge 2$ & $\ne \np{e4}$, $\ne$ NaN & 21 \%
                                                                  & 3
   \end{tabular}
-  \caption{Cas où les composantes out\_cont ou speed\_cont d'un tourbillon sont
-    vides. speed\_cont est nul dans tous ces cas. Colonne fraction :
-    pour la région 5, nombre de cas sur 260 extremums au total. Note 1
-    : pas de bon contour à innermost\_level. Note 2 : out\_cont est
-    près d'une côte, champ de vitesse non défini, max\_speed vaut
-    $10^4$ dans le fichier DBF. Note 3 : pas mieux que out\_cont pour
-    la vitesse maximale.}
+  \caption{Cas où la composante speed\_cont d'un tourbillon est
+    vide. Colonne fraction : pour la région 5, nombre de cas sur 260
+    extremums au total. Note 1 : pas de bon contour à
+    innermost\_level. Note 2 : out\_cont est près d'une côte, champ de
+    vitesse non défini, max\_speed vaut $10^4$ dans le fichier
+    DBF. Note 3 : pas mieux que out\_cont pour la vitesse maximale.}
   \label{tab:null_ssh_contour}
 \end{table}
 Sur la région 5, sur 260 extremums, je trouve les nombres de cas
@@ -106,6 +108,9 @@ suivants pour les 5 lignes du tableau (\ref{tab:null_ssh_contour}) :
 \verb+speed_cont+ d'un tourbillon sont vides, sur 260 tourbillons, soit 81
 \% des tourbillons.
 
+radius4 == 0 si et seulement si non valid. radius4 $\ge 1$ si et
+seulement si valid.
+
 \subsection{Calcul de la vitesse azimutale par rapport à l'extremum}
 
 Notons $C$ la position de l'extremum et $M$ un point où l'on veut
@@ -333,7 +338,9 @@ shx.
 
 valid == 1 dans le fichier
 \verb+extremum_$m.dbf+ signifie qu'un contour extérieur a été trouvé
-et qu'il a une surface suffisante.
+et qu'il a une surface suffisante. Si un contour extérieur est bien
+défini (non vide) alors il a par construction une amplitude supérieure
+à \verb+min_amp+.
 
 Si une forme dans \verb+max_speed_contour_$m.shp+ est non vide alors
 la vitesse correspondante dans \verb+extremum_$m.dbf+ est la vitesse
diff --git a/Tests/Input/Region_4/extraction_eddies_stdout.txt b/Tests/Input/Region_4/extraction_eddies_stdout.txt
new file mode 100644
index 00000000..9bab15fc
--- /dev/null
+++ b/Tests/Input/Region_4/extraction_eddies_stdout.txt
@@ -0,0 +1,16 @@
+ &MAIN_NML
+ MIN_AMP = 1.0000000E-03  ,
+ MAX_RADIUS      = 2*20,
+ MIN_RADIUS      = 25.00000       
+ /
+ Enter namelist main_nml.
+ &MAIN_NML
+ MIN_AMP = 1.0000000E-03  ,
+ MAX_RADIUS      = 2*20,
+ MIN_RADIUS      = 25.00000       
+ /
+ periodic =  F
+ Number of extrema: 5
+ Created shapefile "extremum_1".
+ Created shapefile "outermost_contour_1".
+ Created shapefile "max_speed_contour_1".
diff --git a/Tests/Input/Region_4/extremum.dbf b/Tests/Input/Region_4/extremum.dbf
index 579cfe2987ed5bdcafbb0970e85003581b853deb..f272e0f617342dbd4da5048cd1f222fc1ba9f349 100644
GIT binary patch
delta 92
zcmX@Ye2STyIi6jLm4ShQkx^qJw<{y-#K3*5Mg}IPCX<C3m7y#vI4cp(>R?os2CGv5
T0#gG+0|i3^AkYI51_n|9Ndpno

delta 86
zcmX@be1w^sIi6jLm4ShQkx_jjw<{yd#K3(l3JMB_lSLSnAq*QBBL&9jVN{j`sR3dI
SQv*XVt!toXU}0=vWDEeUL=STS

diff --git a/Tests/Input/Region_4/max_speed_contour.dbf b/Tests/Input/Region_4/max_speed_contour.dbf
index 9ef3accdd7e36c10b06042935fd5d489cfd44354..3d79cb09dfd4aacaac01c1fb912e5ee81aca34ba 100644
GIT binary patch
delta 73
zcmX@cbcTtWIi6jLm4Sg_A%oIH?o>wBiG?1lMg}IPCKLC{LK*j9jQ>!^WCcc9Zm=Q+
IATX5z0Cx%yoB#j-

delta 76
zcmX@Zbc~6cIi6jLm4Sg_A%o&X?o>vWiG?043JMB_6ZgqN=m)Z#AU+T)7*AX%0q2+?
Gb4&qxL=WBo

diff --git a/Tests/Input/Region_4/outermost_contour.dbf b/Tests/Input/Region_4/outermost_contour.dbf
index ed036b0bdab4a325fb5c8a2bd9bc4fcc61593a2f..717bb7966da3e096c0bc5119b3f8db8d4328a906 100644
GIT binary patch
delta 102
zcmeys^o@y|Ii6jLm4SibAcM+8?o>wBiG>-gMg}IPCKK<<Ll~17<u$<~3P512U<9U2
h^vo>HO+icr13e>46C+a-u$swkjPks2)usvtQUG_r6q5h|

delta 97
zcmeyy^nr<+Ii6jLm4SibAcOKm?o>vWiG>*~3JMB_6Yt5ZgJ>XDFakm#Ff-OOGc+~;
jF%=B-jEqe!42^+I!^vul@?bray%^=WK_Va(rV54t+948T

diff --git a/Tests/test_successive_overlap.f90 b/Tests/test_successive_overlap.f90
index 3c3c4709..276bdb67 100644
--- a/Tests/test_successive_overlap.f90
+++ b/Tests/test_successive_overlap.f90
@@ -60,9 +60,9 @@ program test_successive_overlap
        action  = "write")
 
   ! Title lines:
-  write(unit_edgelist, fmt = "(1x, a)") '"predecessor date index" ' &
-       // '"predecessor eddy index" "successor date index" ' &
-       // '"successor eddy index"'
+  write(unit_edgelist, fmt = "(1x, a)") '"predecessor date subscript" ' &
+       // '"predecessor eddy subscript" "successor date subscript" ' &
+       // '"successor eddy subscript"'
   write(unit_edgelist, fmt = *) "k1 i1 k2 i2 weight"
   
   call successive_overlap(unit_edgelist, nlon, nlat, flow, j = 2, k = 2)
diff --git a/derived_types.f90 b/derived_types.f90
index 9170cb4b..c82e669a 100644
--- a/derived_types.f90
+++ b/derived_types.f90
@@ -25,11 +25,11 @@ module derived_types
      logical interpolated
      
      integer:: delta_in = huge(0)
-     ! Maximum difference in time index where there is a direct
+     ! Maximum difference in time subscript where there is a direct
      ! predecessor. huge means there is no direct predecessor.
      
      integer:: delta_out = huge(0)
-     ! Maximum difference in time index where there is a direct
+     ! Maximum difference in time subscript where there is a direct
      ! successor. huge means there is no direct succcessor.
 
      integer radius4
@@ -52,7 +52,7 @@ module derived_types
      ! identification number or this extremum. 0 at other points.
 
      integer, allocatable:: ind_extr(:, :) ! (2, number_vis_extr)
-     ! List of coordinates of ssh extrema in the global grid, in index
+     ! List of coordinates of ssh extrema in the global grid, in subscript
      ! space. The position in this list is the identifying number of
      ! the extremum.
      
-- 
GitLab