Select points and cross sections

cross_section_line(data, start, end)	Obtain an interpolated cross-sectional slice through gridded data.
cross_section_linestring(data, linestring)	Obtain an interpolated cross-sectional slice through gridded data.
points_in_bounds(da, **points)	Returns whether points specified by keyword arguments fall within the bounds of da.
points_values(da, **points)	Get values from specified points.
points_set_values(da, values, **points)	Set values at specified points.
points_indices(da, **points)	Get the indices for points as defined by the arrays x and y.
upper_active_layer(da[, is_ibound, ...])	Function to get the upper active layer from ibound xarray.DataArray