Test `uniform_lon_lat`

Rather than the value of `s%list(1)%extr%coord_proj(1)`.

Overlap/read_snapshot.f90

@@ -12,7 +12,7 @@ contains
 
     use derived_types, only: snapshot, eddy, shpc_slice_handler, shpc_slice_meta
     use read_eddy_m, only: read_eddy
-    use read_grid_m, only: nlon, nlat, corner, step
+    use read_grid_m, only: nlon, nlat, corner, step, uniform_lon_lat
 
     type(snapshot), intent(out):: s
     ! completely defined except s%list%innermost_level
@@ -65,11 +65,9 @@ contains
 
     ! If coord_proj for all eddies were not read from a shapefile then
     ! we define them here:
-    if (s%number_extr /= 0) then
-       if (s%list(1)%extr%coord_proj(1) == huge(0)) &
-            forall (eddy_i = 1:s%number_extr) s%list(eddy_i)%extr%coord_proj &
-            = nint((s%list(eddy_i)%extr%coord - corner) / step + 1.)
-    end if
+    if (s%number_extr /= 0 .and. uniform_lon_lat) &
+         forall (eddy_i = 1:s%number_extr) s%list(eddy_i)%extr%coord_proj &
+         = nint((s%list(eddy_i)%extr%coord - corner) / step + 1.)
 
     ! Define s%extr_map from s%list%extr%coord_proj: