CDO

Hi Guilherme,

in Python CDO can return an xarray Dataset or DataArray. Here, I recommend to use the Dataset because eca_cwd returns two variables, consecutive_wet_days_index_per_time_period and number_of_cwd_periods_with_more_than_5days_per_time_period.

Example: Precipitation data has units "kg m-2 s-1" so we have to multiply it with 86400 to get "mm".

import xarray as xr
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
from cdo import Cdo
cdo = Cdo()

infile = 'your-input-data-file.nc'

ds_cwd = cdo.eca_cwd('1.5', input='-mulc,86400 ' + infile, returnXDataset=True)

fig, ax = plt.subplots(subplot_kw=dict(projection=ccrs.PlateCarree()))
ax.coastlines()
ax.pcolormesh(ds_cwd.lon, ds_cwd.lat, cwd.isel(time=0));