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[, out_of_bounds])	Get values from specified points.
points_set_values(da, values[, out_of_bounds])	Set values at specified points.
points_indices(da[, out_of_bounds])	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
grid_boundary_xy(grid)	Return grid boundary on the xy plane.
active_grid_boundary_xy(active)	Return active boundary cells on the xy plane.