cloudnetpy 1.49.9__py3-none-any.whl → 1.87.3__py3-none-any.whl

cloudnetpy/metadata.py

@@ -4,7 +4,8 @@ from typing import NamedTuple
 
 
 class MetaData(NamedTuple):
-    long_name: str | None = None
+    long_name: str
+    dimensions: tuple[str, ...] | None
     standard_name: str | None = None
     units: str | None = None
     comment: str | None = None
@@ -15,6 +16,7 @@ class MetaData(NamedTuple):
     axis: str | None = None
     calendar: str | None = None
     source: str | None = None
+    source_instrument_pid: str | None = None
 
 
 COMMON_ATTRIBUTES = {
@@ -23,30 +25,38 @@ COMMON_ATTRIBUTES = {
         axis="T",
         standard_name="time",
         calendar="standard",
+        dimensions=("time",),
     ),
     "height": MetaData(
         long_name="Height above mean sea level",
         standard_name="height_above_mean_sea_level",
         units="m",
+        dimensions=("range",),
     ),
     "range": MetaData(
         long_name="Range from instrument",
         axis="Z",
         units="m",
         comment="Distance from instrument to centre of each range bin.",
+        dimensions=("range",),
     ),
     "latitude": MetaData(
         long_name="Latitude of site",
         units="degree_north",
         standard_name="latitude",
+        dimensions=("time",),
     ),
     "longitude": MetaData(
         long_name="Longitude of site",
         units="degree_east",
         standard_name="longitude",
+        dimensions=("time",),
     ),
     "altitude": MetaData(
-        long_name="Altitude of site", standard_name="altitude", units="m"
+        long_name="Altitude of site",
+        standard_name="altitude",
+        units="m",
+        dimensions=("time",),
     ),
     "Zh": MetaData(
         long_name="Radar reflectivity factor",
@@ -54,6 +64,7 @@ COMMON_ATTRIBUTES = {
         comment="Calibrated reflectivity. Calibration convention: in the absence\n"
         "of attenuation, a cloud at 273 K containing one million 100-micron droplets\n"
         "per cubic metre will have a reflectivity of 0 dBZ at all frequencies.",
+        dimensions=("time", "range"),
     ),
     "width": MetaData(
         long_name="Spectral width",
@@ -62,6 +73,7 @@ COMMON_ATTRIBUTES = {
             "This parameter is the standard deviation of the reflectivity-weighted\n"
             "velocities in the radar pulse volume."
         ),
+        dimensions=("time", "range"),
     ),
     "v": MetaData(
         long_name="Doppler velocity",
@@ -70,91 +82,150 @@ COMMON_ATTRIBUTES = {
             "This parameter is the radial component of the velocity, with positive\n"
             "velocities are away from the radar."
         ),
+        dimensions=("time", "range"),
     ),
     "ldr": MetaData(
         long_name="Linear depolarisation ratio",
         units="dB",
         comment="This parameter is the ratio of cross-polar to co-polar reflectivity.",
+        dimensions=("time", "range"),
+    ),
+    "sldr": MetaData(
+        long_name="Slanted linear depolarisation ratio",
+        units="dB",
+        dimensions=("time", "range"),
     ),
-    "sldr": MetaData(long_name="Slanted linear depolarisation ratio", units="dB"),
     "lwp": MetaData(
         long_name="Liquid water path",
         units="kg m-2",
         standard_name="atmosphere_cloud_liquid_water_content",
+        dimensions=("time",),
     ),
     "iwv": MetaData(
         long_name="Integrated water vapour",
         units="kg m-2",
         standard_name="atmosphere_mass_content_of_water_vapor",
+        dimensions=("time",),
     ),
     "kurtosis": MetaData(
         long_name="Kurtosis of spectra",
         units="1",
+        dimensions=("time", "range"),
+    ),
+    "skewness": MetaData(
+        long_name="Skewness of spectra",
+        units="1",
+        dimensions=("time", "range"),
+    ),
+    "nyquist_velocity": MetaData(
+        long_name="Nyquist velocity",
+        units="m s-1",
+        dimensions=("time", "range"),
+    ),
+    "radar_frequency": MetaData(
+        long_name="Radar transmit frequency", units="GHz", dimensions=None
     ),
-    "nyquist_velocity": MetaData(long_name="Nyquist velocity", units="m s-1"),
-    "radar_frequency": MetaData(long_name="Radar transmit frequency", units="GHz"),
     "beta": MetaData(
         long_name="Attenuated backscatter coefficient",
         units="sr-1 m-1",
         comment="SNR-screened attenuated backscatter coefficient.",
+        dimensions=("time", "range"),
     ),
     "beta_raw": MetaData(
         long_name="Attenuated backscatter coefficient",
         units="sr-1 m-1",
         comment="Non-screened attenuated backscatter coefficient.",
+        dimensions=("time", "range"),
     ),
     "beta_smooth": MetaData(
         long_name="Attenuated backscatter coefficient",
         units="sr-1 m-1",
         comment="SNR-screened attenuated backscatter coefficient.\n"
         "Weak background smoothed using Gaussian 2D-kernel.",
+        dimensions=("time", "range"),
     ),
-    "wavelength": MetaData(
-        long_name="Laser wavelength",
-        units="nm",
-    ),
+    "wavelength": MetaData(long_name="Laser wavelength", units="nm", dimensions=None),
     "zenith_angle": MetaData(
         long_name="Zenith angle",
         units="degree",
         standard_name="zenith_angle",
         comment="Angle to the local vertical. A value of zero is directly overhead.",
+        dimensions=("time",),
+    ),
+    "ir_zenith_angle": MetaData(
+        long_name="Infrared sensor zenith angle",
+        units="degree",
+        comment="90=horizon, 0=zenith",
+        dimensions=("time",),
     ),
     "azimuth_angle": MetaData(
         long_name="Azimuth angle",
         standard_name="sensor_azimuth_angle",
         units="degree",
         comment="Angle between North and the line of sight, measured clockwise.",
+        dimensions=("time",),
     ),
-    "temperature": MetaData(
-        long_name="Temperature",
-        units="K",
-    ),
-    "pressure": MetaData(
-        long_name="Pressure",
-        units="Pa",
-    ),
+    "temperature": MetaData(long_name="Temperature", units="K", dimensions=("time",)),
+    "pressure": MetaData(long_name="Pressure", units="Pa", dimensions=("time",)),
     "SNR": MetaData(
-        long_name="Signal-to-noise ratio",
-        units="dB",
+        long_name="Signal-to-noise ratio", units="dB", dimensions=("time", "range")
     ),
     "relative_humidity": MetaData(
         long_name="Relative humidity",
         standard_name="relative_humidity",
         units="1",
+        dimensions=("time",),
+    ),
+    "absolute_humidity": MetaData(
+        long_name="Absolute humidity",
+        standard_name="mass_concentration_of_water_vapor_in_air",
+        units="kg m-3",
+        dimensions=("time", "range"),
     ),
     "wind_speed": MetaData(
         long_name="Wind speed",
         standard_name="wind_speed",
         units="m s-1",
+        dimensions=("time",),
     ),
     "wind_direction": MetaData(
         long_name="Wind direction",
         standard_name="wind_from_direction",
         units="degree",
+        dimensions=("time",),
     ),
     "rainfall_rate": MetaData(
         long_name="Rainfall rate",
         standard_name="rainfall_rate",
         units="m s-1",
+        dimensions=("time",),
+    ),
+    "rainfall_amount": MetaData(
+        long_name="Rainfall amount",
+        standard_name="thickness_of_rainfall_amount",
+        units="m",
+        comment="Cumulated precipitation since 00:00 UTC",
+        dimensions=("time",),
+    ),
+    "air_temperature": MetaData(
+        long_name="Air temperature",
+        standard_name="air_temperature",
+        units="K",
+        dimensions=("time",),
+    ),
+    "air_pressure": MetaData(
+        long_name="Air pressure",
+        standard_name="air_pressure",
+        units="Pa",
+        dimensions=("time",),
+    ),
+    "snr_limit": MetaData(
+        long_name="SNR limit",
+        units="dB",
+        comment="SNR threshold used in data screening.",
+        dimensions=None,
+    ),
+    "synop_WaWa": MetaData(
+        long_name="Synop code WaWa", units="1", dimensions=("time",)
     ),
 }

cloudnetpy/model_evaluation/file_handler.py

@@ -1,7 +1,10 @@
 import os
 from datetime import datetime
+from os import PathLike
+from uuid import UUID
 
 import netCDF4
+import numpy.typing as npt
 
 from cloudnetpy import output, utils
 
@@ -10,12 +13,11 @@ from .metadata import (
     MODEL_ATTRIBUTES,
     MODEL_L3_ATTRIBUTES,
     REGRID_PRODUCT_ATTRIBUTES,
-    MetaData,
 )
 from .products.model_products import ModelManager
 
 
-def update_attributes(model_downsample_variables: dict, attributes: dict):
+def update_attributes(model_downsample_variables: dict, attributes: dict) -> None:
     """Overrides existing Cloudnet-ME Array-attributes.
         Overrides existing attributes using hard-coded values.
         New attributes are added.
@@ -33,11 +35,11 @@ def update_attributes(model_downsample_variables: dict, attributes: dict):
             model_downsample_variables[key].set_attributes(MODEL_ATTRIBUTES[key])
         elif "_".join(key_parts[0:-1]) in REGRID_PRODUCT_ATTRIBUTES:
             model_downsample_variables[key].set_attributes(
-                REGRID_PRODUCT_ATTRIBUTES["_".join(key_parts[0:-1])]
+                REGRID_PRODUCT_ATTRIBUTES["_".join(key_parts[0:-1])],
             )
         elif "_".join(key_parts[0:-2]) in REGRID_PRODUCT_ATTRIBUTES:
             model_downsample_variables[key].set_attributes(
-                REGRID_PRODUCT_ATTRIBUTES["_".join(key_parts[0:-2])]
+                REGRID_PRODUCT_ATTRIBUTES["_".join(key_parts[0:-2])],
             )
         elif (
             "_".join(key_parts[1:]) in MODEL_L3_ATTRIBUTES
@@ -45,29 +47,29 @@ def update_attributes(model_downsample_variables: dict, attributes: dict):
         ):
             try:
                 model_downsample_variables[key].set_attributes(
-                    MODEL_L3_ATTRIBUTES["_".join(key_parts[1:])]
+                    MODEL_L3_ATTRIBUTES["_".join(key_parts[1:])],
                 )
             except KeyError:
                 model_downsample_variables[key].set_attributes(
-                    MODEL_L3_ATTRIBUTES["_".join(key_parts[2:])]
+                    MODEL_L3_ATTRIBUTES["_".join(key_parts[2:])],
                 )
         elif "_".join(key_parts[1:]) in CYCLE_ATTRIBUTES:
             model_downsample_variables[key].set_attributes(
-                CYCLE_ATTRIBUTES["_".join(key_parts[1:])]
+                CYCLE_ATTRIBUTES["_".join(key_parts[1:])],
             )
         elif "_".join(key_parts[2:]) in CYCLE_ATTRIBUTES:
             model_downsample_variables[key].set_attributes(
-                CYCLE_ATTRIBUTES["_".join(key_parts[2:])]
+                CYCLE_ATTRIBUTES["_".join(key_parts[2:])],
             )
 
 
 def save_downsampled_file(
     id_mark: str,
-    file_name: str,
+    file_name: str | PathLike,
     objects: tuple,
-    files: tuple,
-    uuid: str | None,
-):
+    files: tuple[list[str | PathLike], str | PathLike],
+    uuid: UUID,
+) -> None:
     """Saves a standard downsampled day product file.
 
     Args:
@@ -83,35 +85,31 @@ def save_downsampled_file(
     """
     obj = objects[0]
     dimensions = {"time": len(obj.time), "level": len(obj.data["level"][:])}
-    root_group = output.init_file(file_name, dimensions, obj.data, uuid)
-    _augment_global_attributes(root_group)
-    uuid = root_group.file_uuid
-    root_group.cloudnet_file_type = id_mark.split("-")[0]
-    root_group.title = (
-        f"Downsampled {id_mark.capitalize().replace('_', ' of ')} "
-        f"from {obj.dataset.location}"
-    )
-    _add_source(root_group, objects, files)
-    output.copy_global(obj.dataset, root_group, ("location", "day", "month", "year"))
-    try:
-        obj.dataset.day
-    except AttributeError:
-        root_group.year, root_group.month, root_group.day = obj.date
-    output.merge_history(root_group, id_mark, {"l3": obj})
-    root_group.close()
-    return uuid
-
-
-def add_var2ncfile(obj: ModelManager, file_name: str):
-    nc_file = netCDF4.Dataset(file_name, "r+", format="NETCDF4_CLASSIC")
-    _write_vars2nc(nc_file, obj.data)
-    nc_file.close()
-
-
-def _write_vars2nc(rootgrp: netCDF4.Dataset, cloudnet_variables: dict):
+    with output.init_file(file_name, dimensions, obj.data, uuid) as root_group:
+        _augment_global_attributes(root_group)
+        root_group.cloudnet_file_type = "l3-" + id_mark.split("_")[0]
+        root_group.title = (
+            f"Downsampled {id_mark.capitalize().replace('_', ' of ')} "
+            f"from {obj.dataset.location}"
+        )
+        _add_source(root_group, objects, files)
+        output.copy_global(
+            obj.dataset, root_group, ("location", "day", "month", "year")
+        )
+        if not hasattr(obj.dataset, "day"):
+            root_group.year, root_group.month, root_group.day = obj.date
+        output.merge_history(root_group, id_mark, obj)
+
+
+def add_var2ncfile(obj: ModelManager, file_name: str | PathLike) -> None:
+    with netCDF4.Dataset(file_name, "r+", format="NETCDF4_CLASSIC") as nc_file:
+        _write_vars2nc(nc_file, obj.data)
+
+
+def _write_vars2nc(rootgrp: netCDF4.Dataset, cloudnet_variables: dict) -> None:
     """Iterates over Cloudnet-ME instances and write to given rootgrp."""
 
-    def _get_dimensions(array):
+    def _get_dimensions(array: npt.NDArray) -> tuple:
         """Finds correct dimensions for a variable."""
         if utils.isscalar(array):
             return ()
@@ -119,8 +117,8 @@ def _write_vars2nc(rootgrp: netCDF4.Dataset, cloudnet_variables: dict):
         file_dims = rootgrp.dimensions
         array_dims = array.shape
         for length in array_dims:
-            dim = [key for key in file_dims.keys() if file_dims[key].size == length][0]
-            variable_size = variable_size + (dim,)
+            dim = [key for key in file_dims if file_dims[key].size == length][0]  # noqa: RUF015
+            variable_size = (*variable_size, dim)
         return variable_size
 
     for key in cloudnet_variables:
@@ -128,7 +126,10 @@ def _write_vars2nc(rootgrp: netCDF4.Dataset, cloudnet_variables: dict):
         size = _get_dimensions(obj.data)
         try:
             nc_variable = rootgrp.createVariable(
-                obj.name, obj.data_type, size, zlib=True
+                obj.name,
+                obj.data_type,
+                size,
+                zlib=True,
             )
             nc_variable[:] = obj.data
             for attr in obj.fetch_attributes():
@@ -137,13 +138,12 @@ def _write_vars2nc(rootgrp: netCDF4.Dataset, cloudnet_variables: dict):
             continue
 
 
-def _augment_global_attributes(root_group: netCDF4.Dataset):
+def _augment_global_attributes(root_group: netCDF4.Dataset) -> None:
     root_group.Conventions = "CF-1.8"
-    # root_group.cloudnetme_version = version.__version__
 
 
-def _add_source(root_ground: netCDF4.Dataset, objects: tuple, files: tuple):
-    """Generates source info for multiple files"""
+def _add_source(root_ground: netCDF4.Dataset, objects: tuple, files: tuple) -> None:
+    """Generates source info for multiple files."""
     model, obs = objects
     model_files, obs_file = files
     source = f"Observation file: {os.path.basename(obs_file)}"
@@ -154,18 +154,21 @@ def _add_source(root_ground: netCDF4.Dataset, objects: tuple, files: tuple):
         if i < len(model_files) - 1:
             source += "\n"
     root_ground.source = source
-    root_ground.source_file_uuids = output.get_source_uuids(model, obs)
+    root_ground.source_file_uuids = output.get_source_uuids([model, obs])
 
 
 def add_time_attribute(date: datetime) -> dict:
-    """ "Adds time attribute with correct units.
+    """Adds time attribute with correct units.
+
     Args:
         attributes: Attributes of variables.
         date: Date as Y M D 0 0 0.
+
     Returns:
         dict: Same attributes with 'time' attribute added.
     """
-    d = date.strftime("%y.%m.%d")
-    attributes = {}
-    attributes["time"] = MetaData(units=f"hours since {d} 00:00:00", long_name="")
-    return attributes
+    return {
+        "time": MODEL_ATTRIBUTES["time"]._replace(
+            units=f"hours since {date:%Y-%m-%d} 00:00:00 +00:00",
+        ),
+    }

cloudnetpy/model_evaluation/metadata.py

@@ -1,27 +1,24 @@
-from typing import NamedTuple
-
-
-class MetaData(NamedTuple):
-    long_name: str
-    units: str
-    comment: str | None = None
-    standard_name: str | None = None
-    axis: str | None = None
-    positive: str | None = None
-
+from cloudnetpy.metadata import MetaData
 
 MODEL_ATTRIBUTES = {
     "time": MetaData(
-        units="decimal hours since midnight",
+        units="",
         long_name="Time UTC",
         axis="T",
+        calendar="standard",
+        dimensions=("time",),
+    ),
+    "latitude": MetaData(
+        long_name="Latitude of grid point", units="degree_north", dimensions=None
+    ),
+    "longitude": MetaData(
+        long_name="Longitude of grid point", units="degree_east", dimensions=None
     ),
-    "latitude": MetaData(long_name="Latitude of grid point", units="dergees_north"),
-    "longitude": MetaData(long_name="Longitude of grid point", units="degrees_east"),
     "horizontal_resolution": MetaData(
         long_name="Horizontal resolution of model",
         units="km",
         comment="Distance between two grid point",
+        dimensions=None,
     ),
     "level": MetaData(
         long_name="Model level",
@@ -29,6 +26,7 @@ MODEL_ATTRIBUTES = {
         comment="Level 1 describes the highest height from ground.",
         axis="Z",
         positive="down",
+        dimensions=("level",),
     ),
 }
 
@@ -41,6 +39,7 @@ CYCLE_ATTRIBUTES = {
             "which it was taken. Note that the profiles in this file may be taken\n"
             "from more than one forecast."
         ),
+        dimensions=("time",),
     ),
     "height": MetaData(
         long_name="Height above ground",
@@ -50,16 +49,21 @@ CYCLE_ATTRIBUTES = {
             "specific humidity."
         ),
         positive="up",
+        dimensions=("time", "level"),
     ),
-    "pressure": MetaData(long_name="Pressure", units="Pa"),
-    "temperature": MetaData(long_name="Temperature", units="K"),
+    "pressure": MetaData(long_name="Pressure", units="Pa", dimensions=None),
+    "temperature": MetaData(long_name="Temperature", units="K", dimensions=None),
     "uwind": MetaData(
-        long_name="Zonal wind", units="m s-1", standard_name="eastward_wind"
+        long_name="Zonal wind",
+        units="m s-1",
+        standard_name="eastward_wind",
+        dimensions=None,
     ),
     "vwind": MetaData(
         long_name="Meridional wind",
         units="m s-1",
         standard_name="northward_wind",
+        dimensions=None,
     ),
     "wwind": MetaData(
         long_name="Vertical wind",
@@ -69,13 +73,15 @@ CYCLE_ATTRIBUTES = {
             "The vertical wind has been calculated from omega (Pa s-1),\n"
             "height and pressure using: w=omega*dz/dp"
         ),
+        dimensions=None,
     ),
     "omega": MetaData(
         long_name="Vertical wind in pressure coordinates",
         units="PA s-1",
         standard_name="omega",
+        dimensions=None,
     ),
-    "q": MetaData(long_name="Specific humidity", units="1"),
+    "q": MetaData(long_name="Specific humidity", units="1", dimensions=None),
     "rh": MetaData(
         long_name="Relative humidity",
         units="1",
@@ -83,28 +89,36 @@ CYCLE_ATTRIBUTES = {
             "With respect to liquid above 0 degrees C and with respect to ice\n"
             "below 0 degrees C"
         ),
+        dimensions=None,
     ),
 }
 
 MODEL_L3_ATTRIBUTES = {
-    "cf": MetaData(long_name="Cloud fraction of model grid point", units="1"),
+    "cf": MetaData(
+        long_name="Cloud fraction of model grid point",
+        units="1",
+        dimensions=("time", "level"),
+    ),
     "cf_cirrus": MetaData(
         long_name="Cloud fraction of model grid point with filtered cirrus fraction",
         units="1",
         comment="High level cirrus fraction is reduce do to lack if a radar\n"
         "capability to observe correctly particles small and high.",
+        dimensions=("time", "level"),
     ),
     "iwc": MetaData(
         long_name="Ice water content of model grid point",
         units="kg m-3",
         comment="Calculated using model ice water mixing ration, pressure\n"
         "and temperature: qi*P/287*T",
+        dimensions=("time", "level"),
     ),
     "lwc": MetaData(
         long_name="Liquid water content of model grid point",
         units="kg m-3",
         comment="Calculated using model liquid water mixing ration, pressure\n"
         "and temperature:  ql*P/287*T",
+        dimensions=("time", "level"),
     ),
 }
 
@@ -115,8 +129,9 @@ REGRID_PRODUCT_ATTRIBUTES = {
         comment=(
             "Cloud fraction generated from observations and by volume,\n"
             "averaged onto the models grid with height and time. Volume is\n"
-            "space withing four grid points."
+            "space within four grid points."
         ),
+        dimensions=("time", "level"),
     ),
     "cf_A": MetaData(
         long_name="Observed cloud fraction by area",
@@ -126,6 +141,7 @@ REGRID_PRODUCT_ATTRIBUTES = {
             "averaged onto the models grid with height and time. Area is\n"
             "sum of time columns with any cloud fraction."
         ),
+        dimensions=("time", "level"),
     ),
     "cf_V_adv": MetaData(
         long_name="Observed cloud fraction by advection volume",
@@ -135,6 +151,7 @@ REGRID_PRODUCT_ATTRIBUTES = {
             "except that model winds were used to estimate the time taken to advect\n"
             "airflow a distance equivalent to the models horizontal resolution."
         ),
+        dimensions=("time", "level"),
     ),
     "cf_A_adv": MetaData(
         long_name="Observed cloud fraction by advection area",
@@ -144,6 +161,7 @@ REGRID_PRODUCT_ATTRIBUTES = {
             "except that model winds were used to estimate the time taken to advect\n"
             "airflow a distance equivalent to the models horizontal resolution."
         ),
+        dimensions=("time", "level"),
     ),
     "iwc": MetaData(
         long_name="Observed ice water content reshaped to model dimensions by"
@@ -155,6 +173,7 @@ REGRID_PRODUCT_ATTRIBUTES = {
             "The formula  has been applied where the categorization data has\n"
             "diagnosed that the radar echo is due to ice."
         ),
+        dimensions=("time", "level"),
     ),
     "iwc_att": MetaData(
         long_name="Observed ice water content with attenuation reshaped to model grid"
@@ -165,6 +184,7 @@ REGRID_PRODUCT_ATTRIBUTES = {
             "except that profiles with uncorrected attenuation of the radar\n"
             "reflectivity were included."
         ),
+        dimensions=("time", "level"),
     ),
     "iwc_rain": MetaData(
         long_name="Observed ice water content with raining reshaped to model grid"
@@ -175,6 +195,7 @@ REGRID_PRODUCT_ATTRIBUTES = {
             "including attenuation, iwc_att, except that profiles with rain\n"
             "at the surface were also included."
         ),
+        dimensions=("time", "level"),
     ),
     "iwc_adv": MetaData(
         long_name="Observed ice water content reshaped to model advection grid"
@@ -185,6 +206,7 @@ REGRID_PRODUCT_ATTRIBUTES = {
             "except that model winds were used to estimate the time taken to advect\n"
             "the flow a distance equivalent to the models horizontal resolution."
         ),
+        dimensions=("time", "level"),
     ),
     "iwc_att_adv": MetaData(
         long_name="Observed ice water content with attenuation reshaped to model"
@@ -195,6 +217,7 @@ REGRID_PRODUCT_ATTRIBUTES = {
             "iwc_adv,  except that profiles with uncorrected attenuation of the radar\n"
             "reflectivity were included."
         ),
+        dimensions=("time", "level"),
     ),
     "iwc_rain_adv": MetaData(
         long_name="Observed ice water content with raining reshaped to model"
@@ -205,6 +228,7 @@ REGRID_PRODUCT_ATTRIBUTES = {
             "including attenuation, iwc_att_adv, except that profiles with rain\n"
             "at the surface were also included."
         ),
+        dimensions=("time", "level"),
     ),
     "lwc": MetaData(
         long_name="Observed liquid water content reshaped to model grid by averaging",
@@ -214,6 +238,7 @@ REGRID_PRODUCT_ATTRIBUTES = {
             "pixels where the categorization data has diagnosed that liquid water \n"
             "is present, averaged onto the model grid with height and time."
         ),
+        dimensions=("time", "level"),
     ),
     "lwc_adv": MetaData(
         long_name="Observed liquid water content reshaped to model advection grid by "
@@ -224,5 +249,6 @@ REGRID_PRODUCT_ATTRIBUTES = {
             "lwc, except that model winds were used to estimate the time taken to\n"
             "advect flow over the models grid points."
         ),
+        dimensions=("time", "level"),
     ),
 }

cloudnetpy/model_evaluation/model_metadata.py

@@ -28,7 +28,7 @@ MODELS = {
         level=88,
         cycle="1-12, 7-18",
     ),
-    "harmonie": ModelMetaData(
+    "harmonie-fmi-6-11": ModelMetaData(
         model_name="HARMONIE-AROME",
         long_name="the HIRLAM–ALADIN Research on Mesoscale Operational NWP in Euromed",
         level=65,