Assume that the outermost contour is non-null

This must be the case, since commit f7c5ec74.

Inst_eddies/write_eddy.f90

@@ -6,6 +6,9 @@ contains
 
   subroutine write_eddy(e, hshpc, date, i)
 
+    ! We assume and do not check that e%out_cont is not null (and
+    ! therefore e%out_cont%area >= 0).
+
     ! Libraries:
     use jumble, only: pi, rad_to_deg, assert
     use shapelib, only: shpt_polygon, shpfileisnull
@@ -48,12 +51,8 @@ contains
             hshpc%extr_proj_eddy_index, i)
     end if
 
-    if (e%out_cont%n_points == 0) then
-       call shp_append_null_03(ishape, hshpc%outermost)
-    else
-       call shp_append_object_03(ishape, hshpc%outermost, shpt_polygon, &
-            e%out_cont%points * rad_to_deg)
-    end if
+    call shp_append_object_03(ishape, hshpc%outermost, shpt_polygon, &
+         e%out_cont%points * rad_to_deg)
 
     if (e%speed_cont%n_points == 0) then
        call shp_append_null_03(ishape, hshpc%max_speed)
@@ -63,14 +62,8 @@ contains
     end if
 
     ! outermost_contour.dbf:
-    if (e%out_cont%area >= 0) then
-       call dbf_write_attribute_03(hshpc%outermost, ishape, &
-            hshpc%out_r_eq_area, sqrt(e%out_cont%area / 1e6 / pi))
-    else
-       call dbf_write_attribute_03(hshpc%outermost, ishape, &
-            hshpc%out_r_eq_area, - 100.)
-    end if
-
+    call dbf_write_attribute_03(hshpc%outermost, ishape, &
+         hshpc%out_r_eq_area, sqrt(e%out_cont%area / 1e6 / pi))
     call dbf_write_attribute_03(hshpc%outermost, ishape, hshpc%out_ssh, &
          e%out_cont%ssh)
     call dbf_write_attribute_03(hshpc%outermost, ishape, hshpc%out_date, date)