Fstd2nc

* Detect newer surface aggregate fields from GEM 5.1 / 5.2 (assign
"sfctype" categories on those fields).
* Detect SWE (from CalDAS) as a surface aggregate field as well.
* New command-line option --sfc-agg-vars to allow the user to specify
additional fields that should be treated as surface-type fields.
* Fix pickling of Buffer objects.
* Fix compatibility issue when using .to_xarray() on Python 3.10.
* Add support for pathlib.Path objects.
* When invoking this module from a Python script, re-route warning and
error messages through the warnings.warn and Exception interfaces
respectively. This will make it easier for the calling script to
filter the warnings and implement exception handling.
* Traduire en (mauvais) français.

* Allow fstd2nc to be "pip install"ed from an environment that doesn't
yet have numpy installed. Installation no longer fails with a weird
error message.
* Fix buggy behaviour when there's a mask field with values > 1.
E.g., values of "2" sometimes represent partially submerged points.
Now treat such points the same as mask value "1" (keep grid point).
* Handle edge case where one of multiple FST input files does not
contain any records.
* Fix --exclude option to properly remove leadtime/reftime again. These
stopped being excludable after some recent code updates in previous
release.
* Fix various bugs in --fstd-compat option.

* New .to_xarray_list() method to return list of variables wrapped in
individual Datasets.
* New command-line --bounds flag, to output grid cell boundaries.
* Remove "leadtime", "reftime", "surface_type", and "station" from
list of coordinates, since some netCDF viewers have trouble with them
being listed in there.
* Fix crash in xarray interface when interpolation is enabled.

* Allow fstd2nc.Buffer objects to be serialized / deserialized.
This allows the objects to be used with the Python multiprocessing
library. Note: not yet compatible with dask multiprocessing, which
seems to use a different mechanism for serialization?
* New keyword args add_px and add_gz for computing pressure or
geopotential height fields in Buffer.to_xarray(). Not available for
command-line usage.
* Allow asymmetric cropping with the --crop-to-smallest-grid option,
i.e., when the smaller grid is not centered inside the larger grid.
* Fixed an address overflow bug when using the xarray interface on
files larger than 4GB.

* Improved memory usage of dask wrapper. Don't pre-load entire block.
* Some miscellaneous speedups to dask wrapper and netCDF conversion.
* New command-line option --only-diag-level, to select only the near-
surface level data (opposite of --ignore-diag-level).
* Fixed --opdict / --vardict parsing for Python 3.9 (it was using a
deprecated xml function).
* Moved 'interp' module to 'gridhacks', since it contains other
non-interpolation features now.
* Added --yin and --yang options to select one of the two pieces of a
yin-yang grid.
* Added --crop-to-smallest-grid option to handle case where there is
a grid that's slightly nested inside of another grid (such as with
GEM output for dynamics vs. permanent bus?). This option will trim
all the data to the smallest of such grids.
* Experimental --turbo option to speed up conversions in the case where
it needs to extract data from many (thousands) of files.
* Allow progress bar to be turned on, even if there does not appear to
be a TTY available. Useful for interactive work within a Jupyter
notebook.

* Fixed a bug with xarray interface when used with masked fields.
Sometimes it was identifying a directory page as the mask of a field,
with catastrophic consequences.