diff --git a/set_max_speed.f b/set_max_speed.f index 13efc47529f7609d68ef51e4b07a8c29bd59e1c4..8c1c2f1d7ba568108c3ad9e9ca40aae3af6e3729 100644 --- a/set_max_speed.f +++ b/set_max_speed.f @@ -39,7 +39,6 @@ contains ! Local: - logical noise_around integer i ! identifying number of the eddy real, allocatable:: nearby_extr(:, :) ! (2, :) longitude and @@ -49,8 +48,6 @@ contains !--------------------------------------------------------------- - noise_around = any(extr_map < 0) - ! Compute e%radius4 from distance in index space, as a lower bound: e%radius4 = floor(polygon_point_dist_2d(convert_to_ind(e%out_cont%points, & corner, step), real(ind_targ_extr))) + 1 @@ -77,7 +74,8 @@ contains e%radius4) i = extr_map(ind_targ_extr(1), ind_targ_extr(2)) - if (noise_around) then + if (any(extr_map < 0)) then + ! Noise around nearby_extr = convert_to_reg_coord(argwhere(extr_map > 0 & .and. extr_map /= i), corner, step) else