hard-coding them.

Remove component twice of type eddy, no longer needed.

the logic with flat_extr and noise_around is no longer needed.

set_all_outerm: the value of component suff_amp is modified inside the
first batch of get_1_outerm.

>= 0 everywhere outside of procedure local_extrema.

In test Nearby_extr, use output from Get_snapshot_region_1_noise since
we now use suff_amp in extremum shapefile.

than extr_map /= i.

Preparing to remove the coding of sufficient amplitude in extr_map: do
not use anywhere the sign of extr_map, use the component suff_amp of
the type eddy instead.

In program test_successive_overlap, I would rather not use default
values for namelist variables.

Add a larger test for successive_overlap, to have eddies wih
insufficient amplitude in the region.

Stdin to Input. Two input dir provided for in tests.json: $input_dir
and $large_input_dir.

Rename everywhere variable nearby_extr to outside_points to avoid
confusion with function nearby_extr. Rename file "nearby_extr.csv" to
"outside_points.csv".

In procedure nearby_extr, include in the selection extrema with
different sign, even if they do not have sufficient amplitude.

Bug fix in read_eddy: use the value of read attribute to define
e%interpolated.

Bug fix in test_set_max_speed: forgot conversion to rad.

set_all_outerm and in get_snapshot.

define a selection of extrema. Then get their coordinates. This
prepares the way for a more complex selection: taking cylonicity into account.

procedure get_snapshot. The idea is that, in this branch, nearby_extr
will eventually include extrema depending not only on them having
sufficient amplitude but also on their cyclonicity. Both pieces of
information cannot be simply encoded in extr_map. The downside is that
we compute needlessly nearby_extr if e%radius4 == 1 in set_max_speed.

set_max_speed. Abandon the computation with pack_indices, this was
supposed to give the same result than the computation with argwhere
but be faster if not any(extr_map) < 0. I was not sure of that. I
abandon that because I want to eventually remove the coding of
sufficient amplitude in extr_map.

procedure set_all_outerm. The idea is that, in this branch,
nearby_extr will eventually include extrema depending not only on them
having sufficient amplitude but also on their cyclonicity. Both pieces
of information cannot be simply encoded in extr_map.

title lines in "edgelist_1.csv". Add python script reading
"edgelist_1.csv" into a NetworkX graph.

GPC_F and GPC must be used.

read_snapshot, assert that there is no duplicate tuple in
s%ind_extr. Since we are not sure of that, we must use a do construct
and not a forall in order to define s%extr_map.

Bug fix in program test_local_extrema: support missing values. Add
corresponding test.

Output s%ind_extr in program test_read_snapshot.

be useful for testing successive_overlap.

for upcoming procedure read_snapshot.

Bug fix in procedure rossby_number: e%speed_cont could be null.

multi-polygon. This is a useful test because, in successive_overlap,
the intersection of two polygons may give a multi-polygon.

specifying polygon in source code.

init_shapefiles. This allows us to have variables for DBF field
numbers, as module variables. We avoid the burden of passing this
large number of variables as arguments to write_eddy. Safer to have
variables for DBF field numbers than relying on the order of those
fields.

Remove file_list. This was a duplication of information.

Use named constant for missing max speed value. It is used in two procedures.

to 1e4, safer too.

src_test_get_snapshot.

Move field speed from DBF file max_speed_contour to DBF file extremum
because it can be associated to either the contour in
outermost_contour or in max_speed_contour.

instead of defining the component ssh to the ssh of the extremum or to
the ssh of the outermost contour, use a missing value flag
(1e4). Also, when radius4 >= 2 and the max speed contour is the
outermost contour, use missing value flag in the ssh component of max
speed contour instead of ssh on outermost contour. This makes clearer
output.

Add function null_ssh_contour.

In procedure get_1_outerm, simplify the comparison of level_try and
level_good. So we do not need any longer the dummy argument
ssh_extremum.