We want to be able to run `inst_eddies` from a directory other than
the directory containing the SHPC to be appended to. Preparing to
write `grid_nml.txt` and `used_config_nml.txt` in the current
directory. To be able to do this in parallel for several slices, we
need to be running in different directories.

Makes more sense.

So the user can get the usage message more easily.

In program `examine_eddy`.

Change order of command line arguments of program
`examine_eddy`. Preparing to allow a variable number of NetCDF file
arguments.

Do not always specify uv-file when it is the same as ssh-file.

Instead of reading command line arguments in procedure `input_ssh`,
read them outside and pass file names to `input_ssh`. This allows
`input_ssh` to be called by programs with different number of command
line arguments, in different positions. Preparing to make the
specification of uv-file optional.

Lionel GUEZ authored 6 months ago and GUEZ Lionel committed 6 months ago

Reverting part of commit c4410682. We want to extend the program to a
non-uniform or non longitude-latitude grid. Checking equality of whole
longitude and latitude arrays, one or two-dimensional, seems a little
bit too much.

Useless since at least commit bff41cdd.

This was not useful since we do not need an array of snapshots in
`inst_eddies`. It is clearer without the integration in a snapshot
variable. The integration in a snapshot variable came from initial
commit, when there was only one program for detection of instantaneous
eddies and overlapping.

Replace argument `s` of `write_snapshot` by its component `list`.

Replace argument `s` of `set_all_extr` by its components. The idea is
to remove the structure `snapshot` from program `inst_edddies`.

We are circumventing a bug of the Intel Fortran compiler. This was
tested with both ifx `2023.1.0` and ifort `2021.9.0`. The bug occurs
on return from procedure `set_max_speed`: `free(): double free
detected in tcache 2` and crashes execution. The bug is revealed in
commit cc13fd5b. The GNU compiler does not have this bug.

To circumvent the Intel compiler bug, instead of an ssh component of
`cont_list_proj`, we have to use a separate array:
`cont_list_ssh`. This is maybe for the best because the component
`area` of `cont_list_proj` was not defined.

`n_points` can be deduced from previous print, and `n_extr` has been
printed too.

If `l != 0` then `shape_rec.record.r_eq_area >= 0` so `n_points` was
always equal to `n_points_valid`: no need for a separate variable.

Because the loop on shapes may take long.

`MAX_RADIUS_deg` is not general enough for a non longitude latitude
grid. The variable `max_radius` moves from module `input_ssh_m` to
module `config_m`. The variable `MAX_RADIUS_deg` disappears. The
constraint `2 * max_radius_deg(1) < 180` becomes `4 * max_radius(1) <
nlon` but is only applied to the periodic case because we do not have
the general equivalent of 180° in number of points.

In procedure `get_var`.

Promote nlon and nlat to variables of module `input_ssh_m`, avoiding
their repeated computation.

Preparing to allow access to variables of module `input_ssh_m`.

Rename variable `outside_points_proj` to `outside_points`.