Repository: CS-SI/eodag

Issue:
ERA5 download fetches JSON instead of data
Version: '3.0.1.dev7+gfdd3dce8'

Given this code:

```
from eodag import EODataAccessGateway, setup_logging
setup_logging(3)
dag = EODataAccessGateway()
search_result = dag.search(
    start="2019-01-01T00:00:00Z", end="2024-01-01T00:00:00Z",
    provider="dedl", productType="ERA5_SL"
)
product = search_result[0]
product.download()
```

We get data like this:
```
/tmp/ERA5_SL_20190101_20231231_b8d9e155f618bc52665c33913b95dcdb32e52279
└── ERA5_SL_20190101_20231231_b8d9e155f618bc52665c33913b95dcdb32e52279
```
Which is just a JSON file containing something like this:
```
{"description":"Product is not available yet, please try again using given updated location",
"status":"accepted",
"location":"https://hda.data.destination-earth.eu/stac/collections/EO.ECMWF.DAT.REANALYSIS_ERA5_SINGLE_LEVELS/items/ERA5_SL_20190101_20231231_b8d9e155f618bc52665c33913b95dcdb32e52279/download?_dc_qs=%257B%2522data_format%2522%253A%2B%2522grib%2522%252C%2B%2522day%2522%253A%2B%255B%252207%2522%252C%2B%252228%2522%252C%2B%252215%2522%252C%2B%252208%2522%252C%2B%252225%2522%252C%2B%252231%2522%252C%2B%252216%2522%252[...]2022%2522%252C%2B%25222020%2522%252C%2B%25222023%2522%252C%2B%25222021%2522%252C%2B%25222019%2522%255D%257D&processID=reanalysis-era5-single-levels&type=process&jobID=c639c828-c4bd-47c8-9833-c327a23e80d4&status=accepted&created=2024-10-17T16%3A10%3A15.642477&updated=2024-10-17T16%3A10%3A15.642477&links=%7B%27href%27%3A+%27https%3A%2F%2Fcds.climate.copernicus.eu%2Fapi%2Fretrieve%2Fv1%2Fprocesses%2Freanalysis-era5-single-levels%2Fexecution%27%2C+%27rel%27%3A+%27self%27%7D&links=%7B%27href%27%3A+%27https%3A%2F%2Fcds.climate.copernicus.eu%2Fapi%2Fretrieve%2Fv1%2Fjobs%2Fc639c828-c4bd-47c8-9833-c327a23e80d4%27%2C+%27rel%27%3A+%27monitor%27%2C+%27type%27%3A+%27application%2Fjson%27%2C+%27title%27%3A+%27job+status+info%27%7D&metadata=datasetMetadata&orderId=c639c828-c4bd-47c8-9833-c327a23e80d4&provider=copernicus_climate_data_store"}
```

Discussion:

Solution patch:
diff --git a/docs/_static/product_types_information.csv b/docs/_static/product_types_information.csv
index d28c8b0f1..04f6a9485 100644
--- a/docs/_static/product_types_information.csv
+++ b/docs/_static/product_types_information.csv
@@ -227,9 +227,7 @@ SATELLITE_METHANE,"This dataset provides observations of atmospheric methane (CH
 SATELLITE_SEA_ICE_CONCENTRATION,"This dataset provides daily gridded data of sea ice concentration for both hemispheres derived from satellite passive microwave brightness temperatures. Sea ice is an important component of our climate system and a sensitive indicator of climate change. Its presence or its retreat has a strong impact on air-sea interactions, the Earth's energy budget as well as marine ecosystems. It is listed as an Essential Climate Variable by the Global Climate Observing System. Sea ice concentration is defined as the fraction of the ocean surface in a pixel or grid cell that is covered with sea ice. It is one of the parameters commonly used to characterise the sea-ice cover. Other sea ice parameters include sea ice thickness, sea ice edge, and sea ice type, also available in the Climate Data Store.  The dataset consists of two products produced by the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT) Ocean and Sea Ice Satellite Application Facility (OSI SAF) with research & development from European Space Agency Climate Change Initiative projects (ESA CCI):  The Global Sea Ice Concentration Climate Data Record based on measurements from the following sensors: Scanning Multichannel Microwave Radiometer (SMMR; 1978-1987), Special Sensor Microwave/Imager (SSM/I; 1987-2006), and Special Sensor Microwave Imager/Sounder (SSMIS; 2005 onward). This product spans the period from October 1978 to present and is updated daily by an Interim Climate Data Record. In the following, it is referred to as the SSMIS product. The Global Sea Ice Concentration Climate Data Record based on measurements from the Advanced Microwave Scanning Radiometer - Earth Observing System (AMSR-E) sensor (2002-2011) and its successor, AMSR2 (2012-2020). This product spans the 2002-2020 period and is not updated. In the following, it is referred to as the AMSR product. Note, that this product was first produced by the European Space Agency Climate Change Initiative Phase 2 project (ESA CCI) and has been transferred to EUMETSAT OSI SAF since version 3.0. Both products are provided on the same polar projection with a grid resolution of 25 km. However, the AMSR product has a true spatial resolution (as resolved by the sensor) of about 15-25 km versus 30-60 km for the SSMIS product. Therefore, the AMSR product provides a much more detailed view of the sea ice cover than the SSMIS product, especially in the marginal ice zone, the transitional zone between open water and the dense sea ice pack. On the other hand, the clear strength of the SSMIS product is its more than 40-year long and consistent record with daily updates.  The two products share the same algorithm baseline, which is both a continuation of the EUMETSAT OSI SAF approach and a series of innovations contributed by ESA CCI activities. For both products, the underlying algorithm makes use of a combination of the same three temperature channels near 19 GHz and 37 GHz. The data also share a common data format, that allows expert users to revert some of the filtering steps and access the raw output of the SIC algorithms. Both are level-4 products in the sense that gaps are filled by temporal and spatial interpolation. However, gap filling is not applied to fill in days when no input satellite data are available.  Further details about each product can be found below as well as in the Documentation section. ",,,,,"ECMWF,CDS,C3S,sea,ice",ATMOSPHERIC,proprietary,Sea ice concentration daily gridded data from 1978 to present derived from satellite observations,1978-10-25T00:00:00Z,SATELLITE_SEA_ICE_CONCENTRATION,,,,available,,,,,,,,,,,,,,,,,,,,,,,,available,
 SATELLITE_SEA_ICE_EDGE_TYPE,"This dataset provides daily gridded data of sea ice edge and sea ice type derived from brightness temperatures measured by satellite passive microwave radiometers. Sea ice is an important component of our climate system and a sensitive indicator of climate change. Its presence or its retreat has a strong impact on air-sea interactions, the Earth's energy budget as well as marine ecosystems. It is recognized by the Global Climate Observing System as an Essential Climate Variable. Sea ice edge and type are some of the parameters used to characterise sea ice. Other parameters include sea ice concentration and sea ice thickness, also available in the Climate Data Store. Sea ice edge and type are defined as follows: Sea ice edge classifies the sea surface into open water, open ice, and closed ice depending on the amount of sea ice present in each grid cell. This variable is provided for both the Northern and Southern Hemispheres. Note that a sea ice concentration threshold of 30% is used to distinguish between open water and open ice, which differs from the 15% threshold commonly used for other sea ice products such as sea ice extent. Sea ice type classifies ice-covered areas into two categories based on the age of the sea ice: multiyear ice versus seasonal first-year ice. This variable is currently only available for the Northern Hemisphere and limited to the extended boreal winter months (mid-October through April). Sea ice type classification during summer is difficult due to the effect of melting at the ice surface which disturbs the passive microwave signature. Both sea ice products are based on measurements from the series of Scanning Multichannel Microwave Radiometer (SMMR), Special Sensor Microwave/Imager (SSM/I), and Special Sensor Microwave Imager/Sounder (SSMIS) sensors and share the same algorithm baseline. However, sea ice edge makes use of two lower frequencies near 19 GHz and 37 GHz and a higher frequency near 90 GHz whereas sea ice type only uses the two lower frequencies. This dataset combines Climate Data Records (CDRs), which are intended to have sufficient length, consistency, and continuity to assess climate variability and change, and Interim Climate Data Records (ICDRs), which provide regular temporal extensions to the CDRs and where consistency with the CDRs is expected but not extensively checked. For this dataset, both the CDR and ICDR parts of each product were generated using the same software and algorithms. The CDRs of sea ice edge and type currently extend from 25 October 1978 to 31 December 2020 whereas the corresponding ICDRs extend from January 2021 to present (with a 16-day latency behind real time). All data from the current release of the datasets (version 2.0) are Level-4 products, in which data gaps are filled by temporal and spatial interpolation. For product limitations and known issues, please consult the Product User Guide. This dataset is produced on behalf of Copernicus Climate Change Service (C3S), with heritage from the operational products generated by EUMETSAT Ocean and Sea Ice Satellite Application Facility (OSI SAF). Variables in the dataset/application are: Sea ice edge, Sea ice type Variables in the dataset/application are: Status flag, Uncertainty platform: ",,,,,"ECMWF,CDS,C3S,sea,ice",ATMOSPHERIC,proprietary,Sea ice edge and type daily gridded data from 1978 to present derived from satellite observations,1978-10-25T00:00:00Z,SATELLITE_SEA_ICE_EDGE_TYPE,,,,available,,,,,,available,,,,,,,,,,,,,,,,,,available,
 SATELLITE_SEA_ICE_THICKNESS,"This dataset provides monthly gridded data of sea ice thickness for the Arctic region based on satellite radar altimetry observations. Sea ice is an important component of our climate system and a sensitive indicator of climate change. Its presence or its retreat has a strong impact on air-sea interactions, the Earth's energy budget as well as marine ecosystems. It is recognized by the Global Climate Observing System as an Essential Climate Variable. Sea ice thickness is one of the parameters commonly used to characterise sea ice, alongside sea ice concentration, sea ice edge, and sea ice type, also available in the Climate Data Store. Satellite radar altimeters provide measurements of the sea ice freeboard, which is the difference between the height of the surface of sea ice and the surface of water in open leads (areas of open water within the sea ice). Because of the buoyancy of ice in water, typically about 90% of the ice thickness remains under water and thus the total ice thickness is about 10 times the freeboard. However, snow on top of sea ice changes this ratio and complicates the estimation of the ice thickness, requiring the use of auxiliary information about snow depth and density. The retrieval of ice thickness uses the narrow radar swath at the nadir of the satellite at full resolution of approximately 1-10 km and a point spacing of 300 meters. This Level-2 sea-ice thickness products (not provided here) is then gridded for a period of a month to obtain full coverage of a north polar grid at a resolution of 25 km. The algorithm used was developed as part of the European Space Agency Climate Change Initiative (ESA CCI) on Sea Ice. The data provided here are Level-3 Collated (L3C) products: they contain monthly gridded values from orbit data from a single platform (Envisat or CryoSat-2) without interpolation or any other form of gap filling. The files also contain estimates of the algorithm uncertainty as well as a quality status flag indicating potential issues with the retrieval not captured in the algorithm uncertainty. Sources of uncertainty in the algorithm are related to the auxiliary data and to the use of different radar altimeter concepts in Envisat (pulse-limited) and CryoSat-2 (synthetic aperture radar). This dataset combines a Climate Data Record (CDR), which has sufficient length, consistency, and continuity to be used to assess climate variability and change, and an Interim Climate Data Record (ICDR), which provides regular temporal extensions to the CDR and where consistency with the CDR is expected but not extensively checked. Here, the CDR is based on measurements from the RA-2 altimeter on Envisat (October 2002 to October 2010) and the SIRAL altimeter on CryoSat-2 (November 2010 to April 2020). The ICDR is based on observations from CryoSat-2 only (from April 2015 onward) and is updated monthly with a one-month delay behind real time. Users should note that the quality and accuracy of the data record are higher during the CryoSat-2 period than during the Envisat period. As a result, care should be taken when combining the two missions to assess long-term changes and trends. More information can be found in the Product User Guide and Product Quality Assessment Report. This dataset is currently limited spatially to the Arctic region and temporally to the winter months of October through April due to unresolved bias originating from melting snow or open melt ponds in the remaining five months. For a similar reason, no sea-ice thickness data with sufficient quality exist for the Southern Hemisphere. The extension of the CDR/ICDR to other periods, regions, and radar altimeter missions is under development in the extension of the ESA CCI Sea Ice project (ESA CCI+). This dataset is produced on behalf of the Copernicus Climate Change Service (C3S). ",,,,,"ECMWF,CDS,C3S,sea,ice",ATMOSPHERIC,proprietary,Sea ice thickness monthly gridded data for the Arctic from 2002 to present derived from satellite observations,2002-10-01T00:00:00Z,SATELLITE_SEA_ICE_THICKNESS,,,,available,,,,,,,,,,,,,,,,,,,,,,,,available,
-SATELLITE_SEA_LEVEL_BLACK_SEA,"Sea level anomaly is the height of water over the mean sea surface in a given time and region. Up-to-date altimeter standards are used to estimate the sea level anomalies with a mapping algorithm dedicated to the Black sea region. Anomalies are computed with respect to a twenty-year mean reference period (1993-2012). The steady number of reference satellite used in the production of this dataset contributes to the long-term stability of the sea level record. Improvements of the accuracy, sampling of meso-scale processes and of the high-latitude coverage were achieved by using a few additional satellite missions. New data are provided with a delay of about 4-5 months relatively to near-real time or interim sea level products. This delay is mainly due to the timeliness of the input data, the centred processing temporal window and the validation process. However, this processing and validation adds stability and accuracy to the sea level variables and make them adapted to climate applications. This dataset includes uncertainties for each grid cell. More details about the sea level retrieval, additional filters, optimisation procedures, and the error estimation are given in the Documentation section. Variables in the dataset/application are: Absolute dynamic topography, Absolute geostrophic velocity meridian component, Absolute geostrophic velocity zonal component, Geostrophic velocity anomalies meridian component, Geostrophic velocity anomalies zonal component, Sea level anomaly ",,,,,"Climate,ECMWF,CDS,C3S,sea,level,Black Sea",HYDROLOGICAL,proprietary,Sea level daily gridded data from satellite observations for the Black Sea from 1993 to 2020,1993-01-01T00:00:00Z,SATELLITE_SEA_LEVEL_BLACK_SEA,,,,,,,,,,available,,,,,,,,,,,,,,,,,,available,
 SATELLITE_SEA_LEVEL_GLOBAL,"This data set provides gridded daily global estimates of sea level anomaly based on satellite altimetry measurements. The rise in global mean sea level in recent decades has been one of the most important and well-known consequences of climate warming, putting a large fraction of the world population and economic infrastructure at greater risk of flooding. However, changes in the global average sea level mask regional variations that can be one order of magnitude larger. Therefore, it is essential to measure changes in sea level over the world's oceans as accurately as possible. Sea level anomaly is defined as the height of water over the mean sea surface in a given time and region. In this dataset sea level anomalies are computed with respect to a twenty-year mean reference period (1993-2012) using up-to-date altimeter standards. In the past, the altimeter sea level datasets were distributed on the CNES AVISO altimetry portal until their production was taken over by the Copernicus Marine Environment Monitoring Service (CMEMS) and the Copernicus Climate Change Service (C3S) in 2015 and 2016 respectively. The sea level data set provided here by C3S is climate-oriented, that is, dedicated to the monitoring of the long-term evolution of sea level and the analysis of the ocean/climate indicators, both requiring a homogeneous and stable sea level record. To achieve this, a steady two-satellite merged constellation is used at all time steps in the production system: one satellite serves as reference and ensures the long-term stability of the data record; the other satellite (which varies across the record) is used to improve accuracy, sample mesoscale processes and provide coverage at high latitudes. The C3S sea level data set is used to produce Ocean Monitoring Indicators (e.g. global and regional mean sea level evolution), available in the CMEMS catalogue. The CMEMS sea level dataset has a more operational focus as it is dedicated to the retrieval of mesoscale signals in the context of ocean modeling and analysis of the ocean circulation on a global or regional scale. Such applications require the most accurate sea level estimates at each time step with the best spatial sampling of the ocean with all satellites available, with less emphasis on long-term stability and homogeneity. This data set is updated three times a year with a delay of about 6 months relative to present time. This delay is mainly due to the timeliness of the input data, the centred processing temporal window and the validation process. However, these processing and validation steps are essential to enhance the stability and accuracy of the sea level products and make them suitable for climate applications. This dataset includes estimates of sea level anomaly and absolute dynamic topography together with the corresponding geostrophic velocities. More details about the sea level retrieval algorithms, additional filters, optimisation procedures, and the error estimation are given in the Documentation tab. Variables in the dataset/application are: Absolute dynamic topography, Absolute geostrophic velocity meridian component, Absolute geostrophic velocity zonal component, Geostrophic velocity anomalies meridian component, Geostrophic velocity anomalies zonal component, Sea level anomaly ",,,,,"Climate,ECMWF,CDS,C3S,sea,level,global",HYDROLOGICAL,proprietary,Sea level gridded data from satellite observations for the global ocean,1993-01-01T00:00:00Z,SATELLITE_SEA_LEVEL_GLOBAL,,,,available,,,,,,available,,,,,,,,,,,,,,,,,,available,
-SATELLITE_SEA_LEVEL_MEDITERRANEAN,"Sea level anomaly is the height of water over the mean sea surface in a given time and region. In this dataset sea level anomalies are computed with respect to a twenty-year mean reference period (1993-2012). Up-to-date altimeter standards are used to estimate the sea level anomalies with a mapping algorithm specifically dedicated to the Mediterranean Sea. The steady number of reference satellite used in the production of this dataset contributes to the long-term stability of the sea level record. Improvements of the accuracy, sampling of meso-scale processes and of the high-latitude coverage were achieved by using a few additional satellite missions. New data are provided with a delay of about 4-5 months relatively to near-real time or interim sea level products. This delay is mainly due to the timeliness of the input data, the centred processing temporal window and the validation process. However, this processing and validation adds stability and accuracy to the sea level variables and make them adapted to climate applications. This dataset includes uncertainties for each grid cell. More details about the sea level retrieval, additional filters, optimisation procedures, and the error estimation are given in the Documentation section. Variables in the dataset/application are: Absolute dynamic topography, Absolute geostrophic velocity meridian component, Absolute geostrophic velocity zonal component, Geostrophic velocity anomalies meridian component, Geostrophic velocity anomalies zonal component, Sea level anomaly ",,,,,"Climate,ECMWF,CDS,C3S,sea,level,mediterranean",HYDROLOGICAL,proprietary,Sea level daily gridded data from satellite observations for the Mediterranean Sea,1993-01-01T00:00:00Z,SATELLITE_SEA_LEVEL_MEDITERRANEAN,,,,,,,,,,available,,,,,,,,,,,,,,,,,,available,
 SEASONAL_MONTHLY_PL,"This entry covers pressure-level data aggregated on a monthly time resolution. \nSeasonal forecasts provide a long-range outlook of changes in the Earth system over periods of a few weeks or months, as a result of predictable changes in some of the slow-varying components of the system. For example, ocean temperatures typically vary slowly, on timescales of weeks or months; as the ocean has an impact on the overlaying atmosphere, the variability of its properties (e.g. temperature) can modify both local and remote atmospheric conditions. Such modifications of the 'usual' atmospheric conditions are the essence of all long-range (e.g. seasonal) forecasts. This is different from a weather forecast, which gives a lot more precise detail - both in time and space - of the evolution of the state of the atmosphere over a few days into the future. Beyond a few days, the chaotic nature of the atmosphere limits the possibility to predict precise changes at local scales. This is one of the reasons long-range forecasts of atmospheric conditions have large uncertainties. To quantify such uncertainties, long-range forecasts use ensembles, and meaningful forecast products reflect a distributions of outcomes.\nGiven the complex, non-linear interactions between the individual components of the Earth system, the best tools for long-range forecasting are climate models which include as many of the key components of the system and possible; typically, such models include representations of the atmosphere, ocean and land surface. These models are initialised with data describing the state of the system at the starting point of the forecast, and used to predict the evolution of this state in time.\nWhile uncertainties coming from imperfect knowledge of the initial conditions of the components of the Earth system can be described with the use of ensembles, uncertainty arising from approximations made in the models are very much dependent on the choice of model. A convenient way to quantify the effect of these approximations is to combine outputs from several models, independently developed, initialised and operated.\nTo this effect, the C3S provides a multi-system seasonal forecast service, where data produced by state-of-the-art seasonal forecast systems developed, implemented and operated at forecast centres in several European countries is collected, processed and combined to enable user-relevant applications. The composition of the C3S seasonal multi-system and the full content of the database underpinning the service are described in the documentation. The data is grouped in several catalogue entries (CDS datasets), currently defined by the type of variable (single-level or multi-level, on pressure surfaces) and the level of post-processing applied (data at original time resolution, processing on temporal aggregation and post-processing related to bias adjustment).\nThe variables available in this data set are listed in the table below. The data includes forecasts created in real-time (since 2017) and retrospective forecasts (hindcasts) initialised at equivalent intervals during the period 1993-2016.\n\nVariables in the dataset/application are:\nGeopotential, Specific humidity, Temperature, U-component of wind, V-component of wind ",,,,,"ECMWF,CDS,C3S,seasonal,forecast,monthly,pressure,levels",ATMOSPHERIC,proprietary,Seasonal forecast monthly statistics on pressure levels,1981-01-01T00:00:00Z,SEASONAL_MONTHLY_PL,,,,available,,,,,,available,,,,,,,,,,,,,,,,,,available,
 SEASONAL_MONTHLY_SL,"This entry covers single-level data aggregated on a monthly time resolution. \nSeasonal forecasts provide a long-range outlook of changes in the Earth system over periods of a few weeks or months, as a result of predictable changes in some of the slow-varying components of the system. For example, ocean temperatures typically vary slowly, on timescales of weeks or months; as the ocean has an impact on the overlaying atmosphere, the variability of its properties (e.g. temperature) can modify both local and remote atmospheric conditions. Such modifications of the 'usual' atmospheric conditions are the essence of all long-range (e.g. seasonal) forecasts. This is different from a weather forecast, which gives a lot more precise detail - both in time and space - of the evolution of the state of the atmosphere over a few days into the future. Beyond a few days, the chaotic nature of the atmosphere limits the possibility to predict precise changes at local scales. This is one of the reasons long-range forecasts of atmospheric conditions have large uncertainties. To quantify such uncertainties, long-range forecasts use ensembles, and meaningful forecast products reflect a distributions of outcomes.\nGiven the complex, non-linear interactions between the individual components of the Earth system, the best tools for long-range forecasting are climate models which include as many of the key components of the system and possible; typically, such models include representations of the atmosphere, ocean and land surface. These models are initialised with data describing the state of the system at the starting point of the forecast, and used to predict the evolution of this state in time.\nWhile uncertainties coming from imperfect knowledge of the initial conditions of the components of the Earth system can be described with the use of ensembles, uncertainty arising from approximations made in the models are very much dependent on the choice of model. A convenient way to quantify the effect of these approximations is to combine outputs from several models, independently developed, initialised and operated.\nTo this effect, the C3S provides a multi-system seasonal forecast service, where data produced by state-of-the-art seasonal forecast systems developed, implemented and operated at forecast centres in several European countries is collected, processed and combined to enable user-relevant applications. The composition of the C3S seasonal multi-system and the full content of the database underpinning the service are described in the documentation. The data is grouped in several catalogue entries (CDS datasets), currently defined by the type of variable (single-level or multi-level, on pressure surfaces) and the level of post-processing applied (data at original time resolution, processing on temporal aggregation and post-processing related to bias adjustment).\nThe variables available in this data set are listed in the table below. The data includes forecasts created in real-time (since 2017) and retrospective forecasts (hindcasts) initialised at equivalent intervals during the period 1993-2016.\n\nVariables in the dataset/application are:\n10m u-component of wind, 10m v-component of wind, 10m wind gust since previous post-processing, 10m wind speed, 2m dewpoint temperature, 2m temperature, East-west surface stress rate of accumulation, Evaporation, Maximum 2m temperature in the last 24 hours, Mean sea level pressure, Mean sub-surface runoff rate, Mean surface runoff rate, Minimum 2m temperature in the last 24 hours, North-south surface stress rate of accumulation, Runoff, Sea surface temperature, Sea-ice cover, Snow density, Snow depth, Snowfall, Soil temperature level 1, Solar insolation rate of accumulation, Surface latent heat flux, Surface sensible heat flux, Surface solar radiation, Surface solar radiation downwards, Surface thermal radiation, Surface thermal radiation downwards, Top solar radiation, Top thermal radiation, Total cloud cover, Total precipitation ",,,,,"ECMWF,CDS,C3S,seasonal,forecast,monthly,single,levels",ATMOSPHERIC,proprietary,Seasonal forecast monthly statistics on single levels,1981-01-01T00:00:00Z,SEASONAL_MONTHLY_SL,,,,available,,,,,,available,,,,,,,,,,,,,,,,,,available,
 SEASONAL_ORIGINAL_PL,"his entry covers pressure-level data at the original time resolution (once every 12 hours). \nSeasonal forecasts provide a long-range outlook of changes in the Earth system over periods of a few weeks or months, as a result of predictable changes in some of the slow-varying components of the system. For example, ocean temperatures typically vary slowly, on timescales of weeks or months; as the ocean has an impact on the overlaying atmosphere, the variability of its properties (e.g. temperature) can modify both local and remote atmospheric conditions. Such modifications of the 'usual' atmospheric conditions are the essence of all long-range (e.g. seasonal) forecasts. This is different from a weather forecast, which gives a lot more precise detail - both in time and space - of the evolution of the state of the atmosphere over a few days into the future. Beyond a few days, the chaotic nature of the atmosphere limits the possibility to predict precise changes at local scales. This is one of the reasons long-range forecasts of atmospheric conditions have large uncertainties. To quantify such uncertainties, long-range forecasts use ensembles, and meaningful forecast products reflect a distributions of outcomes.\nGiven the complex, non-linear interactions between the individual components of the Earth system, the best tools for long-range forecasting are climate models which include as many of the key components of the system and possible; typically, such models include representations of the atmosphere, ocean and land surface. These models are initialised with data describing the state of the system at the starting point of the forecast, and used to predict the evolution of this state in time.\nWhile uncertainties coming from imperfect knowledge of the initial conditions of the components of the Earth system can be described with the use of ensembles, uncertainty arising from approximations made in the models are very much dependent on the choice of model. A convenient way to quantify the effect of these approximations is to combine outputs from several models, independently developed, initialised and operated.\nTo this effect, the C3S provides a multi-system seasonal forecast service, where data produced by state-of-the-art seasonal forecast systems developed, implemented and operated at forecast centres in several European countries is collected, processed and combined to enable user-relevant applications. The composition of the C3S seasonal multi-system and the full content of the database underpinning the service are described in the documentation. The data is grouped in several catalogue entries (CDS datasets), currently defined by the type of variable (single-level or multi-level, on pressure surfaces) and the level of post-processing applied (data at original time resolution, processing on temporal aggregation and post-processing related to bias adjustment).\nThe variables available in this data set are listed in the table below. The data includes forecasts created in real-time (since 2017) and retrospective forecasts (hindcasts) initialised at equivalent intervals during the period 1993-2016.\n\nVariables in the dataset/application are:\nGeopotential, Specific humidity, Temperature, U-component of wind, V-component of wind ",,,,,"ECMWF,CDS,C3S,seasonal,forecast,subdaily,pressure,levels",ATMOSPHERIC,proprietary,Seasonal forecast subdaily data on pressure levels,1981-01-01T00:00:00Z,SEASONAL_ORIGINAL_PL,,,,available,,,,,,available,,,,,,,,,,,,,,,,,,available,
diff --git a/docs/notebooks/api_user_guide/2_providers_products_available.ipynb b/docs/notebooks/api_user_guide/2_providers_products_available.ipynb
index 696ff46de..4a3eb3eb4 100644
--- a/docs/notebooks/api_user_guide/2_providers_products_available.ipynb
+++ b/docs/notebooks/api_user_guide/2_providers_products_available.ipynb
@@ -485,9 +485,7 @@
        " 'SATELLITE_SEA_ICE_CONCENTRATION',\n",
        " 'SATELLITE_SEA_ICE_EDGE_TYPE',\n",
        " 'SATELLITE_SEA_ICE_THICKNESS',\n",
-       " 'SATELLITE_SEA_LEVEL_BLACK_SEA',\n",
        " 'SATELLITE_SEA_LEVEL_GLOBAL',\n",
-       " 'SATELLITE_SEA_LEVEL_MEDITERRANEAN',\n",
        " 'SEASONAL_MONTHLY_PL',\n",
        " 'SEASONAL_MONTHLY_SL',\n",
        " 'SEASONAL_ORIGINAL_PL',\n",
diff --git a/docs/notebooks/tutos/tuto_wekeo.ipynb b/docs/notebooks/tutos/tuto_wekeo.ipynb
index bc6e8b5cd..995edb807 100644
--- a/docs/notebooks/tutos/tuto_wekeo.ipynb
+++ b/docs/notebooks/tutos/tuto_wekeo.ipynb
@@ -148,9 +148,7 @@
        " 'SATELLITE_SEA_ICE_CONCENTRATION',\n",
        " 'SATELLITE_SEA_ICE_EDGE_TYPE',\n",
        " 'SATELLITE_SEA_ICE_THICKNESS',\n",
-       " 'SATELLITE_SEA_LEVEL_BLACK_SEA',\n",
        " 'SATELLITE_SEA_LEVEL_GLOBAL',\n",
-       " 'SATELLITE_SEA_LEVEL_MEDITERRANEAN',\n",
        " 'SEASONAL_MONTHLY_PL',\n",
        " 'SEASONAL_MONTHLY_SL',\n",
        " 'SEASONAL_ORIGINAL_PL',\n",
diff --git a/eodag/resources/product_types.yml b/eodag/resources/product_types.yml
index 5424e48ba..686b100d4 100644
--- a/eodag/resources/product_types.yml
+++ b/eodag/resources/product_types.yml
@@ -3286,32 +3286,6 @@ SEASONAL_POSTPROCESSED_PL:
   title: Seasonal forecast anomalies on pressure levels
   missionStartDate: "2017-09-01T00:00:00Z"
 
-SATELLITE_SEA_LEVEL_BLACK_SEA:
-  abstract: |
-    Sea level anomaly is the height of water over the mean sea surface in a given time and region. Up-to-date
-    altimeter standards are used to estimate the sea level anomalies with a mapping algorithm dedicated to the
-    Black sea region. Anomalies are computed with respect to a twenty-year mean reference period (1993-2012).
-    The steady number of reference satellite used in the production of this dataset contributes to the long-term
-    stability of the sea level record. Improvements of the accuracy, sampling of meso-scale processes and of the
-    high-latitude coverage were achieved by using a few additional satellite missions. New data are provided with
-    a delay of about 4-5 months relatively to near-real time or interim sea level products. This delay is mainly
-    due to the timeliness of the input data, the centred processing temporal window and the validation process.
-    However, this processing and validation adds stability and accuracy to the sea level variables and make them
-    adapted to climate applications. This dataset includes uncertainties for each grid cell. More details about
-    the sea level retrieval, additional filters, optimisation procedures, and the error estimation are given in
-    the Documentation section. Variables in the dataset/application are: Absolute dynamic topography, Absolute
-    geostrophic velocity meridian component, Absolute geostrophic velocity zonal component, Geostrophic velocity
-    anomalies meridian component, Geostrophic velocity anomalies zonal component, Sea level anomaly
-  platform:
-  instrument:
-  platformSerialIdentifier:
-  processingLevel:
-  keywords: Climate,ECMWF,CDS,C3S,sea,level,Black Sea
-  sensorType: HYDROLOGICAL
-  license: proprietary
-  title: Sea level daily gridded data from satellite observations for the Black Sea from 1993 to 2020
-  missionStartDate: "1993-01-01T00:00:00Z"
-
 SATELLITE_SEA_LEVEL_GLOBAL:
   abstract: |
     This data set provides gridded daily global estimates of sea level anomaly based on satellite altimetry
@@ -3355,33 +3329,6 @@ SATELLITE_SEA_LEVEL_GLOBAL:
   missionStartDate: "1993-01-01T00:00:00Z"
   missionEndDate: "2022-08-04T23:59:59Z"
 
-SATELLITE_SEA_LEVEL_MEDITERRANEAN:
-  abstract: |
-    Sea level anomaly is the height of water over the mean sea surface in a given time and region. In this dataset sea
-    level anomalies are computed with respect to a twenty-year mean reference period (1993-2012). Up-to-date altimeter
-    standards are used to estimate the sea level anomalies with a mapping algorithm specifically dedicated to the
-    Mediterranean Sea. The steady number of reference satellite used in the production of this dataset contributes to
-    the long-term stability of the sea level record. Improvements of the accuracy, sampling of meso-scale processes and
-    of the high-latitude coverage were achieved by using a few additional satellite missions. New data are provided with
-    a delay of about 4-5 months relatively to near-real time or interim sea level products. This delay is mainly due to
-    the timeliness of the input data, the centred processing temporal window and the validation process. However, this
-    processing and validation adds stability and accuracy to the sea level variables and make them adapted to climate
-    applications. This dataset includes uncertainties for each grid cell. More details about the sea level retrieval,
-    additional filters, optimisation procedures, and the error estimation are given in the Documentation section.
-    Variables in the dataset/application are: Absolute dynamic topography, Absolute geostrophic velocity meridian
-    component, Absolute geostrophic velocity zonal component, Geostrophic velocity anomalies meridian component,
-    Geostrophic velocity anomalies zonal component, Sea level anomaly
-  platform:
-  instrument:
-  platformSerialIdentifier:
-  processingLevel:
-  keywords: Climate,ECMWF,CDS,C3S,sea,level,mediterranean
-  sensorType: HYDROLOGICAL
-  license: proprietary
-  title: Sea level daily gridded data from satellite observations for the Mediterranean Sea
-  missionStartDate: "1993-01-01T00:00:00Z"
-  missionEndDate: "2018-11-01T23:59:59Z"
-
 SEASONAL_POSTPROCESSED_SL:
   abstract: |
     This entry covers single-level data post-processed for bias adjustment on a monthly time resolution.
diff --git a/eodag/resources/providers.yml b/eodag/resources/providers.yml
index 76ba3ce4b..a117f4ce2 100644
--- a/eodag/resources/providers.yml
+++ b/eodag/resources/providers.yml
@@ -2624,9 +2624,6 @@
         output_extension: .grib
       UERRA_EUROPE_SL:
         output_extension: .grib
-      SATELLITE_SEA_LEVEL_BLACK_SEA:
-        output_extension: .zip
-        extract: true
       GLACIERS_DIST_RANDOLPH:
         output_extension: .zip
         extract: true
@@ -4772,18 +4769,6 @@
         <<: *day_month_year
         defaultGeometry: 'POLYGON((180 -90, 180 90, -180 90, -180 -90, 180 -90))'
         orderLink: 'https://gateway.prod.wekeo2.eu/hda-broker/api/v1/dataaccess/download?{{"location": "{downloadLink}","product_id":"{id}", "dataset_id": "EO:ECMWF:DAT:SATELLITE_SEA_ICE_CONCENTRATION"}}'
-    SATELLITE_SEA_LEVEL_BLACK_SEA:
-      productType: EO:ECMWF:DAT:SATELLITE_SEA_LEVEL_BLACK_SEA
-      variable: all
-      format: zip
-      metadata_mapping:
-        id: '$.id'
-        <<: *day_month_year
-        variable:
-          - '{{"variable": "{variable}"}}'
-          - '$.null'
-        defaultGeometry: 'POLYGON((180 -90, 180 90, -180 90, -180 -90, 180 -90))'
-        orderLink: 'https://gateway.prod.wekeo2.eu/hda-broker/api/v1/dataaccess/download?{{"location": "{downloadLink}","product_id":"{id}", "dataset_id": "EO:ECMWF:DAT:SATELLITE_SEA_LEVEL_BLACK_SEA"}}'
     SATELLITE_SEA_LEVEL_GLOBAL:
       productType: EO:ECMWF:DAT:SATELLITE_SEA_LEVEL_GLOBAL
       variable:
@@ -4799,18 +4784,6 @@
         <<: *variable_list
         defaultGeometry: 'POLYGON((180 -90, 180 90, -180 90, -180 -90, 180 -90))'
         orderLink: 'https://gateway.prod.wekeo2.eu/hda-broker/api/v1/dataaccess/download?{{"location": "{downloadLink}","product_id":"{id}", "dataset_id": "EO:ECMWF:DAT:SATELLITE_SEA_LEVEL_GLOBAL"}}'
-    SATELLITE_SEA_LEVEL_MEDITERRANEAN:
-      productType: EO:ECMWF:DAT:SATELLITE_SEA_LEVEL_MEDITERRANEAN
-      variable: all
-      format: zip
-      metadata_mapping:
-        id: '$.id'
-        <<: *day_month_year
-        variable:
-          - '{{"variable": "{variable}"}}'
-          - '$.null'
-        defaultGeometry: 'POLYGON((180 -90, 180 90, -180 90, -180 -90, 180 -90))'
-        orderLink: 'https://gateway.prod.wekeo2.eu/hda-broker/api/v1/dataaccess/download?{{"location": "{downloadLink}","product_id":"{id}", "dataset_id": "EO:ECMWF:DAT:SATELLITE_SEA_LEVEL_MEDITERRANEAN"}}'
     SEASONAL_ORIGINAL_SL:
       productType: EO:ECMWF:DAT:SEASONAL_ORIGINAL_SINGLE_LEVELS
       variable:
@@ -6215,48 +6188,84 @@
     # ECMWF - COPERNICUS CDS
     ERA5_SL:
       productType: EO.ECMWF.DAT.REANALYSIS_ERA5_SINGLE_LEVELS
+      metadata_mapping:
+        <<: *orderable_mm
     ERA5_SL_MONTHLY:
       productType: EO.ECMWF.DAT.REANALYSIS_ERA5_SINGLE_LEVELS_MONTHLY_MEANS
+      metadata_mapping:
+        <<: *orderable_mm
     ERA5_PL:
       productType: EO.ECMWF.DAT.ERA5_HOURLY_VARIABLES_ON_PRESSURE_LEVELS
+      metadata_mapping:
+        <<: *orderable_mm
     ERA5_PL_MONTHLY:
       productType: EO.ECMWF.DAT.ERA5_MONTHLY_MEANS_VARIABLES_ON_PRESSURE_LEVELS
+      metadata_mapping:
+        <<: *orderable_mm
     ERA5_LAND:
       productType: EO.ECMWF.DAT.ERA5_LAND_HOURLY
+      metadata_mapping:
+        <<: *orderable_mm
     ERA5_LAND_MONTHLY:
       productType: EO.ECMWF.DAT.ERA5_LAND_MONTHLY
+      metadata_mapping:
+        <<: *orderable_mm
     UERRA_EUROPE_SL:
       productType: EO.ECMWF.DAT.REANALYSIS_UERRA_EUROPE_SINGLE_LEVELS
+      metadata_mapping:
+        <<: *orderable_mm
     GRIDDED_GLACIERS_MASS_CHANGE:
-      productType: EO.ECMWF.DAT.GLACIERS_ELEVATION_AND_MASS_CHANGE_DATA_1850_PRESENT
+      productType: EO.ECMWF.DAT.DERIVED_GRIDDED_GLACIER_MASS_CHANGE
+      metadata_mapping:
+        <<: *orderable_mm
     GLACIERS_DIST_RANDOLPH:
       productType: EO.ECMWF.DAT.GLACIERS_DISTRIBUTION_DATA_FROM_RANDOLPH_GLACIER_INVENTORY_2000
+      metadata_mapping:
+        <<: *orderable_mm
     SATELLITE_CARBON_DIOXIDE:
       productType: EO.ECMWF.DAT.CO2_DATA_FROM_SATELLITE_SENSORS_2002_PRESENT
+      metadata_mapping:
+        <<: *orderable_mm
     SATELLITE_METHANE:
       productType: EO.ECMWF.DAT.METHANE_DATA_SATELLITE_SENSORS_2002_PRESENT
+      metadata_mapping:
+        <<: *orderable_mm
     SEASONAL_POSTPROCESSED_PL:
       productType: EO.ECMWF.DAT.SEASONAL_FORECAST_ANOMALIES_ON_PRESSURE_LEVELS_2017_PRESENT
-    SATELLITE_SEA_LEVEL_BLACK_SEA:
-      productType: EO.ECMWF.DAT.SEA_LEVEL_DAILY_GRIDDED_DATA_FOR_BLACK_SEA_1993_PRESENT
+      metadata_mapping:
+        <<: *orderable_mm
     SATELLITE_SEA_LEVEL_GLOBAL:
       productType: EO.ECMWF.DAT.SEA_LEVEL_DAILY_GRIDDED_DATA_FOR_GLOBAL_OCEAN_1993_PRESENT
-    SATELLITE_SEA_LEVEL_MEDITERRANEAN:
-      productType: EO.ECMWF.DAT.SEA_LEVEL_DAILY_GRIDDED_DATA_FOR_MEDITERRANEAN_SEA_1993_PRESENT
+      metadata_mapping:
+        <<: *orderable_mm
     SATELLITE_SEA_ICE_EDGE_TYPE:
       productType: EO.ECMWF.DAT.SATELLITE_SEA_ICE_EDGE_TYPE
+      metadata_mapping:
+        <<: *orderable_mm
     SEASONAL_POSTPROCESSED_SL:
       productType: EO.ECMWF.DAT.SEASONAL_FORECAST_ANOMALIES_ON_SINGLE_LEVELS_2017_PRESENT
+      metadata_mapping:
+        <<: *orderable_mm
     SEASONAL_ORIGINAL_SL:
       productType: EO.ECMWF.DAT.SEASONAL_FORECAST_DAILY_DATA_ON_SINGLE_LEVELS_2017_PRESENT
+      metadata_mapping:
+        <<: *orderable_mm
     SEASONAL_ORIGINAL_PL:
       productType: EO.ECMWF.DAT.SEASONAL_FORECAST_DAILY_DATA_ON_PRESSURE_LEVELS_2017_PRESENT
+      metadata_mapping:
+        <<: *orderable_mm
     SEASONAL_MONTHLY_PL:
       productType: EO.ECMWF.DAT.SEASONAL_FORECAST_MONTHLY_STATISTICS_ON_PRESSURE_LEVELS_2017_PRESENT
+      metadata_mapping:
+        <<: *orderable_mm
     SEASONAL_MONTHLY_SL:
       productType: EO.ECMWF.DAT.SEASONAL_FORECAST_MONTHLY_STATISTICS_ON_SINGLE_LEVELS_2017_PRESENT
+      metadata_mapping:
+        <<: *orderable_mm
     SIS_HYDRO_MET_PROJ:
       productType: EO.ECMWF.DAT.SIS_HYDROLOGY_METEOROLOGY_DERIVED_PROJECTIONS
+      metadata_mapping:
+        <<: *orderable_mm
     # ECMWF - CEMS
     FIRE_HISTORICAL:
       productType: EO.ECMWF.DAT.CEMS_FIRE_HISTORICAL
@@ -6283,30 +6292,56 @@
     # COPERNICUS ADS
     CAMS_GAC_FORECAST:
       productType: EO.ECMWF.DAT.CAMS_GLOBAL_ATMOSHERIC_COMPO_FORECAST
+      metadata_mapping:
+        <<: *orderable_mm
     CAMS_EU_AIR_QUALITY_FORECAST:
       productType: EO.ECMWF.DAT.CAMS_EUROPE_AIR_QUALITY_FORECASTS
+      metadata_mapping:
+        <<: *orderable_mm
     CAMS_GFE_GFAS:
       productType: EO.ECMWF.DAT.CAMS_GLOBAL_FIRE_EMISSIONS_GFAS
+      metadata_mapping:
+        <<: *orderable_mm
     CAMS_SOLAR_RADIATION:
       productType: EO.ECMWF.DAT.CAMS_SOLAR_RADIATION_TIMESERIES
+      metadata_mapping:
+        <<: *orderable_mm
     CAMS_GREENHOUSE_INVERSION:
       productType: EO.ECMWF.DAT.CAMS_GREENHOUSE_GAS_FLUXES
+      metadata_mapping:
+        <<: *orderable_mm
     CAMS_EAC4_MONTHLY:
       productType: EO.ECMWF.DAT.CAMS_GLOBAL_REANALYSIS_EAC4_MONTHLY_AV_FIELDS
+      metadata_mapping:
+        <<: *orderable_mm
     CAMS_EU_AIR_QUALITY_RE:
       productType: EO.ECMWF.DAT.CAMS_EUROPE_AIR_QUALITY_REANALYSES
+      metadata_mapping:
+        <<: *orderable_mm
     CAMS_EAC4:
       productType: EO.ECMWF.DAT.CAMS_GLOBAL_REANALYSIS_EAC4
+      metadata_mapping:
+        <<: *orderable_mm
     CAMS_GRF_AUX:
       productType: EO.ECMWF.DAT.CAMS_GLOBAL_RADIATIVE_FORCING_AUX
+      metadata_mapping:
+        <<: *orderable_mm
     CAMS_GRF:
       productType: EO.ECMWF.DAT.CAMS_GLOBAL_RADIATIVE_FORCING
+      metadata_mapping:
+        <<: *orderable_mm
     CAMS_GREENHOUSE_EGG4_MONTHLY:
       productType: EO.ECMWF.DAT.CAMS_GLOBAL_GREENHOUSE_GAS_REANALYSIS_MONTHLY_AV_FIELDS
+      metadata_mapping:
+        <<: *orderable_mm
     CAMS_GREENHOUSE_EGG4:
       productType: EO.ECMWF.DAT.CAMS_GLOBAL_GREENHOUSE_GAS_REANALYSIS
+      metadata_mapping:
+        <<: *orderable_mm
     CAMS_GLOBAL_EMISSIONS:
       productType: EO.ECMWF.DAT.CAMS_GLOBAL_EMISSION_INVENTORIES
+      metadata_mapping:
+        <<: *orderable_mm
     # COPERNICUS ADS -  Digital Elevation Model
     COP_DEM_GLO30_DGED:
       productType: EO.DEM.DAT.COP-DEM_GLO-30-DGED


Test patch:
diff --git a/tests/units/test_core.py b/tests/units/test_core.py
index 64ba254b9..10026c064 100644
--- a/tests/units/test_core.py
+++ b/tests/units/test_core.py
@@ -548,10 +548,8 @@ class TestCore(TestCoreBase):
         "SATELLITE_CARBON_DIOXIDE": ["cop_cds", "dedl", "wekeo_ecmwf"],
         "SATELLITE_FIRE_BURNED_AREA": ["cop_cds", "wekeo_ecmwf"],
         "SATELLITE_METHANE": ["cop_cds", "dedl", "wekeo_ecmwf"],
-        "SATELLITE_SEA_LEVEL_BLACK_SEA": ["dedl", "wekeo_ecmwf"],
         "SATELLITE_SEA_ICE_EDGE_TYPE": ["cop_cds", "dedl", "wekeo_ecmwf"],
         "SATELLITE_SEA_LEVEL_GLOBAL": ["cop_cds", "dedl", "wekeo_ecmwf"],
-        "SATELLITE_SEA_LEVEL_MEDITERRANEAN": ["dedl", "wekeo_ecmwf"],
         "SATELLITE_SEA_ICE_CONCENTRATION": ["cop_cds", "wekeo_ecmwf"],
         "SATELLITE_SEA_ICE_THICKNESS": ["cop_cds", "wekeo_ecmwf"],
         "SEASONAL_MONTHLY_PL": ["cop_cds", "dedl", "wekeo_ecmwf"],