Polish.

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

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, &

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,