Interpolate speed in projection coordinates

Inst_eddies/mean_speed.f90

@@ -20,6 +20,8 @@ contains
     use numer_rec_95, only: bilinear_interp2_reg
     use jumble, only: rad_to_deg
 
+    use input_ssh_m, only: corner_whole => corner
+
     real, intent(in):: u(:, :), v(:, :) ! velocity
 
     type(polyline), intent(in):: p
@@ -62,7 +64,8 @@ contains
     end if
 
     ni = p%n_points - 1
-    call bilinear_interp2_reg(corner, step, u, v, p%points(:, :ni), ui, vi)
+    call bilinear_interp2_reg((corner - corner_whole) / step + 1., [1., 1.], &
+         u, v, points_xy(:, :ni), ui, vi)
 
     do j = 1, ni
        x = cos(center(2)) * (p%points(1, j) - center(1))