climdata 0.1.2__py2.py3-none-any.whl → 0.1.3__py2.py3-none-any.whl

climdata/utils/utils_download.py

@@ -1,9 +1,5 @@
 import pandas as pd
 import numpy as np
-from wetterdienst import Settings
-from wetterdienst.provider.dwd.observation import DwdObservationRequest
-import geemap
-import ee
 import geopandas as gpd
 from omegaconf import DictConfig
 import os
@@ -15,8 +11,6 @@ from datetime import datetime, timedelta
 import xarray as xr
 import hydra
 
-import pint
-import pint_pandas
 
 from google.oauth2 import service_account
 from googleapiclient.discovery import build
@@ -34,257 +28,6 @@ import concurrent.futures
 
 warnings.filterwarnings("ignore", category=Warning)
 
-def fetch_dwd_loc(cfg: DictConfig):
-    
-    param_mapping = cfg.mappings
-    provider = cfg.dataset.lower()
-    parameter_key = cfg.weather.parameter
-    # Validate provider and parameter
-    
-    if provider not in param_mapping:
-        raise ValueError(f"Provider '{provider}' not found in parameter map.")
-    if parameter_key not in param_mapping[provider]['variables']:
-        raise ValueError(f"Parameter '{parameter_key}' not defined for provider '{provider}'.")
-
-    param_info = param_mapping[provider]['variables'][parameter_key]
-    resolution = param_info["resolution"]
-    dataset = param_info["dataset"]
-    variable_name = param_info["name"]
-    units = param_info.get("unit", None)
-
-    lat = cfg.location.lat
-    lon = cfg.location.lon
-    distance = cfg.location.buffer_km
-    start_date = cfg.time_range.start_date
-    end_date = cfg.time_range.end_date
-    output_file = cfg.output.filename
-
-    settings = Settings(
-        ts_shape="long",
-        ts_humanize=True,
-        # ts_si_units=False
-        )
-
-    request = DwdObservationRequest(
-        parameters=(resolution, dataset, variable_name),
-        start_date=start_date,
-        end_date=end_date,
-        settings=settings
-    ).filter_by_distance(
-        latlon=(lat, lon),
-        distance=distance,
-        unit="km"
-    )
-
-    df = request.values.all().df.to_pandas()
-    
-    df['date'] = pd.to_datetime(df['date'])
-    df = df.groupby(['date']).agg({
-        'value': 'mean',
-        'station_id': lambda x: x.mode().iloc[0] if not x.mode().empty else None,
-        'resolution': lambda x: x.mode().iloc[0] if not x.mode().empty else None,
-        'dataset': lambda x: x.mode().iloc[0] if not x.mode().empty else None,
-        'parameter': lambda x: x.mode().iloc[0] if not x.mode().empty else None,
-        'quality': lambda x: x.mode().iloc[0] if not x.mode().empty else None,
-    }).reset_index()
-
-    df.set_index("date", inplace=True)
-    df.reset_index(inplace=True)
-    
-    # Standardize column names
-    df = df.rename(columns={
-        "date": "time",
-        "value": "value",
-        "station_id": "frequent_station",
-
-    })
-    df["variable"] = parameter_key
-    df["latitude"] = lat
-    df["longitude"] = lon
-    df['source'] = 'DWD'
-    df['units'] = units
-    df = df[["latitude", "longitude", "time", "source", "variable", "value","units"]]
-    
-    out_dir = hydra.utils.to_absolute_path(cfg.output.out_dir)
-    os.makedirs(out_dir, exist_ok=True)
-    out_path = os.path.join(out_dir, cfg.output.filename)
-
-    df.to_csv(out_path, index=False)
-    print(f"✅ Saved time series to: {out_path}")
-    return df
-def fetch_ee_loc(cfg: DictConfig):
-    ee.Initialize(project='earthengine-462007')
-
-    provider = cfg.dataset.lower()
-    variable_name = cfg.weather.parameter
-    ee_image_collection = cfg.mappings[provider].params.collection
-    
-    # Prepare the image collection
-    sd = cfg.time_range.start_date
-    ed = cfg.time_range.end_date
-    var_name = cfg.mappings[provider].variables[variable_name].name
-    units = cfg.mappings[provider].variables[variable_name].unit
-    if provider=='gddp':
-        model = cfg.mappings[provider].params.model
-        scenario = cfg.mappings[provider].params.scenario
-        dataset = ee.ImageCollection(ee_image_collection)\
-                    .filter(ee.Filter.date(sd, ed))\
-                    .filter(ee.Filter.eq('model', model))\
-                    .filter(ee.Filter.eq('scenario', scenario))
-    elif provider=='era5-land':
-        dataset = ee.ImageCollection(ee_image_collection)\
-                    .filter(ee.Filter.date(sd, ed))
-    else:
-        raise ValueError(f"Provider '{provider}' is not supported for Earth Engine data fetching.")    
-    image_var = dataset.select(var_name)
-
-    
-    lat = cfg.location.lat
-    lon = cfg.location.lon
-    # identifier = cfg.location.id
-    out_dir = cfg.output.out_dir
-    # buffer = cfg.location.buffer_km
-    buffer = None
-    scale = cfg.mappings[provider].params.scale
-    # retry_delay = cfg.download.retry_delay
-
-    os.makedirs(out_dir, exist_ok=True)
-
-    df = pd.DataFrame([{ "lat": lat, "lon": lon, "id": 0}])
-    gdf = gpd.GeoDataFrame(df, geometry=gpd.points_from_xy(df["lon"], df["lat"]), crs="EPSG:4326")
-
-    try:
-        gdf_ee = geemap.gdf_to_ee(gdf)
-
-        # if buffer:
-        #     pixel_values = gdf_ee.map(
-        #         lambda f: f.set('ts', image_var.getRegion(
-        #             f.buffer(buffer*1e3).bounds().geometry(), scale))
-        #     )
-        # else:
-        pixel_values = gdf_ee.map(
-            lambda f: f.set('ts', image_var.getRegion(
-                f.geometry(), scale))
-        )
-
-        pixel_values_info = pixel_values.getInfo()
-
-        for feature in pixel_values_info['features']:
-            data = feature['properties']['ts']
-            data_id = feature['properties']['id']
-
-            if data:
-                columns = data[0]
-                rows = data[1:]
-                df_out = pd.DataFrame(rows, columns=columns)
-
-                out_dir = hydra.utils.to_absolute_path(cfg.output.out_dir)
-                os.makedirs(out_dir, exist_ok=True)
-                out_path = os.path.join(out_dir, cfg.output.filename)
-                
-                df_out["variable"] = variable_name
-                df_out["latitude"] = lat
-                df_out["longitude"] = lon
-                df_out['source'] = provider.upper()
-                df_out['units'] = units
-                df_out['time'] = pd.to_datetime(df_out['time'], unit='ms')
-                df_out.rename(columns={variable_name: 'value'}, inplace=True)
-                df_out = df_out[["latitude", "longitude", "time", "source", "variable", "value","units"]]
-
-                df_out.to_csv(out_path, index=False)
-                print(f"[\u2713] Saved: {out_path}")
-
-                return df_out
-            else:
-                print(f"[!] No data for ID {data_id}")
-
-    except Exception as e:
-        print(f"[\u2717] Error: {e}")
-        # time.sleep(retry_delay)
-        raise RuntimeError("Failed to download data.")
-
-def fetch_ee_loc_mod(cfg: DictConfig):
-    # Initialize Earth Engine
-    ee.Initialize(project='earthengine-462007')
-
-    provider = cfg.dataset.lower()
-    variable_name = cfg.weather.parameter
-    ee_image_collection = cfg.mappings[provider].params.collection
-
-    sd = cfg.time_range.start_date
-    ed = cfg.time_range.end_date
-    var_name = cfg.mappings[provider].variables[variable_name].name
-    units = cfg.mappings[provider].variables[variable_name].unit
-    scale = cfg.mappings[provider].params.scale
-    out_dir = cfg.output.out_dir
-
-    lat = cfg.location.lat
-    lon = cfg.location.lon
-
-    # Handle model/scenario if needed
-    if provider == 'gddp':
-        model = cfg.mappings[provider].params.model
-        scenario = cfg.mappings[provider].params.scenario
-        dataset = ee.ImageCollection(ee_image_collection) \
-            .filter(ee.Filter.date(sd, ed)) \
-            .filter(ee.Filter.eq('model', model)) \
-            .filter(ee.Filter.eq('scenario', scenario))
-    elif provider == 'era5-land':
-        dataset = ee.ImageCollection(ee_image_collection) \
-            .filter(ee.Filter.date(sd, ed))
-    else:
-        raise ValueError(f"Provider '{provider}' is not supported.")
-
-    image_var = dataset.select(var_name)
-    point = ee.Geometry.Point(lon, lat)
-
-    os.makedirs(out_dir, exist_ok=True)
-    results = []
-
-    print(f"[i] Fetching time series for point: ({lat}, {lon})")
-
-    # Use a client-side list of images
-    image_list = image_var.toList(image_var.size())
-    n_images = image_var.size().getInfo()
-
-    for i in tqdm(range(n_images), desc="Processing images"):
-        try:
-            img = ee.Image(image_list.get(i))
-            date = img.date().format('YYYY-MM-dd').getInfo()
-
-            value = img.reduceRegion(
-                reducer=ee.Reducer.first(),
-                geometry=point,
-                scale=scale,
-                bestEffort=True
-            ).get(var_name)
-
-            value = value.getInfo() if value else None
-            results.append({"date": date, var_name: value})
-        except Exception as e:
-            print(f"[!] Skipping image {i} due to error: {e}")
-            continue
-
-    df_out = pd.DataFrame(results)
-    out_dir = hydra.utils.to_absolute_path(cfg.output.out_dir)
-    os.makedirs(out_dir, exist_ok=True)
-    out_path = os.path.join(out_dir, cfg.output.filename)
-
-    df_out["variable"] = variable_name
-    df_out["latitude"] = lat
-    df_out["longitude"] = lon
-    df_out['units'] = units
-    df_out['source'] = provider.upper()
-    df_out.rename(columns={var_name: 'value', "date": 'time'}, inplace=True)
-    df_out = df_out[["latitude", "longitude", "time", "source", "variable", "value",'units']]
-
-    # ureg = pint.UnitRegistry()
-    # pint_pandas.PintType.ureg = ureg
-    # df_out['temperature'] = df['temperature'].astype('pint[C]')
-    
-    df_out.to_csv(out_path, index=False)
-    print(f"[✓] Saved timeseries to: {out_path}")
-    return df_out
 def list_drive_files(folder_id, service):
     """
     List all files in a Google Drive folder, handling pagination.
@@ -320,79 +63,12 @@ def download_drive_file(file_id, local_path, service):
         while not done:
             status, done = downloader.next_chunk()
             print(f"   → Download {int(status.progress() * 100)}% complete")
-def fetch_MSWX(var_cfg):
-    param_mapping = var_cfg.mappings
-    provider = var_cfg.dataset.lower()
-    parameter_key = var_cfg.weather.parameter
-
-    param_info = param_mapping[provider]['variables'][parameter_key]
-    folder_id = param_info["folder_id"]
-
-    start_date = var_cfg.time_range.start_date
-    end_date = var_cfg.time_range.end_date
-
-    # === 1) Generate expected filenames ===
-    start = datetime.fromisoformat(start_date)
-    end = datetime.fromisoformat(end_date)
-
-    expected_files = []
-    current = start
-    while current <= end:
-        doy = current.timetuple().tm_yday
-        basename = f"{current.year}{doy:03d}.nc"
-        expected_files.append(basename)
-        current += timedelta(days=1)
-
-    output_dir = var_cfg.data_dir
-    local_files = []
-    missing_files = []
-
-    for basename in expected_files:
-        local_path = os.path.join(output_dir, provider, parameter_key, basename)
-        if os.path.exists(local_path):
-            local_files.append(basename)
-        else:
-            missing_files.append(basename)
-
-    if not missing_files:
-        print(f"✅ All {len(expected_files)} files already exist locally. No download needed.")
-        return local_files
-
-    print(f"📂 {len(local_files)} exist, {len(missing_files)} missing — fetching from Drive...")
-
-    # === 2) Connect to Drive ===
-    SCOPES = ['https://www.googleapis.com/auth/drive.readonly']
-    creds = service_account.Credentials.from_service_account_file(
-        param_mapping[provider].params.google_service_account, scopes=SCOPES
-    )
-    service = build('drive', 'v3', credentials=creds)
 
-    # === 3) List all Drive files ===
-    drive_files = list_drive_files(folder_id, service)
-    valid_filenames = set(missing_files)
-
-    files_to_download = [f for f in drive_files if f['name'] in valid_filenames]
-
-    if not files_to_download:
-        print(f"⚠️ None of the missing files found in Drive. Check folder & date range.")
-        return local_files
-
-    # === 4) Download missing ===
-    for file in files_to_download:
-        filename = file['name']
-        local_path = os.path.join(output_dir, provider, parameter_key, filename)
-        print(f"⬇️ Downloading {filename} ...")
-        download_drive_file(file['id'], local_path, service)
-        local_files.append(filename)
-
-    return local_files
-
-
-def fetch_dwd(var_cfg):
+def fetch_dwd(var_cfg,var):
     """Download HYRAS data for one variable and a list of years."""
     param_mapping = var_cfg.mappings
     provider = var_cfg.dataset.lower()
-    parameter_key = var_cfg.weather.parameter
+    parameter_key = var
     # Validate provider and parameter
 
     param_info = param_mapping[provider]['variables'][parameter_key]
@@ -524,452 +200,42 @@ def extract_ts_dwd(cfg: DictConfig):
     ts_all.to_csv(out_path, index=False)
     print(f"✅ Saved time series to: {out_path}")
 
-    return ts_all
-def extract_ts_MSWX(cfg: DictConfig):
-    parameter = cfg.weather.parameter
-    param_mapping = cfg.mappings
-    provider = cfg.dataset.lower()
-    parameter_key = cfg.weather.parameter
-    # Validate provider and parameter
-
-    param_info = param_mapping[provider]['variables'][parameter_key]
-
-    base_dir = cfg.data_dir
-
-    target_lat = cfg.location.lat
-    target_lon = cfg.location.lon
-
-    start_date = pd.to_datetime(cfg.time_range.start_date)
-    end_date = pd.to_datetime(cfg.time_range.end_date)
-
-    # === 1) Rebuild exact basenames ===
-    current = start_date
-    basenames = []
-    while current <= end_date:
-        doy = current.timetuple().tm_yday
-        basename = f"{current.year}{doy:03d}.nc"
-        basenames.append(basename)
-        current += timedelta(days=1)
-
-    # === 2) Process only those files ===
-    ts_list = []
-    missing = []
-
-    for basename in basenames:
-        file_path = os.path.join(base_dir, provider, parameter, basename)
-
-        if not os.path.exists(file_path):
-            missing.append(basename)
-            continue
-
-        print(f"📂 Opening: {file_path}")
-        ds = xr.open_dataset(file_path)
-
-        time_name = [x for x in ds.coords if "time" in x.lower()][0]
-        data_var = [v for v in ds.data_vars][0]
-
-        ts = ds[data_var].sel(
-            lat=target_lat,
-            lon=target_lon,
-            method='nearest'
-        )
-
-        df = ts.to_dataframe().reset_index()[[time_name, data_var]]
-        ts_list.append(df)
-
-    if missing:
-        print(f"⚠️ Warning: {len(missing)} files were missing and skipped:")
-        for m in missing:
-            print(f"   - {m}")
-
-    if not ts_list:
-        raise RuntimeError("❌ No valid files were found. Cannot extract time series.")
-
-    # === 3) Combine and slice (for safety) ===
-    ts_all = pd.concat(ts_list).sort_values(by=time_name).reset_index(drop=True)
-
-    ts_all[time_name] = pd.to_datetime(ts_all[time_name])
-    ts_all = ts_all[
-        (ts_all[time_name] >= start_date) &
-        (ts_all[time_name] <= end_date)
-    ].reset_index(drop=True)
-
-    out_dir = hydra.utils.to_absolute_path(cfg.output.out_dir)
-    os.makedirs(out_dir, exist_ok=True)
-    out_path = os.path.join(out_dir, cfg.output.filename)
-
-    ts_all["variable"] = param_info['name']
-    ts_all["latitude"] = target_lat
-    ts_all["longitude"] = target_lon
-    ts_all['source'] = provider.upper()
-    ts_all['units'] = ts.attrs['units']
-    ts_all.rename(columns={param_info['name']: 'value'}, inplace=True)
-    ts_all = ts_all[["latitude", "longitude", "time", "source", "variable", "value",'units']]
-    
-    ts_all.to_csv(out_path, index=False)
-    print(f"✅ Saved MSWX time series to: {out_path}")
-
     return ts_all
 
 import os
 from omegaconf import DictConfig
 
-def build_output_filename(cfg: DictConfig) -> str:
-    """Generate full output file path from pattern and config."""
-    provider = cfg.dataset.lower()
-    parameter = cfg.weather.parameter
-    lat = cfg.location.lat
-    lon = cfg.location.lon
-    start = cfg.time_range.start_date
-    end = cfg.time_range.end_date
-
-    pattern = cfg.output.get("filename", "{provider}_{parameter}_{start}_{end}.csv")
-    filename = pattern.format(
-        provider=provider,
-        parameter=parameter,
-        lat=lat,
-        lon=lon,
-        start=start,
-        end=end
-    )
-
-    out_dir = cfg.output.out_dir
-    fmt = cfg.output.fmt  # format is a reserved word in Python, so use 'fmt'
-
-    # return os.path.join(out_dir, fmt, filename)
-    return filename
-
-# SPDX-FileCopyrightText: Copyright (c) 2023 - 2024 NVIDIA CORPORATION & AFFILIATES.
-# SPDX-FileCopyrightText: All rights reserved.
-# SPDX-License-Identifier: Apache-2.0
-#
-# Licensed under the Apache License, Version 2.0 (the "License");
-# you may not use this file except in compliance with the License.
-# You may obtain a copy of the License at
-#
-#     http://www.apache.org/licenses/LICENSE-2.0
-#
-# Unless required by applicable law or agreed to in writing, software
-# distributed under the License is distributed on an "AS IS" BASIS,
-# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-# See the License for the specific language governing permissions and
-# limitations under the License.
-'''
-import os
-import tempfile
-import cdsapi
-import xarray as xr
-import datetime
-import json
-import dask
-import calendar
-from dask.diagnostics import ProgressBar
-from typing import List, Tuple, Dict, Union
-import urllib3
-import logging
-import numpy as np
-import fsspec
-
-urllib3.disable_warnings(urllib3.exceptions.InsecureRequestWarning)
-
-
-class ERA5Mirror:
+def get_output_filename(cfg, output_type="nc", lat=None, lon=None):
     """
-    A class to manage downloading ERA5 datasets. The datasets are downloaded from the Copernicus Climate Data Store (CDS) and stored in Zarr format.
-
-    Attributes
-    ----------
-    base_path : Path
-        The path to the Zarr dataset.
-    fs : fsspec.AbstractFileSystem
-        The filesystem to use for the Zarr dataset. If None, the local filesystem will be used.
+    Generate output filename based on config, output type, and extraction mode.
+    output_type: "nc", "csv", or "zarr"
     """
-
-    def __init__(self, base_path: str, fs: fsspec.AbstractFileSystem = None):
-        # Get parameters
-        self.base_path = base_path
-        if fs is None:
-            fs = fsspec.filesystem("file")
-        self.fs = fs
-
-        # Create the base path if it doesn't exist
-        if not self.fs.exists(self.base_path):
-            self.fs.makedirs(self.base_path)
-
-        # Create metadata that will be used to track which chunks have been downloaded
-        self.metadata_file = os.path.join(self.base_path, "metadata.json")
-        self.metadata = self.get_metadata()
-
-    def get_metadata(self):
-        """Get metadata"""
-        if self.fs.exists(self.metadata_file):
-            with self.fs.open(self.metadata_file, "r") as f:
-                try:
-                    metadata = json.load(f)
-                except json.decoder.JSONDecodeError:
-                    metadata = {"chunks": []}
-        else:
-            metadata = {"chunks": []}
-        return metadata
-
-    def save_metadata(self):
-        """Save metadata"""
-        with self.fs.open(self.metadata_file, "w") as f:
-            json.dump(self.metadata, f)
-
-    def chunk_exists(self, variable, year, month, hours, pressure_level):
-        """Check if chunk exists"""
-        for chunk in self.metadata["chunks"]:
-            if (
-                chunk["variable"] == variable
-                and chunk["year"] == year
-                and chunk["month"] == month
-                and chunk["hours"] == hours
-                and chunk["pressure_level"] == pressure_level
-            ):
-                return True
-        return False
-
-    def download_chunk(
-        self,
-        variable: str,
-        year: int,
-        month: int,
-        hours: List[int],
-        pressure_level: int = None,
-    ):
-        """
-        Download ERA5 data for the specified variable, date range, hours, and pressure levels.
-
-        Parameters
-        ----------
-        variable : str
-            The ERA5 variable to download, e.g. 'tisr' for solar radiation or 'z' for geopotential.
-        year : int
-            The year to download.
-        month : int
-            The month to download.
-        hours : List[int]
-            A list of hours (0-23) for which data should be downloaded.
-        pressure_level : int, optional
-            A pressure level to include in the download, by default None. If None, the single-level data will be downloaded.
-
-        Returns
-        -------
-        xr.Dataset
-            An xarray Dataset containing the downloaded data.
-        """
-
-        with tempfile.TemporaryDirectory() as tmpdir:
-            # Get all days in the month
-            days_in_month = calendar.monthrange(year, month)[1]
-
-            # Make tmpfile to store the data
-            output_file = os.path.join(
-                tmpdir,
-                f"{variable}_{year}_{month:02d}_{str(hours)}_{str(pressure_level)}.nc",
-            )
-
-            # start the CDS API client (maybe need to move this outside the loop?)
-            c = cdsapi.Client(quiet=True)
-
-            # Setup the request parameters
-            request_params = {
-                "product_type": "reanalysis",
-                "variable": variable,
-                "year": str(year),
-                "month": str(month),
-                "day": [f"{day:02d}" for day in range(1, days_in_month + 1)],
-                "time": [f"{hour:02d}:00" for hour in hours],
-                "format": "netcdf",
-            }
-            if pressure_level:
-                request_params["pressure_level"] = [str(pressure_level)]
-                dataset_name = "reanalysis-era5-pressure-levels"
-            else:
-                dataset_name = "reanalysis-era5-single-levels"
-
-            # Download the data
-            c.retrieve(
-                dataset_name,
-                request_params,
-                output_file,
-            )
-
-            # Open the downloaded data
-            ds = xr.open_dataset(output_file)
-        return ds
-
-    def variable_to_zarr_name(self, variable: str, pressure_level: int = None):
-        """convert variable to zarr name"""
-        # create zarr path for variable
-        zarr_path = f"{self.base_path}/{variable}"
-        if pressure_level:
-            zarr_path += f"_pressure_level_{pressure_level}"
-        zarr_path += ".zarr"
-        return zarr_path
-
-    def download_and_upload_chunk(
-        self,
-        variable: str,
-        year: int,
-        month: int,
-        hours: List[int],
-        pressure_level: int = None,
-    ):
-        """
-        Downloads a chunk of ERA5 data for a specific variable and date range, and uploads it to a Zarr array.
-        This downloads a 1-month chunk of data.
-
-        Parameters
-        ----------
-        variable : str
-            The variable to download.
-        year : int
-            The year to download.
-        month : int
-            The month to download.
-        hours : List[int]
-            A list of hours to download.
-        pressure_level : int, optional
-            Pressure levels to download, if applicable.
-        """
-
-        # Download the data
-        ds = self.download_chunk(variable, year, month, hours, pressure_level)
-        if "valid_time" in ds.dims:
-            ds = ds.rename({"valid_time": "time"})
-
-        # Create the Zarr path
-        zarr_path = self.variable_to_zarr_name(variable, pressure_level)
-
-        # Specify the chunking options
-        chunking = {"time": 1, "latitude": 721, "longitude": 1440}
-        if "level" in ds.dims:
-            chunking["level"] = 1
-
-        # Re-chunk the dataset
-        ds = ds.chunk(chunking)
-
-        # Check if the Zarr dataset exists
-        if self.fs.exists(zarr_path):
-            mode = "a"
-            append_dim = "time"
-            create = False
-        else:
-            mode = "w"
-            append_dim = None
-            create = True
-
-        # Upload the data to the Zarr dataset
-        mapper = self.fs.get_mapper(zarr_path, create=create)
-        ds.to_zarr(mapper, mode=mode, consolidated=True, append_dim=append_dim)
-
-        # Update the metadata
-        self.metadata["chunks"].append(
-            {
-                "variable": variable,
-                "year": year,
-                "month": month,
-                "hours": hours,
-                "pressure_level": pressure_level,
-            }
+    if output_type == "csv":
+        template = cfg.output.filename_csv
+    elif output_type == "zarr":
+        template = cfg.output.filename_zarr
+    else:
+        template = cfg.output.filename_nc
+
+    # If lat/lon are provided, use point template
+    if lat is not None and lon is not None:
+        filename = template.format(
+            provider=cfg.dataset,
+            parameter="surface",
+            lat=f"{lat}",
+            lon=f"{lon}",
+            start=cfg.time_range.start_date.replace("-", ""),
+            end=cfg.time_range.end_date.replace("-", ""),
+        )
+    else:
+        # Use region bounds
+        region_bounds = cfg.bounds[cfg.region]
+        filename = template.format(
+            provider=cfg.dataset,
+            parameter="surface",
+            lat_range=f"{region_bounds['lat_min']}-{region_bounds['lat_max']}",
+            lon_range=f"{region_bounds['lon_min']}-{region_bounds['lon_max']}",
+            start=cfg.time_range.start_date.replace("-", ""),
+            end=cfg.time_range.end_date.replace("-", ""),
         )
-        self.save_metadata()
-
-    def download(
-        self,
-        variables: List[Union[str, Tuple[str, int]]],
-        date_range: Tuple[datetime.date, datetime.date],
-        hours: List[int],
-    ):
-        """
-        Start the process of mirroring the specified ERA5 variables for the given date range and hours.
-
-        Parameters
-        ----------
-        variables : List[Union[str, Tuple[str, List[int]]]]
-            A list of variables to mirror, where each element can either be a string (single-level variable)
-            or a tuple (variable with pressure level).
-        date_range : Tuple[datetime.date, datetime.date]
-            A tuple containing the start and end dates for the data to be mirrored. This will download and store every month in the range.
-        hours : List[int]
-            A list of hours for which to download the data.
-
-        Returns
-        -------
-        zarr_paths : List[str]
-            A list of Zarr paths for each of the variables.
-        """
-
-        start_date, end_date = date_range
-
-        # Reformat the variables list so all elements are tuples
-        reformated_variables = []
-        for variable in variables:
-            if isinstance(variable, str):
-                reformated_variables.append(tuple([variable, None]))
-            else:
-                reformated_variables.append(variable)
-
-        # Start Downloading
-        with ProgressBar():
-            # Round dates to months
-            current_date = start_date.replace(day=1)
-            end_date = end_date.replace(day=1)
-
-            while current_date <= end_date:
-                # Create a list of tasks to download the data
-                tasks = []
-                for variable, pressure_level in reformated_variables:
-                    if not self.chunk_exists(
-                        variable,
-                        current_date.year,
-                        current_date.month,
-                        hours,
-                        pressure_level,
-                    ):
-                        task = dask.delayed(self.download_and_upload_chunk)(
-                            variable,
-                            current_date.year,
-                            current_date.month,
-                            hours,
-                            pressure_level,
-                        )
-                        tasks.append(task)
-                    else:
-                        print(
-                            f"Chunk for {variable} {pressure_level} {current_date.year}-{current_date.month} already exists. Skipping."
-                        )
-
-                # Execute the tasks with Dask
-                print(f"Downloading data for {current_date.year}-{current_date.month}")
-                if tasks:
-                    dask.compute(*tasks)
-
-                # Update the metadata
-                self.save_metadata()
-
-                # Update the current date
-                days_in_month = calendar.monthrange(
-                    year=current_date.year, month=current_date.month
-                )[1]
-                current_date += datetime.timedelta(days=days_in_month)
-
-        # Return the Zarr paths
-        zarr_paths = []
-        for variable, pressure_level in reformated_variables:
-            zarr_path = self.variable_to_zarr_name(variable, pressure_level)
-            zarr_paths.append(zarr_path)
-
-        # Check that Zarr arrays have correct dt for time dimension
-        for zarr_path in zarr_paths:
-            ds = xr.open_zarr(zarr_path)
-            time_stamps = ds.time.values
-            dt = time_stamps[1:] - time_stamps[:-1]
-            assert np.all(
-                dt == dt[0]
-            ), f"Zarr array {zarr_path} has incorrect dt for time dimension. An error may have occurred during download. Please delete the Zarr array and try again."
-
-        return zarr_paths
-'''
+    return filename