disdrodb 0.1.4__py3-none-any.whl → 0.2.0__py3-none-any.whl

disdrodb/__init__.py

@@ -39,8 +39,6 @@ from disdrodb.configs import (
     define_configs,
     get_data_archive_dir,
     get_metadata_archive_dir,
-    get_model_options,
-    get_product_options,
     get_scattering_table_dir,
 )
 from disdrodb.data_transfer.download_data import download_archive, download_station
@@ -145,8 +143,6 @@ __all__ = [
     "generate_l2m",
     "get_data_archive_dir",
     "get_metadata_archive_dir",
-    "get_model_options",
-    "get_product_options",
     "get_reader",
     "get_scattering_table_dir",
     "get_station_reader",
@@ -177,7 +173,7 @@ __all__ = [
 ]
 
 
-__root_path__ = os.path.dirname(os.path.dirname(os.path.realpath(__file__)))
+package_dir = os.path.dirname(os.path.realpath(__file__))
 
 
 def is_pytmatrix_available():

disdrodb/_version.py

@@ -28,7 +28,7 @@ version_tuple: VERSION_TUPLE
 commit_id: COMMIT_ID
 __commit_id__: COMMIT_ID
 
-__version__ = version = '0.1.4'
-__version_tuple__ = version_tuple = (0, 1, 4)
+__version__ = version = '0.2.0'
+__version_tuple__ = version_tuple = (0, 2, 0)
 
 __commit_id__ = commit_id = None

disdrodb/accessor/methods.py

@@ -89,6 +89,12 @@ class DISDRODB_Base_Accessor:
 
         return align(self._obj, *args)
 
+    def plot_spectrum(self, **kwargs):
+        """Plot spectrum."""
+        from disdrodb.viz.plots import plot_spectrum
+
+        plot_spectrum(self._obj, **kwargs)
+
 
 @xr.register_dataset_accessor("disdrodb")
 class DISDRODB_Dataset_Accessor(DISDRODB_Base_Accessor):
@@ -97,7 +103,7 @@ class DISDRODB_Dataset_Accessor(DISDRODB_Base_Accessor):
     def __init__(self, xarray_obj):
         super().__init__(xarray_obj)
 
-    def resample(self, accumulation_interval, rolling=True):
+    def resample(self, temporal_resolution):
         """Resample a L1 or L2 DISDRODB Product."""
         from disdrodb.l1.resampling import resample_dataset
 
@@ -105,8 +111,7 @@ class DISDRODB_Dataset_Accessor(DISDRODB_Base_Accessor):
         ds = resample_dataset(
             self._obj,
             sample_interval=sample_interval,
-            accumulation_interval=accumulation_interval,
-            rolling=rolling,
+            temporal_resolution=temporal_resolution,
         )
         return ds
 
@@ -116,6 +121,12 @@ class DISDRODB_Dataset_Accessor(DISDRODB_Base_Accessor):
 
         plot_nd(self._obj, **kwargs)
 
+    def plot_raw_and_filtered_spectra(self, **kwargs):
+        """Plot the raw and filtered spectra."""
+        from disdrodb.viz.plots import plot_raw_and_filtered_spectra
+
+        plot_raw_and_filtered_spectra(self._obj, **kwargs)
+
 
 @xr.register_dataarray_accessor("disdrodb")
 class DISDRODB_DataArray_Accessor(DISDRODB_Base_Accessor):

disdrodb/api/checks.py

@@ -154,6 +154,16 @@ def check_rolling(rolling):
         raise TypeError("'rolling' must be a boolean.")
 
 
+def check_temporal_resolution(temporal_resolution):
+    """Check temporal resolution validity."""
+    from disdrodb.utils.time import get_sampling_information
+
+    if not isinstance(temporal_resolution, str):
+        raise TypeError("'temporal_resolution' must be a string.")
+    # If correct, the follow should not raise error
+    sample_interval, rolling = get_sampling_information(temporal_resolution)
+
+
 def check_folder_partitioning(folder_partitioning):
     """
     Check if the given folder partitioning scheme is valid.

disdrodb/api/create_directories.py

@@ -85,8 +85,6 @@ def create_l0_directory_structure(
     ``product = "L0A"`` will call ``run_l0a``.
     ``product = "L0B"`` will call ``run_l0b_nc``.
     """
-    from disdrodb.configs import get_data_archive_dir, get_metadata_archive_dir
-
     # Retrieve the DISDRODB Metadata Archive directory
     data_archive_dir = get_data_archive_dir(data_archive_dir)
     metadata_archive_dir = get_metadata_archive_dir(metadata_archive_dir)

disdrodb/api/io.py

@@ -40,6 +40,7 @@ from disdrodb.api.path import (
     define_station_dir,
 )
 from disdrodb.l0.l0_reader import define_readers_directory
+from disdrodb.utils.dict import extract_product_kwargs
 from disdrodb.utils.directories import list_files
 from disdrodb.utils.logger import (
     log_info,
@@ -117,7 +118,7 @@ def filter_by_time(filepaths, start_time=None, end_time=None):
     start_time, end_time = check_start_end_time(start_time, end_time)
 
     # -------------------------------------------------------------------------.
-    # - Retrieve start_time and end_time of GPM granules
+    # - Retrieve files start_time and end_time
     l_start_time, l_end_time = get_start_end_time_from_filepaths(filepaths)
 
     # -------------------------------------------------------------------------.
@@ -165,12 +166,9 @@ def find_files(
 
     Other Parameters
     ----------------
-    sample_interval : int, optional
-        The sampling interval in seconds of the product.
-        It must be specified only for product L2E and L2M !
-    rolling : bool, optional
-        Whether the dataset has been resampled by aggregating or rolling.
-        It must be specified only for product L2E and L2M !
+    temporal_resolution : str, optional
+        The temporal resolution of the product (e.g., "1MIN", "10MIN", "1H").
+        It must be specified only for product L1, L2E and L2M !
     model_name : str
         The model name of the statistical distribution for the DSD.
         It must be specified only for product L2M !
@@ -411,21 +409,19 @@ def open_dataset(
         The name of the station.
     product : str
         The name DISDRODB product.
-    sample_interval : int, optional
-        The sampling interval in seconds of the product.
-        It must be specified only for product L2E and L2M !
-    rolling : bool, optional
-        Whether the dataset has been resampled by aggregating or rolling.
-        It must be specified only for product L2E and L2M !
-    model_name : str
-        The model name of the statistical distribution for the DSD.
-        It must be specified only for product L2M !
     debugging_mode : bool, optional
         If ``True``, it select maximum 3 files for debugging purposes.
         The default value is ``False``.
     data_archive_dir : str, optional
         The base directory of DISDRODB, expected in the format ``<...>/DISDRODB``.
         If not specified, the path specified in the DISDRODB active configuration will be used.
+    **product_kwargs : optional
+        DISDRODB product options
+        It must be specified only for product L1, L2E and L2M products !
+        For L1, L2E and L2M products, temporal_resolution is required
+        FOr L2M product, model_name is required.
+    **open_kwargs : optional
+        Additional keyword arguments passed to ``xarray.open_mfdataset()``.
 
     Returns
     -------
@@ -434,7 +430,8 @@ def open_dataset(
     """
     from disdrodb.l0.l0a_processing import read_l0a_dataframe
 
-    product_kwargs = product_kwargs if product_kwargs else {}
+    # Extract product kwargs from open_kwargs
+    product_kwargs = extract_product_kwargs(open_kwargs, product=product)
 
     # List product files
     filepaths = find_files(

disdrodb/api/path.py

@@ -336,13 +336,13 @@ def define_issue_filepath(
 
 def define_config_dir(product):
     """Define the config directory path of a given DISDRODB product."""
-    from disdrodb import __root_path__
+    from disdrodb import package_dir
 
     if product.upper() in ["RAW", "L0A", "L0B"]:
         dir_name = "l0"
     else:
         raise NotImplementedError(f"Product {product} not implemented.")
-    config_dir = os.path.join(__root_path__, "disdrodb", dir_name, "configs")
+    config_dir = os.path.join(package_dir, dir_name, "configs")
     return config_dir
 
 
@@ -448,12 +448,9 @@ def define_product_dir_tree(
     ----------
     product : str
        The DISDRODB product. See ``disdrodb.available_products()``.
-    sample_interval : int, optional
-        The sampling interval in seconds of the product.
-        It must be specified only for product L2E and L2M !
-    rolling : bool, optional
-        Whether the dataset has been resampled by aggregating or rolling.
-        It must be specified only for product L2E and L2M !
+    temporal_resolution : str, optional
+        The temporal resolution of the product.
+        It must be specified only for product L1, L2E and L2M !
     model_name : str
         The custom model name of the fitted statistical distribution.
         It must be specified only for product L2M !
@@ -463,28 +460,23 @@ def define_product_dir_tree(
     data_dir : str
         Station data directory path
     """
-    from disdrodb.api.checks import check_product, check_product_kwargs, check_rolling, check_sample_interval
+    from disdrodb.api.checks import check_product, check_product_kwargs, check_temporal_resolution
 
     product = check_product(product)
     product_kwargs = check_product_kwargs(product, product_kwargs)
     if product.upper() == "RAW":
         return ""
-    if product.upper() in ["L0A", "L0B", "L0C", "L1"]:
+    if product.upper() in ["L0A", "L0B", "L0C"]:
         return ""
-    if product == "L2E":
-        rolling = product_kwargs.get("rolling")
-        sample_interval = product_kwargs.get("sample_interval")
-        check_rolling(rolling)
-        check_sample_interval(sample_interval)
-        temporal_resolution = define_temporal_resolution(seconds=sample_interval, rolling=rolling)
-        return os.path.join(temporal_resolution)
+    if product in ["L1", "L2E"]:
+        temporal_resolution = product_kwargs.get("temporal_resolution")
+        check_temporal_resolution(temporal_resolution)
+        return temporal_resolution
     # L2M if product == "L2M":
-    rolling = product_kwargs.get("rolling")
-    sample_interval = product_kwargs.get("sample_interval")
+    temporal_resolution = product_kwargs.get("temporal_resolution")
+    check_temporal_resolution(temporal_resolution)
     model_name = product_kwargs.get("model_name")
-    check_rolling(rolling)
-    check_sample_interval(sample_interval)
-    temporal_resolution = define_temporal_resolution(seconds=sample_interval, rolling=rolling)
+
     return os.path.join(model_name, temporal_resolution)
 
 
@@ -611,12 +603,9 @@ def define_data_dir(
         If not specified, the path specified in the DISDRODB active configuration will be used.
     check_exists : bool, optional
         Whether to check if the directory exists. The default value is ``False``.
-    sample_interval : int, optional
-        The sampling interval in seconds of the product.
-        It must be specified only for product L2E and L2M !
-    rolling : bool, optional
-        Whether the dataset has been resampled by aggregating or rolling.
-        It must be specified only for product L2E and L2M !
+    temporal_resolution : str, optional
+        The temporal resolution of the product.
+        It must be specified only for product L1, L2E and L2M !
     model_name : str
         The name of the fitted statistical distribution for the DSD.
         It must be specified only for product L2M !
@@ -701,11 +690,8 @@ def define_filename(
     end_time : datetime.datatime, optional
         End time.
         Required if add_time_period = True.
-    sample_interval : int, optional
-        The sampling interval in seconds of the product.
-        It must be specified only for product L0C, L1, L2E and L2M !
-    rolling : bool, optional
-        Whether the dataset has been resampled by aggregating or rolling.
+    temporal_resolution : str, optional
+        The temporal resolution of the product.
         It must be specified only for product L1, L2E and L2M !
     model_name : str
         The model name of the fitted statistical distribution for the DSD.
@@ -716,7 +702,7 @@ def define_filename(
     str
         L0B file name.
     """
-    from disdrodb.api.checks import check_product, check_product_kwargs
+    from disdrodb.api.checks import check_product, check_product_kwargs, check_temporal_resolution
 
     product = check_product(product)
     product_kwargs = check_product_kwargs(product, product_kwargs)
@@ -729,16 +715,9 @@ def define_filename(
     product_name = f"{product}"
 
     # L0C ... sample interval known only per-file
-    # L1 ... in future known a priori
-    # if product in ["L1"]:
-    #     # TODO: HACK FOR CURRENT L0C and L1 log files in create_product_logs
-    #     sample_interval = product_kwargs.get("sample_interval",  0)
-    #     temporal_resolution = define_temporal_resolution(seconds=sample_interval, rolling=False)
-    #     product_name = f"{product}.{temporal_resolution}"
-    if product in ["L2E", "L2M"]:
-        rolling = product_kwargs.get("rolling")
-        sample_interval = product_kwargs.get("sample_interval")
-        temporal_resolution = define_temporal_resolution(seconds=sample_interval, rolling=rolling)
+    if product in ["L1", "L2E", "L2M"]:
+        temporal_resolution = product_kwargs.get("temporal_resolution")
+        check_temporal_resolution(temporal_resolution)
         product_name = f"{product}.{temporal_resolution}"
     if product in ["L2M"]:
         model_name = product_kwargs.get("model_name")
@@ -828,7 +807,7 @@ def define_l0b_filename(ds, campaign_name: str, station_name: str) -> str:
 
 def define_l0c_filename(ds, campaign_name: str, station_name: str) -> str:
     """Define L0C file name."""
-    # TODO: add sample_interval as function argument    s
+    # TODO: add sample_interval as function argument
     sample_interval = int(ensure_sample_interval_in_seconds(ds["sample_interval"]).data.item())
     temporal_resolution = define_temporal_resolution(sample_interval, rolling=False)
     starting_time, ending_time = get_file_start_end_time(ds)
@@ -839,37 +818,26 @@ def define_l0c_filename(ds, campaign_name: str, station_name: str) -> str:
     return filename
 
 
-def define_l1_filename(ds, campaign_name, station_name: str) -> str:
+def define_l1_filename(ds, campaign_name, station_name: str, temporal_resolution: str) -> str:
     """Define L1 file name."""
-    # TODO: add sample_interval and rolling as function argument
-
     starting_time, ending_time = get_file_start_end_time(ds)
-    sample_interval = int(ensure_sample_interval_in_seconds(ds["sample_interval"]).data.item())
-    temporal_resolution = define_temporal_resolution(sample_interval, rolling=False)
-    starting_time, ending_time = get_file_start_end_time(ds)
-    starting_time = starting_time.strftime("%Y%m%d%H%M%S")
-    ending_time = ending_time.strftime("%Y%m%d%H%M%S")
-    version = ARCHIVE_VERSION
-    filename = f"L1.{temporal_resolution}.{campaign_name}.{station_name}.s{starting_time}.e{ending_time}.{version}.nc"
-
-    # filename = define_filename(
-    #     product="L1",
-    #     campaign_name=campaign_name,
-    #     station_name=station_name,
-    #     # Filename options
-    #     start_time=starting_time,
-    #     end_time=ending_time,
-    #     add_version=True,
-    #     add_time_period=True,
-    #     add_extension=True,
-    #     # Product options
-    #     # sample_interval=sample_interval,
-    #     # rolling=rolling,
-    # )
+    filename = define_filename(
+        product="L1",
+        campaign_name=campaign_name,
+        station_name=station_name,
+        # Filename options
+        start_time=starting_time,
+        end_time=ending_time,
+        add_version=True,
+        add_time_period=True,
+        add_extension=True,
+        # Product options
+        temporal_resolution=temporal_resolution,
+    )
     return filename
 
 
-def define_l2e_filename(ds, campaign_name: str, station_name: str, sample_interval: int, rolling: bool) -> str:
+def define_l2e_filename(ds, campaign_name: str, station_name: str, temporal_resolution: str) -> str:
     """Define L2E file name."""
     starting_time, ending_time = get_file_start_end_time(ds)
     filename = define_filename(
@@ -883,8 +851,7 @@ def define_l2e_filename(ds, campaign_name: str, station_name: str, sample_interv
         add_time_period=True,
         add_extension=True,
         # Product options
-        sample_interval=sample_interval,
-        rolling=rolling,
+        temporal_resolution=temporal_resolution,
     )
     return filename
 
@@ -893,8 +860,7 @@ def define_l2m_filename(
     ds,
     campaign_name: str,
     station_name: str,
-    sample_interval: int,
-    rolling: bool,
+    temporal_resolution: str,
     model_name: str,
 ) -> str:
     """Define L2M file name."""
@@ -910,8 +876,7 @@ def define_l2m_filename(
         add_time_period=True,
         add_extension=True,
         # Product options
-        sample_interval=sample_interval,
-        rolling=rolling,
+        temporal_resolution=temporal_resolution,
         model_name=model_name,
     )
     return filename

disdrodb/api/search.py

@@ -17,6 +17,7 @@ from disdrodb.api.path import (
 )
 from disdrodb.configs import get_data_archive_dir, get_metadata_archive_dir
 from disdrodb.constants import PRODUCTS_REQUIREMENTS
+from disdrodb.utils.dict import extract_product_kwargs
 from disdrodb.utils.directories import contains_files, contains_netcdf_or_parquet_files, list_directories, list_files
 from disdrodb.utils.yaml import read_yaml
 
@@ -173,7 +174,15 @@ def list_station_names(
 #### Filtering utilities for available_stations
 
 
-def _finalize_output(list_info, return_tuple):
+def _finalize_output(list_info, return_tuple, metadata_archive_dir, filter_kwargs):
+    # Filter stations if metadata filtering values are specified
+    if len(filter_kwargs) != 0:
+        list_info = select_stations_matching_metadata_values(
+            metadata_archive_dir=metadata_archive_dir,
+            list_info=list_info,
+            filter_kwargs=filter_kwargs,
+        )
+
     # - Return the (data_source, campaign_name, station_name) tuple
     if return_tuple:
         return list_info
@@ -192,7 +201,7 @@ def is_disdrodb_data_url_specified(metadata_filepath):
     return isinstance(disdrodb_data_url, str) and len(disdrodb_data_url) > 1
 
 
-def keep_list_info_with_disdrodb_data_url(metadata_archive_dir, list_info):
+def select_stations_with_disdrodb_data_url(metadata_archive_dir, list_info):
     """Keep only the stations with disdrodb_data_url specified in the metadata file."""
     list_info_with_data = []
     for data_source, campaign_name, station_name in list_info:
@@ -209,7 +218,41 @@ def keep_list_info_with_disdrodb_data_url(metadata_archive_dir, list_info):
     return list_info_with_data
 
 
-def keep_list_info_elements_with_product_directory(data_archive_dir, product, list_info):
+def _matches(metadata_value, expected_value):
+    """Return True if metadata_value matches expected_value."""
+    # Case 1: both lists → check any intersection
+    if isinstance(metadata_value, list) and isinstance(expected_value, list):
+        return any(v in metadata_value for v in expected_value)
+    # Case 2: metadata is list → check membership
+    if isinstance(metadata_value, list):
+        return expected_value in metadata_value
+    # Case 3: expected is list → check membership
+    if isinstance(expected_value, list):
+        return metadata_value in expected_value
+    # Case 4: both scalars → direct equality
+    return metadata_value == expected_value
+
+
+def select_stations_matching_metadata_values(metadata_archive_dir, list_info, filter_kwargs):
+    """Keep only the stations with the specified metadata key matching the specified value."""
+    list_info_valid = []
+    for data_source, campaign_name, station_name in list_info:
+        # Define metadata filepath
+        metadata_filepath = define_metadata_filepath(
+            metadata_archive_dir=metadata_archive_dir,
+            data_source=data_source,
+            campaign_name=campaign_name,
+            station_name=station_name,
+        )
+        # Read metadata
+        metadata = read_yaml(metadata_filepath)
+        if np.all([_matches(metadata.get(k), v) for k, v in filter_kwargs.items()]):
+            list_info_valid.append((data_source, campaign_name, station_name))
+
+    return list_info_valid
+
+
+def select_stations_with_product_directory(data_archive_dir, product, list_info):
     """Keep only the stations with the product directory."""
     list_info_with_product_directory = []
     for data_source, campaign_name, station_name in list_info:
@@ -228,7 +271,7 @@ def keep_list_info_elements_with_product_directory(data_archive_dir, product, li
     return list_info_with_product_directory
 
 
-def keep_list_info_elements_with_product_data(data_archive_dir, product, list_info, **product_kwargs):
+def select_stations_with_product_data(data_archive_dir, product, list_info, **product_kwargs):
     """Keep only the stations with product data."""
     # Define file checking function
     checking_function = contains_files if product == "RAW" else contains_netcdf_or_parquet_files
@@ -265,7 +308,7 @@ def available_stations(
     invalid_fields_policy="raise",
     data_archive_dir=None,
     metadata_archive_dir=None,
-    **product_kwargs,
+    **filter_kwargs,
 ):
     """
     Return stations information for which metadata or product data are available on disk.
@@ -289,7 +332,7 @@ def available_stations(
 
         If the DISDRODB product is specified,
         it lists the stations present in the local DISDRODB Data Archive given the specified filtering criteria.
-        The default is is None.
+        The default is None.
 
     data_sources : str or sequence of str, optional
         One or more data source identifiers to filter stations by.
@@ -342,15 +385,15 @@ def available_stations(
         If None, the default metadata base directory is used. Default is None.
     **product_kwargs : dict, optional
         Additional arguments required for some products.
-        For example, for the "L2E" product, you need to specify ``rolling`` and
-        ``sample_interval``. For the "L2M" product, you need to specify also
-        the ``model_name``.
+        It must be specified only for product L1, L2E and L2M products !
+        For L1, L2E and L2M products, ``temporal_resolution`` is required.
+        FOr L2M product, ``model_name`` is required.
 
     Returns
     -------
     list
         If ``return_tuple=True``, return a list of tuples ``(data_source, campaign_name, station_name)``.
-        If ``return_tuple=True``,, return a list of station names.
+        If ``return_tuple=True``, return a list of station names.
 
     Examples
     --------
@@ -361,15 +404,18 @@ def available_stations(
     >>> # List stations with raw data available in the local DISDRODB Data Archive
     >>> raw_stations = available_stations(product="RAW", available_data=True)
     >>> # List stations of specific data sources
-    >>> stations = available_stations(data_sources=["GPM", "EPFL"])
+    >>> stations = available_stations(data_sources=["NASA", "EPFL"])
     """
     # Retrieve DISDRODB Data and Metadata Archive directories
     metadata_archive_dir = get_metadata_archive_dir(metadata_archive_dir)
     product = check_product(product) if product is not None else None
     invalid_fields_policy = check_invalid_fields_policy(invalid_fields_policy)
 
+    # Extract product_kwargs from filter_kwargs
+    product_kwargs = extract_product_kwargs(filter_kwargs, product=product) if product is not None else {}
+
     # Retrieve available stations from the Metadata Archive
-    # - Raise error if no stations availables !
+    # - Raise error if no stations available !
     list_info = list_station_names(
         metadata_archive_dir=metadata_archive_dir,
         data_sources=data_sources,
@@ -386,24 +432,34 @@ def available_stations(
             raise_error_if_empty=raise_error_if_empty,
             msg="No station available in the DISDRODB Metadata Archive.",
         )
-        return _finalize_output(list_info, return_tuple=return_tuple)
+        return _finalize_output(
+            list_info,
+            return_tuple=return_tuple,
+            metadata_archive_dir=metadata_archive_dir,
+            filter_kwargs=filter_kwargs,
+        )
 
     # Return stations in the Metadata Archive with specified disdrodb_data_url
     if product is None and available_data:
-        list_info = keep_list_info_with_disdrodb_data_url(metadata_archive_dir, list_info)
+        list_info = select_stations_with_disdrodb_data_url(metadata_archive_dir, list_info)
         _raise_an_error_if_no_stations(
             list_info,
             raise_error_if_empty=raise_error_if_empty,
             msg="No station has the disdrodb_data_url specified in the metadata.",
         )
-        return _finalize_output(list_info, return_tuple=return_tuple)
+        return _finalize_output(
+            list_info,
+            return_tuple=return_tuple,
+            metadata_archive_dir=metadata_archive_dir,
+            filter_kwargs=filter_kwargs,
+        )
 
     # If product is specified, select stations available in the local DISDRODB Data Archive
     # - If available_data=False, search for station with the existing product directory (do not check for data)
     data_archive_dir = get_data_archive_dir(data_archive_dir)
     product = check_product(product)
     if not available_data:
-        list_info = keep_list_info_elements_with_product_directory(
+        list_info = select_stations_with_product_directory(
             data_archive_dir=data_archive_dir,
             product=product,
             list_info=list_info,
@@ -413,11 +469,16 @@ def available_stations(
             raise_error_if_empty=raise_error_if_empty,
             msg=f"No station product {product} directory available in the local DISDRODB Data Archive.",
         )
-        return _finalize_output(list_info, return_tuple=return_tuple)
+        return _finalize_output(
+            list_info,
+            return_tuple=return_tuple,
+            metadata_archive_dir=metadata_archive_dir,
+            filter_kwargs=filter_kwargs,
+        )
 
     # - If available_data=True, search for station with product data
     product_kwargs = check_product_kwargs(product, product_kwargs)
-    list_info = keep_list_info_elements_with_product_data(
+    list_info = select_stations_with_product_data(
         data_archive_dir=data_archive_dir,
         product=product,
         list_info=list_info,
@@ -429,7 +490,12 @@ def available_stations(
         raise_error_if_empty=raise_error_if_empty,
         msg=f"No station has {product} {product_kwargs} data available in the local DISDRODB Data Archive.",
     )
-    return _finalize_output(list_info, return_tuple=return_tuple)
+    return _finalize_output(
+        list_info,
+        return_tuple=return_tuple,
+        metadata_archive_dir=metadata_archive_dir,
+        filter_kwargs=filter_kwargs,
+    )
 
 
 def available_data_sources(
@@ -441,7 +507,7 @@ def available_data_sources(
     invalid_fields_policy="raise",
     data_archive_dir=None,
     metadata_archive_dir=None,
-    **product_kwargs,
+    **kwargs,
 ):
     """Return data sources for which stations are available."""
     list_info = available_stations(
@@ -455,7 +521,7 @@ def available_data_sources(
         invalid_fields_policy=invalid_fields_policy,
         data_archive_dir=data_archive_dir,
         metadata_archive_dir=metadata_archive_dir,
-        **product_kwargs,
+        **kwargs,
     )
     data_sources = [info[0] for info in list_info]
     data_sources = np.unique(data_sources).tolist()
@@ -471,7 +537,7 @@ def available_campaigns(
     invalid_fields_policy="raise",
     data_archive_dir=None,
     metadata_archive_dir=None,
-    **product_kwargs,
+    **kwargs,
 ):
     """Return campaigns names for which stations are available."""
     list_info = available_stations(
@@ -485,7 +551,7 @@ def available_campaigns(
         invalid_fields_policy=invalid_fields_policy,
         data_archive_dir=data_archive_dir,
         metadata_archive_dir=metadata_archive_dir,
-        **product_kwargs,
+        **kwargs,
     )
     campaign_names = [info[1] for info in list_info]
     campaign_names = np.unique(campaign_names).tolist()