Skip to content
Snippets Groups Projects
Commit 30283c27 authored by Lionel GUEZ's avatar Lionel GUEZ Committed by GUEZ Lionel
Browse files

Search extrema on the whole domain

This is a waste of computing time but we need to do this because we
are going to remove corner from the dummy arguments of `set_all_extr`.
parent 656a4bc9
No related branches found
No related tags found
No related merge requests found
......@@ -173,9 +173,9 @@ program examine_eddy
end do
n2 = count(in_window)
call set_all_extr(s, step, periodic = .false., &
ssh = ssh(llc(1):urc(1), llc(2):urc(2)), corner = corner_window)
selection = pack(s%extr_map, s%extr_map /= 0)
call set_all_extr(s, step, periodic, ssh, corner)
selection = pack(s%extr_map(llc(1):urc(1), llc(2):urc(2)), &
s%extr_map(llc(1):urc(1), llc(2):urc(2)) /= 0)
selection = pack(selection, s%list(selection)%cyclone .neqv. cyclone)
n1 = size(selection)
allocate(outside_points(2, n1 + n2))
......@@ -189,8 +189,7 @@ program examine_eddy
end if
end do
i = s%extr_map(e%extr%coord_proj(1) - llc(1) + 1, &
e%extr%coord_proj(2) - llc(2) + 1)
i = s%extr_map(e%extr%coord_proj(1), e%extr%coord_proj(2))
e%innermost_level = s%list(i)%innermost_level
call create_cont_list
call set_contours(e%out_cont, e%speed_cont, e%max_speed, cyclone, e%extr, &
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment