diff --git a/Tests/test_max_speed_contour_ssh.f b/Tests/test_max_speed_contour_ssh.f
index 10a18029ea8e1f102d8770cf2b5f3305971a9b79..90f91304a1589eb5c5b6732194ab465a659be772 100644
--- a/Tests/test_max_speed_contour_ssh.f
+++ b/Tests/test_max_speed_contour_ssh.f
@@ -20,8 +20,11 @@ program test_max_speed_contour_ssh
   integer ncid, varid
   real, allocatable:: ssh(:, :) ! (n_lon, n_lat) sea-surface height, in m
   real, allocatable:: u(:, :), v(:, :) ! (n_lon, n_lat) wind, in m s-1
-  integer:: ind_extr(2) = [19, 11] ! indices in the two dimensions of
-                                   ! the extremum
+
+  integer:: ind_extr(2) = [19, 11]
+  ! indices in the two dimensions of the extremum, relative to lower
+  ! left corner
+  
   integer:: radius = 3
 
   namelist /main_nml/ ilon_llc, ilat_llc, ilon_urc, ilat_urc, ind_extr, radius
diff --git a/get_snapshot.f b/get_snapshot.f
index 48515b1a5e19bae2b4e6ec1f44771c67c4ce10b6..d7a4ab3f85c343973a849af29af0c4a7a999c860 100644
--- a/get_snapshot.f
+++ b/get_snapshot.f
@@ -89,8 +89,9 @@ contains
             new_fill_value = ieee_value(0., IEEE_QUIET_NAN))
        call get_var(periodic, max_radius(1), v, ncid, nlon, k, name = "v", &
             new_fill_value = ieee_value(0., IEEE_QUIET_NAN))
-       ! (We will need quiet NaNs rather the original fill values for u
-       ! and v when we compute the max-speed contours.)
+       ! (We will need quiet NaNs rather the original fill values for
+       ! u and v when we compute the max-speed contours and when we
+       ! search the ssh of max-speed contours.)
        call nf95_close(ncid)
 
        call set_all_outerm(s, min_amp, max_radius, step, periodic, ssh, &
diff --git a/get_var.f b/get_var.f
index 57452beb5611a399dd7b6f06e1bd38ef7c9c7b96..9b578ddd77ca86d34e824a11be29371890f855cb 100644
--- a/get_var.f
+++ b/get_var.f
@@ -8,7 +8,7 @@ contains
        new_fill_value)
 
     ! Read a NetCDF variable, change the missing value and extend it
-    ! in longitude.
+    ! in longitude if periodic.
 
     ! Libraries:
     use netcdf95, only: nf95_inq_varid, nf95_get_var, nf95_get_att
@@ -16,7 +16,7 @@ contains
     logical, intent(in):: periodic ! grid is periodic in longitude
     
     integer, intent(in):: max_rad_lon ! maximum radius of an eddy in
-    ! longitude, in number of grid points
+    ! longitude, in number of grid points, used only if periodic
 
     real, intent(out):: values(1 - merge(max_rad_lon, 0, periodic):, :)
     ! (1 - merge(max_rad_lon, 0, periodic):nlon + merge(max_rad_lon, 0,