def basin_means(data, lat1=2.5, lat2=-2.5):
    """
    Computes the mean in different basins of the equatorial Pacific

    :param data: The data for which the Basins means shall be computed with
    dimension (time, lat, lon).

    :param lat1, lat2: The latidual bounds

    :returns: The mean in the west Pacific (120E- 160E), the central Pacifc (160E-180E),
    east Pacifc  (180E- 240W).
    """
    data_WP = data.loc[dict(lat=slice(lat1, lat2), lon=slice(120, 160))]
    data_WP_mean = data_WP.mean(dim='lat', skipna=True).mean(dim='lon', skipna=True)

    data_CP = data.loc[dict(lat=slice(lat1, lat2), lon=slice(160, 180))]
    data_CP_mean = data_CP.mean(dim='lat', skipna=True).mean(dim='lon', skipna=True)

    data_EP = data.loc[dict(lat=slice(lat1, lat2), lon=slice(180, 240))]
    data_EP_mean = data_EP.mean(dim='lat', skipna=True).mean(dim='lon', skipna=True)

    return data_WP_mean, data_CP_mean, data_EP_mean