diff --git a/Inst_eddies/inst_eddies.f90 b/Inst_eddies/inst_eddies.f90 index ff2acd5845e792d83bf7639c39e7bb43d44a53aa..83d32d0e1295500d6f8c73b7228da0e96f42d54c 100644 --- a/Inst_eddies/inst_eddies.f90 +++ b/Inst_eddies/inst_eddies.f90 @@ -148,22 +148,20 @@ program inst_eddies loop_extr: do l = 1, s%number_extr i = sorted_extr(l) - associate (e => s%list(i)) - ! Define the geographical window around each eddy extremum: + ! Define the geographical window around each eddy extremum: - llc = e%extr%coord_proj - max_radius - urc = e%extr%coord_proj + max_radius + llc = s%list(i)%extr%coord_proj - max_radius + urc = s%list(i)%extr%coord_proj + max_radius - llc(2) = max(llc(2), 1) - urc(2) = min(urc(2), nlat) + llc(2) = max(llc(2), 1) + urc(2) = min(urc(2), nlat) - if (.not. periodic) then - llc(1) = max(llc(1), 1) - urc(1) = min(urc(1), nlon) - end if + if (.not. periodic) then + llc(1) = max(llc(1), 1) + urc(1) = min(urc(1), nlon) + end if - corner_window = corner + (llc - 1) * step - end associate + corner_window = corner + (llc - 1) * step outside_points = nearby_extr(s%extr_map(llc(1):urc(1), llc(2):urc(2)), & s%list, i)