diff --git a/Inst_eddies/Documentation_texfol/Graphiques/nearby_extr.odg b/Inst_eddies/Documentation_texfol/Graphiques/nearby_extr.odg
index c93afc97b85c52dc73d2fe136469c114c69c0933..a070129abd19fbc16e3a2d43d068ee0a2e36e71d 100644
Binary files a/Inst_eddies/Documentation_texfol/Graphiques/nearby_extr.odg and b/Inst_eddies/Documentation_texfol/Graphiques/nearby_extr.odg differ
diff --git a/Inst_eddies/Documentation_texfol/documentation.tex b/Inst_eddies/Documentation_texfol/documentation.tex
index d5dde9ab0ce84d4e3f24b1ca3e51b152d308fda4..d6f70da365a70a84eea1c698fe1c99002070928b 100644
--- a/Inst_eddies/Documentation_texfol/documentation.tex
+++ b/Inst_eddies/Documentation_texfol/documentation.tex
@@ -1256,6 +1256,12 @@ puisqu'on a initialisé le champ closed à vrai. Cf. figure
   \label{fig:nearby_extr}
 \end{figure}
 
+Nota bene : les numéros d'identification des tourbillons stockés dans
+\verb+extr_map+ ne sont a priori pas les mêmes que les numéros
+intra-date stockés dans les shapefiles. Un numéro dans \verb+extr_map+
+identifie un tourbillon dans l'ensemble des tourbillons des deux
+orientations.
+
 \subsection{set\_contours}
 \label{sec:set_contours}
 
diff --git a/Inst_eddies/Tests/test_nearby_extr.f90 b/Inst_eddies/Tests/test_nearby_extr.f90
index 447e672d72d55500e9cd9696b58c97a875f1aa99..fb3f853550fe0734c4ff8841e13bde9eb8940ac6 100644
--- a/Inst_eddies/Tests/test_nearby_extr.f90
+++ b/Inst_eddies/Tests/test_nearby_extr.f90
@@ -1,7 +1,7 @@
 program test_nearby_extr
 
   ! Libraries:
-  use jumble, only: get_command_arg_dyn, read_opcol, rad_to_deg
+  use jumble, only: get_command_arg_dyn, read_opcol
   use shapelib_03, only: dbf_read_attribute_03
 
   use derived_types, only: snapshot, shpc_slice_handler, shpc_slice_meta
@@ -19,8 +19,8 @@ program test_nearby_extr
   integer l
   type(shpc_slice_meta) ssm
 
-  real, allocatable:: nearby(:, :) ! (2, :) longitude and
-  ! latitude, in rad, of extrema near the target extremum
+  integer, allocatable:: nearby(:, :) ! (2, :)
+  ! coordinates in projection space of extrema near the target extremum
 
   integer:: i = 3
   namelist /main_nml/ i
@@ -42,7 +42,7 @@ program test_nearby_extr
   print *, "nearby:"
 
   do l = 1, size(nearby, 2)
-     print *, nearby(:, l) * rad_to_deg
+     print *, nearby(:, l)
   end do
 
 end program test_nearby_extr
diff --git a/Inst_eddies/inst_eddies.f90 b/Inst_eddies/inst_eddies.f90
index 4f3c36d7d08c1206a2feef3aae1519309c99a78a..cc173a182bd4c0cfab4503f3e22cf8d6d456925f 100644
--- a/Inst_eddies/inst_eddies.f90
+++ b/Inst_eddies/inst_eddies.f90
@@ -7,7 +7,7 @@ program inst_eddies
   use, intrinsic:: ISO_FORTRAN_ENV, only: ERROR_UNIT
 
   ! Libraries:
-  use contour_531, only: convert_to_ind
+  use contour_531, only: convert_to_ind, convert_to_reg_coord
   use jumble, only: assert, new_unit, argwhere, twopi
   use numer_rec_95, only: indexx
 
@@ -160,8 +160,9 @@ program inst_eddies
 
      ! Done defining geographical window
 
-     outside_points = nearby_extr(s%extr_map(llc(1):urc(1), llc(2):urc(2)), &
-          s%list, i)
+     outside_points &
+          = convert_to_reg_coord(nearby_extr(s%extr_map(llc(1):urc(1), &
+          llc(2):urc(2)), s%list, i), corner, step)
 
      ! Shift the longitude of each point to a value close to the
      ! longitude of the target extremum:
diff --git a/Inst_eddies/nearby_extr.f90 b/Inst_eddies/nearby_extr.f90
index 052c3a0ef6b9c446bc31e165be4b9983c63d7cea..c028d1f7af338f78d4f10740801a5da1754147b8 100644
--- a/Inst_eddies/nearby_extr.f90
+++ b/Inst_eddies/nearby_extr.f90
@@ -11,8 +11,8 @@ contains
 
     use derived_types, only: eddy
 
-    real, allocatable:: nearby_extr(:, :) ! (2, :) longitude and
-    ! latitude, in rad, of extrema
+    integer, allocatable:: nearby_extr(:, :) ! (2, :)
+    ! coordinates in projection space of extrema
 
     integer, intent(in):: extr_map(:, :)
     ! At a point of extremum SSH: identification number or this
@@ -39,7 +39,8 @@ contains
          .or. list(selection)%out_cont%closed))
     n_select = size(selection)
     allocate(nearby_extr(2, n_select))
-    forall (l = 1:n_select) nearby_extr(:, l) = list(selection(l))%extr%coord
+    forall (l = 1:n_select) nearby_extr(:, l) &
+         = list(selection(l))%extr%coord_proj
 
   end function nearby_extr