import xarray as xr
import numpy as np

# Abrir arquivo
ds = xr.open_dataset("test2.nc")
# Selecionar a variável
var = ds["mdepth"]
lon = ds["mplon"]
lat = ds["mplat"]

# Extrair somente os pontos que atendem à condição
valores = var.where(var<=0, drop=True)

print(valores)


idx = np.argwhere(var.values <= 0)
# imprime nesta ordem yc,xc
print(idx)

# nesta ordem yc,xc

print(lon[458,206])
print(lat[458,206])
print(var[458,206])

print(lon[553,169])
print(lat[553,169])
print(var[553,169])