exemplo4.py

import xarray as xr
import matplotlib.pyplot as plt

""" Plotando a média mensal da ET0 para algumas cidades.
"""

# pegando variavel
path_var = '/home/alexandre/Dropbox/ParaUbuntu/netcdfgrid3/'
ds = xr.open_mfdataset(path_var + 'ETo_daily_UT_Brazil_v2_*.nc', combine='by_coords')
var = ds['ETo']

# cidades e coordenadas
cityNames = ['Santa Maria-RS', 'Manaus-AM',
             'Petrolina-PE', 'Alegre-ES']
cityCoord = [[-29.7, -53.7],
             [-3., -60.],
             [-9.4, -40.5],
             [-20.7, -41.5]]

# calculando a media mensal
varMean = var.resample(time='M').mean('time')

# plotando
for n, city in enumerate(cityNames):
    varMean.sel(latitude=cityCoord[n][0], longitude=cityCoord[n][1],
                method='nearest').plot(label=city)

plt.ylim(1, 8)
plt.title('')
plt.legend(ncol=2)
plt.show()