Polish

Inst_eddies/get_1_outerm.f90

@@ -55,8 +55,9 @@ contains
     ! points given by longitude, latitude, in rad
 
     type(ssh_contour), intent(out):: cont_list(:) ! (n_max_cont)
-    ! Contour list. Defined only for subscripts 1:n_cont. If n_cont >=
-    ! 1, the outermost contour is element with subscript n_cont. The
+    ! Contour list, with points given by longitude and latitude, in
+    ! radians. Defined only for subscripts 1:n_cont. If n_cont >= 1,
+    ! the outermost contour is element with subscript n_cont. The
     ! contours are in monotonic order of ssh. For i <= n_cont,
     ! cont_list(i)%area has its default initialization value and the
     ! order of points in cont_list(i)%polyline (clockwise or

Inst_eddies/good_contour.f90

@@ -16,6 +16,13 @@ contains
     ! counter-clockwise) is not specified. The contour is in
     ! projection (x, y) coordinates.
 
+    ! Formally, the return value good_contour obeys:
+
+    ! good_contour%closed and good_contour%n_points >= 4 or
+    ! good_contour == null_polyline
+
+    ! good_contour%closed eqv good_contour%n_points /= 0
+
     ! Libraries:
     use contour_531, only: polyline, find_contours_reg_grid, null_polyline
     use geometry, only: polygon_contains_point