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

In main program unit `test_successive_overlap`, do not print eddies

for which `delta_out` == 1 and eddies for which `delta_in` == 1. This
is redundant with information in file `edgelist.csv`. I have checked
in test `Successive_overlap_global` that eddies for which `delta_out`
== 1 correspond to column i1 (sorted, unique) and eddies for which
`delta_in` == 1 correspond to column i2 (sorted, unique).
parent 8b656514
No related branches found
No related tags found
No related merge requests found
...@@ -104,30 +104,15 @@ program test_successive_overlap ...@@ -104,30 +104,15 @@ program test_successive_overlap
print *, 'Created file "edgelist.csv".' print *, 'Created file "edgelist.csv".'
print *, snapshot_1, ":" print *, snapshot_1, ":"
print *, "Eddies for which delta_out == 1:"
do i = 1, flow(1)%number_vis_extr
if (flow(1)%list_vis(i)%delta_out /= huge(0)) &
write(unit = *, fmt = "(i0, 1x)", advance = "no") i
end do
print *
print *, "Valid isolated eddies:" print *, "Valid isolated eddies:"
do i = 1, flow(1)%number_vis_extr do i = 1, flow(1)%number_vis_extr
if (flow(1)%list_vis(i)%valid .and. flow(1)%list_vis(i)%delta_out & if (flow(1)%list_vis(i)%valid .and. flow(1)%list_vis(i)%delta_out &
== huge(0)) write(unit = *, fmt = "(i0, 1x)", advance = "no") i == huge(0)) write(unit = *, fmt = "(i0, 1x)", advance = "no") i
end do end do
print * print *
print *, Snapshot_2, ":"
print *, "Eddies for which delta_in == 1:"
do i = 1, flow(2)%number_vis_extr
if (flow(2)%list_vis(i)%delta_in /= huge(0)) &
write(unit = *, fmt = "(i0, 1x)", advance = "no") i
end do
print * print *
print *, Snapshot_2, ":"
print *, "Valid isolated eddies:" print *, "Valid isolated eddies:"
do i = 1, flow(2)%number_vis_extr do i = 1, flow(2)%number_vis_extr
......
...@@ -18,7 +18,8 @@ contains ...@@ -18,7 +18,8 @@ contains
! extremum. 0 at other points. ! extremum. 0 at other points.
type(eddy), intent(in):: list_vis(:) type(eddy), intent(in):: list_vis(:)
! Visible eddies at a given date. We need components ? to be defined. ! Visible eddies at a given date. We need components valid and
! cyclone to be defined.
logical, intent(in):: cyclone ! cyclonicity of the target extremum logical, intent(in):: cyclone ! cyclonicity of the target extremum
......
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