Simplifiy computation of `corner_proj`

Inst_eddies/Tests/test_set_max_speed.f90

@@ -36,13 +36,9 @@ program test_set_max_speed
   character(len = :), allocatable:: path
 
   ! Window around the extremum:
-
   integer llc(2) ! indices in global grid of lower left corner
   integer urc(2) ! indices in global grid of upper right corner
 
-  real corner_window(2) ! longitude and latitude of the window around
-  ! the extremum, in rad
-
   integer, parameter:: n_max_cont = 31
 
   type(polyline) cont_list(n_max_cont)
@@ -132,14 +128,11 @@ program test_set_max_speed
      urc(2) = min(urc(2), nlat)
      if (.not. periodic) urc(1) = min(urc(1), nlon)
 
-     corner_window = corner + (llc - 1) * step
-
      call create_cont_list
      call set_max_speed(e%speed_cont, e%max_speed, cont_list, cont_list_proj, &
           n_cont, e%extr%coord, real(e%extr%coord_proj), &
           ssh(llc(1):urc(1), llc(2):urc(2)), u(llc(1):urc(1), llc(2):urc(2)), &
-          v(llc(1):urc(1), llc(2):urc(2)), &
-          (corner_window - corner) / step + 1., step)
+          v(llc(1):urc(1), llc(2):urc(2)), real(llc), step)
      call close_cont_list
      call shpc_create(hshpc, shpc_dir = "SHPC", cyclone = cyclone, slice = 0, &
           grid_lon_lat = .true.)