Correct date subscript in Tests/Input/Region_4.

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

Tests/Input/Region_4/extraction_eddies_stdout.txt

+ &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".

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)

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.