Numerical issues in ERA5 data about atlite HOT 4 CLOSED

import cdsapi
import xarray as xr

product = 'reanalysis-era5-single-levels'
request = {'product_type': 'reanalysis',
           'format': 'netcdf',
           'day': ['01', '02', '03', '04', '05', '06', '07', '08', '09', '10',
                   '11', '12', '13', '14', '15', '16', '17', '18', '19', '20', '21',
                   '22', '23', '24', '25', ' 26', '27', '28', '29', '30', '31'],
           'time': ['00:00', '01:00', '02:00', '03:00', '04:00', '05:00',
                    '06:00', '07:00', '08:00', '09:00', '10:00', '11:00', '12:00',
                    '13:00', '14:00', '15:00', '16:00', '17:00', '18:00', '19:00',
                    '20:00', '21:00', '22:00', '23:00'],
           'area': [71.0, -10.5, 34.5, 34.5],
           'year': 2016,
           'month': 1,
           'variable': ['100m_u_component_of_wind',
                        '100m_v_component_of_wind',
                        '2m_temperature',
                        'runoff',
                        'soil_temperature_level_4',
                        'surface_net_solar_radiation',
                        'surface_pressure',
                        'surface_solar_radiation_downwards',
                        'toa_incident_solar_radiation',
                        'total_sky_direct_solar_radiation_at_surface']}

result = cdsapi.Client().retrieve(
    product,
    request
)

result.download("test.nc")

ds = xr.open_dataset("test.nc")
ds['ro'].to_series().describe()

## Returns
#
# count    1.979561e+07
# mean     1.617490e-05
# std      6.164426e-05
# min     -4.656613e-10
# 25%     -4.656613e-10
# 50%      3.972091e-07
# 75%      1.391862e-05
# max      1.303101e-02
# Name: ro, dtype: float64