3D interpolation Qhull error

[browniea]

Using Interp.griddata ran into a Qhull error when interpolating. This code seemed to fix the issue:

        bad = ~np.isfinite(sigma) | (np.abs(sigma) > 1e6)
        if bad.any():
            sigma[bad] = -1.0
            _ds["mesh2d_interface_sigma"] = xr.DataArray(
                sigma,
                dims=_ds["mesh2d_interface_sigma"].dims,
                attrs=_ds["mesh2d_interface_sigma"].attrs,
            )

[browniea]

I'm investigating this further, I'll plan to add a fix to PR #432

[browniea]

"mesh2d_interface_sigma" should be between 0 (at the water surface) and -1 (at the sea bed). Our data threw an error at the seabed.

[browniea]

I'm assuming the errors would occur at the boundaries 0 and -1, in this case the outliers should be replaced with -1. This is the first example I have of this issue, so I could be wrong with the future errors that this could throw. I've also never had an issue with the with the interface sigma vales in the old naming convention. It works for this case though and can be improved upon the future