Create procedure set_max_speed (with a large part coming from removed

procedure get_eddy), with corresponding test.

Take into account the fact that derived type polygon became polyline
in library Contour_531.

Use new procedure null_polyline from library Contour_531.

Correct a sneaky bug in mean_speed: the last point in the polygon
array is a duplicate of the first point, do not take the last point
into account when averaging (and interpolating). Discovered this bug
because the mean speed was not exactly the same for equivalent
polygons with circular permutation.

Change philosophy for the tests: read whole NetCDF variables, prepare
different NetCDF files for different geographical regions if needed;
read all data from file or namelist; read names of files on command
line.

Use generic test_program.sh (which is out of this project) instead of
the specific tests_detection_eddies.sh.