Use function `shp_read_point`

Overlap/read_eddy.f90

@@ -13,9 +13,8 @@ contains
     use contour_531, only: null_polyline
     use gpc_f, only: shp_read_pol, polygon
     use jumble, only: deg_to_rad, pi
-    use shapelib, only: shpobject, shpdestroyobject, shpfileisnull
-    use shapelib_03, only: shp_read_object_03, dbf_read_attribute_03, &
-         shp_read_point
+    use shapelib, only: shpfileisnull
+    use shapelib_03, only: dbf_read_attribute_03, shp_read_point
 
     use derived_types, only: eddy, shpc_slice_handler, missing_speed
 
@@ -31,9 +30,8 @@ contains
     integer, intent(in):: ishape ! 0-based
 
     ! Local:
-    TYPE(shpobject) psobject
     type(polygon) p
-    real attr
+    real attr, point(2)
 
     !---------------------------------------------------------------------
 
@@ -53,9 +51,8 @@ contains
 
     ! extr_proj:
     if (.not. shpfileisnull(hshp%extr_proj)) then
-       call shp_read_object_03(hshp%extr_proj, ishape, psobject)
-       e%extr%coord_proj = nint([psobject%padfx(1), psobject%padfy(1)])
-       call shpdestroyobject(psobject)
+       call shp_read_point(point, hshp%extr_proj, ishape)
+       e%extr%coord_proj = nint(point)
     end if
 
     ! outermost_contour: