diff --git a/Inst_eddies/Analysis/distribution_function.py b/Inst_eddies/Analysis/distribution_function.py
index 4d52898715c8832d268653ae1f1ae3692f26e38b..bb384fb77cbe674e3c153d6776f67e8d611051a6 100755
--- a/Inst_eddies/Analysis/distribution_function.py
+++ b/Inst_eddies/Analysis/distribution_function.py
@@ -129,7 +129,8 @@ def plot_all(dict_list, label = None):
     
 if __name__ == "__main__":
     import sys
-    
+
+    if len(sys.argv) != 2: sys.exit("Required argument: directory")
     d = read(sys.argv[1])
     plot_all([d])
     plt.show()
diff --git a/Inst_eddies/Documentation_texfol/documentation.tex b/Inst_eddies/Documentation_texfol/documentation.tex
index 01b521f944b423c51f0a5c40e03762f9e3f5c381..e309088374293a712131cdd49a5ba3e60ef2e024 100644
--- a/Inst_eddies/Documentation_texfol/documentation.tex
+++ b/Inst_eddies/Documentation_texfol/documentation.tex
@@ -16,7 +16,7 @@
 
 \usepackage{hyperref}
 
-\hypersetup{pdftitle={Documentation}, pdfauthor={Lionel
+\hypersetup{pdftitle={Instantaneous eddies}, pdfauthor={Lionel
     Guez}, hypertexnames=false, pdfborderstyle={/S/U/W 1}}
 
 \usepackage{algorithm}
@@ -837,23 +837,22 @@ arguments.
 \section{Tests}
 
 Pour faire des tests : données au 1\ier{} janvier 2006. Tests sur
-différents domaines. Cf. figure (\ref{fig:regions}).
+différents domaines. Cf. figure (\ref{fig:regions}) et
+\verb+domains.ods+.
 \begin{figure}[htbp]
   \centering
   \includegraphics[width=\textwidth]{regions}
   \caption{Régions pour les tests.}
   \label{fig:regions}
 \end{figure}
-Régions 1, 2 et 3 : 29 $\times$ 17, 57 $\times$ 33; 120 $\times$ 120,
-contenant respectivement 8, 27 et 216 extremums.
 
 Pour la région 1, coin en :
 \begin{equation*}
   \lambda = \np{5.125}^\circ, \phi = -\np{36.375}^\circ
 \end{equation*}
 soit $(\np[rad]{0.0894}, - \np[rad]{0.6349})$. Indices base 1 du coin
-: 21, 215. Pas de grille : \np{0.25}°, soit \np[rad]{0.0044}. 5
-minimums et 3 maximums dans ce cadre. Cf. figure (\ref{fig:extrema}).
+: 21, 215. 5 minimums et 3 maximums dans ce cadre. Cf. figure
+(\ref{fig:extrema}).
 \begin{figure}[htbp]
   \centering
   \includegraphics[width=\textwidth]{extrema}
diff --git a/Overlap/Documentation_texfol/documentation.tex b/Overlap/Documentation_texfol/documentation.tex
index 9283f7d2efb437062c368828c30f9237d1cc924d..3ec054042f2fe1c8a421aceeed1120e3cefbca18 100644
--- a/Overlap/Documentation_texfol/documentation.tex
+++ b/Overlap/Documentation_texfol/documentation.tex
@@ -16,7 +16,7 @@
 
 \usepackage{hyperref}
 
-\hypersetup{pdftitle={Documentation}, pdfauthor={Lionel
+\hypersetup{pdftitle={Overlap}, pdfauthor={Lionel
     Guez}, hypertexnames=false, pdfborderstyle={/S/U/W 1}}
 
 \usepackage{algorithm}
@@ -166,9 +166,6 @@ Cf. discussion dans
   du programme de détection des tourbillons instantanés}, § \og Résumé
 du problème de géométrie algorithmique\fg{}.
 
-Compilation et exécution testées avec NAG Fortran Compiler Release
-6.2.
-
 Dans le programme de recouvrement, idée d'admettre des tourbillons
 dont les numéros à une date donnée n'iraient pas forcément de 1 au
 nombre de tourbillons visibles. L'intérêt serait de pouvoir traiter
@@ -950,4 +947,12 @@ intervenir les tourbillons interpolés le poids associé aux tourbillons
 visibles entre lesquels on interpole. Normalement, ce choix doit
 donner un poids plus faible à ces arcs.
 
+\section{Tests}
+
+Compilation et exécution testées avec NAG Fortran Compiler Release
+6.2.
+
+Révision ed60f10, ifort \verb+19.0+, sur ciclad, 1 processus MPI, test
+EGg. Temps d'exécution en mode debug : 34 s, en mode release : 9 s.
+
 \end{document}
diff --git a/Overlap/eddy_graph_in.sh b/Overlap/eddy_graph_in.sh
index 26298fc0efafde74b6a7d84436c7d37e7c0d8211..f3876bbf67a06c5d818555b1923fa2123e5a52c5 100755
--- a/Overlap/eddy_graph_in.sh
+++ b/Overlap/eddy_graph_in.sh
@@ -1,5 +1,8 @@
 #!/bin/bash
 
+# This program creates edgelists and interpolated eddies. You run it
+# only once and it processes both anticyclones and cyclones.
+
 # Note: no need to sort number_eddies.csv, already sorted by construction.
 
 if (($# != 2))
diff --git a/Overlap/write_edge.f90 b/Overlap/write_edge.f90
index 1093e25a25753f72ec45054bbb60f9df8972408b..69538262fec14b92220859229bc24056453a71ea 100644
--- a/Overlap/write_edge.f90
+++ b/Overlap/write_edge.f90
@@ -12,13 +12,13 @@ contains
     integer, intent(in):: k1, i1, k2, i2
     real, intent(in):: w
 
+    ! Local:
+    integer unit_edge
+
     !--------------------------------------------------------------------
 
-    if (cyclone) then
-       write(unit_edge_cyclo, fmt = *) k1, i1, k2, i2, w
-    else
-       write(unit_edge_anti, fmt = *) k1, i1, k2, i2, w
-    end if
+    unit_edge = merge(unit_edge_cyclo, unit_edge_anti, cyclone)
+    write(unit_edge, fmt = *) k1, i1, k2, i2, w
     
   end subroutine write_edge