import xarray as xr
import matplotlib.pyplot as plt

# abre o netcdf
ds = xr.open_dataset("remap_t_woa23_01_01.nc")

print(ds)  # veja as variáveis e dimensões disponíveis

# supondo que sua variável seja 't'
temp = ds['t_an'].isel(time=0, depth=0)  # seleciona 1 tempo e 1 nível vertical

# plota
plt.figure(figsize=(8,4))
temp.plot(cmap="coolwarm")
#plt.title("Temperatura superficial WOA23")
plt.show()
