From 40c357c6eb6c21629d34340029059e315cabb9cd Mon Sep 17 00:00:00 2001
From: bczernecki <bczernecki@gmail.com>
Date: Tue, 15 Sep 2026 14:30:15 +0200
Subject: [PATCH] fix: broken col_names and attributes for IMGW data

---
 .Rbuildignore                               |   1 +
 .github/copilot-instructions.md             |   2 +-
 .gitignore                                  |   1 +
 CLAUDE.md                                   |   2 +-
 DESCRIPTION                                 |   2 +-
 NEWS.md                                     |   8 +++
 R/hydro_imgw.R                              |  15 +++--
 R/hydro_imgw_daily.R                        |   2 +
 R/hydro_imgw_monthly.R                      |   2 +
 R/hydro_shortening_imgw.R                   |  60 ++++++++++----------
 R/meteo_imgw.R                              |  13 ++---
 R/meteo_imgw_daily.R                        |  12 +---
 R/meteo_imgw_hourly.R                       |  15 ++---
 R/meteo_imgw_monthly.R                      |  15 +----
 R/meteo_shortening_imgw.R                   |  52 +++++++++--------
 R/utils.R                                   |   7 ++-
 README.md                                   |   9 ++-
 data-raw/hydro_parametry_skroty.csv         |   1 +
 data-raw/parametry_skrot.csv                |  32 +++++++++++
 data/imgw_hydro_abbrev.rda                  | Bin 1010 -> 1123 bytes
 data/imgw_meteo_abbrev.rda                  | Bin 5237 -> 5230 bytes
 man/hydro_imgw.Rd                           |   1 +
 man/hydro_shortening_imgw.Rd                |  17 ++----
 man/meteo_imgw.Rd                           |   5 +-
 man/meteo_imgw_daily.Rd                     |   6 --
 man/meteo_imgw_hourly.Rd                    |   5 --
 man/meteo_imgw_monthly.Rd                   |  11 ----
 man/meteo_shortening_imgw.Rd                |  12 ++--
 tests/testthat/test-hydro_imgw.R            |   5 +-
 tests/testthat/test-hydro_shortening_imgw.R |  36 ++++++++++--
 tests/testthat/test-imgw_meteo_abbrev.R     |   2 +-
 tests/testthat/test-meteo_imgw.R            |   4 +-
 tests/testthat/test-meteo_imgw_daily.R      |  40 +++++++++++++
 vignettes/articles/pl.Rmd                   |  20 +++----
 vignettes/getstarted.Rmd                    |   6 ++
 35 files changed, 250 insertions(+), 171 deletions(-)

diff --git a/.Rbuildignore b/.Rbuildignore
index 7c282feb..eb2bd1ed 100644
--- a/.Rbuildignore
+++ b/.Rbuildignore
@@ -21,3 +21,4 @@ vignettes/articles/usecase.Rmd
 ^\.positai$
 ^\.claude$
 CLAUDE.md$
+^\.posit/assistant$
diff --git a/.github/copilot-instructions.md b/.github/copilot-instructions.md
index b21f957b..a3ff330a 100644
--- a/.github/copilot-instructions.md
+++ b/.github/copilot-instructions.md
@@ -28,7 +28,7 @@ Run commands from the package root.
   - **OGIMET**: HTML is scraped with `XML::readHTMLTable`; station identity is based on WMO IDs. Hourly precipitation post-processing is handled by `precip_split()`.
   - **NOAA / Wyoming**: direct file or page downloads for ISH hourly data, Mauna Loa CO2, and Wyoming soundings.
 
-- IMGW column renaming is a distinct normalization layer. Most IMGW functions accept `col_names = "short" | "full" | "polish"` and pass results through `meteo_shortening_imgw()` or `hydro_shortening_imgw()`. The mapping tables live in built-in datasets backed by `data-raw/`.
+- IMGW column renaming is a distinct normalization layer. IMGW functions return short English names and attach the original full parameter label as a per-column `label` attribute through `meteo_shortening_imgw()` or `hydro_shortening_imgw()`. The mapping tables live in built-in datasets backed by `data-raw/`.
 
 - Package data and docs follow standard R package patterns:
   - exported code in `R/`
diff --git a/.gitignore b/.gitignore
index 17050e89..287ea2d0 100644
--- a/.gitignore
+++ b/.gitignore
@@ -16,3 +16,4 @@ pkgdown
 test-out.txt
 .positai
 .aider*
+.posit/assistant
diff --git a/CLAUDE.md b/CLAUDE.md
index 7379afe7..66180ce8 100644
--- a/CLAUDE.md
+++ b/CLAUDE.md
@@ -29,7 +29,7 @@ Run from the package root in R:
 - **OGIMET**: HTML scraping via `XML::readHTMLTable` from `ogimet.com`. Stations are identified by WMO ID. `precip_split` / `R/precip_split.R` handles 6/12/24h precipitation disaggregation for hourly data.
 - **NOAA / Wyoming**: direct file downloads (ISH gzipped fixed-width, CO2 text, sounding HTML).
 
-**Column-name shortening layer.** Most IMGW download functions accept `col_names = "short" | "full" | "polish"` and pass the raw frame through `meteo_shortening_imgw()` / `hydro_shortening_imgw()` (in `R/*_shortening_imgw.R`). Full and short names are looked up against `imgw_meteo_abbrev` / `imgw_hydro_abbrev` (built-in data). When you add a new IMGW column, update both the abbrev table (`data-raw/`) and the shortener.
+**Column-name shortening layer.** IMGW download functions return short English names and attach the original full IMGW parameter label as a per-column `label` attribute. Names are looked up against `imgw_meteo_abbrev` / `imgw_hydro_abbrev` (built-in data). When you add a new IMGW column, update both the abbrev table (`data-raw/`) and the shortener.
 
 **Graceful network failure** is required for CRAN. Use `test_url()` (`R/test_url.R`) to gate downloads, and follow the existing `allow_failure = TRUE` pattern: wrap the real worker (`*_bp` "best practice" inner function) in `tryCatch` so user-facing functions return `NULL`/`invisible()` with a `message()` instead of erroring. Tests follow the same convention — every network test starts with `if (!curl::has_internet()) return(invisible(NULL))`. Don't add tests that fail when offline.
 
diff --git a/DESCRIPTION b/DESCRIPTION
index 023fb520..4a5974d0 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: climate
 Title: Interface to Download Meteorological (and Hydrological) Datasets
-Version: 1.4.0
+Version: 1.4.1
 Authors@R: c(person(given = "Bartosz",
            family = "Czernecki",
            role = c("aut", "cre"),
diff --git a/NEWS.md b/NEWS.md
index 3fd7bc66..e78a2215 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -1,3 +1,11 @@
+# climate 1.4.1
+
+* fixed IMGW meteorological and hydrological parameter mappings for current metadata labels
+* removed the `col_names` choice between short, full, and Polish column names to simplify usage
+* original full IMGW parameter names are available through each column's `label` attribute
+* updated documentation, README, and vignettes to describe the simplified naming scheme
+
+
 # climate 1.4.0
 
 * adding the `synop_parser()` function for reading raw SYNOP messages
diff --git a/R/hydro_imgw.R b/R/hydro_imgw.R
index 551b29ce..3225f903 100644
--- a/R/hydro_imgw.R
+++ b/R/hydro_imgw.R
@@ -16,6 +16,7 @@
 #' @return A data.frame with columns describing the hydrological parameters
 #' (e.g. flow, water level) where each row represent a measurement,
 #' depending on the interval, at a given hour, month or year.
+#' IMGW parameter columns carry a `label` attribute with the original metadata label.
 #' If `coords = TRUE` additional two columns with geographic coordinates are added.
 #' @examples
 #' \donttest{
@@ -30,12 +31,18 @@ hydro_imgw = function(interval,
 
   if (interval == "daily") {
     # dobowe
-    calosc = hydro_imgw_daily(year = year, station = station, ...)
+    calosc = hydro_imgw_daily(
+      year = year,
+      station = station,
+      ...
+    )
   } else if (interval == "monthly") {
     # miesieczne
-    calosc = hydro_imgw_monthly(year = year,
-                                station = station,
-                                ...)
+    calosc = hydro_imgw_monthly(
+      year = year,
+      station = station,
+      ...
+    )
   } else{
     stop("Wrong `interval` value. It should be either 'daily' or 'monthly'", call. = FALSE)
   }
diff --git a/R/hydro_imgw_daily.R b/R/hydro_imgw_daily.R
index 0a8a1bf2..c68c301e 100644
--- a/R/hydro_imgw_daily.R
+++ b/R/hydro_imgw_daily.R
@@ -174,6 +174,8 @@ hydro_imgw_daily_bp = function(year,
       attr(all_data[[cols]], "label") = meta$label[ind]
     }
   }
+  all_data = imgw_rename_params_to_labels(all_data, meta)
+  all_data = hydro_shortening_imgw(all_data, ...)
   
   return(all_data)
 }
\ No newline at end of file
diff --git a/R/hydro_imgw_monthly.R b/R/hydro_imgw_monthly.R
index c3009e2a..39da0803 100644
--- a/R/hydro_imgw_monthly.R
+++ b/R/hydro_imgw_monthly.R
@@ -127,6 +127,8 @@ hydro_imgw_monthly_bp = function(year,
       attr(all_data[[cols]], "label") = meta$label[ind]
       }
   }
+  all_data = imgw_rename_params_to_labels(all_data, meta)
+  all_data = hydro_shortening_imgw(all_data, ...)
 
   return(all_data)
 }
diff --git a/R/hydro_shortening_imgw.R b/R/hydro_shortening_imgw.R
index 7f58ebbc..80ab175e 100644
--- a/R/hydro_shortening_imgw.R
+++ b/R/hydro_shortening_imgw.R
@@ -3,58 +3,58 @@
 #' Shortening column names of hydrological parameters to improve the readability of downloaded dataset from
 #' the danepubliczne.imgw.pl collection and removing duplicated column names
 #' @param data downloaded dataset with original column names
-#' @param col_names three types of column names possible: "short" - default,
-#' values with shorten names,
-#' "full" - full English description,
-#' "polish" - original names in the dataset
 #' @param remove_duplicates whether to remove duplicated column names
 #' (default TRUE - i.e., columns with duplicated names are deleted)
 #' @export
-#' @returns data.frame with shorten names of hydrological parameters
+#' @returns data.frame with short English names of hydrological parameters.
+#' Existing column attributes, including the original IMGW `label` metadata,
+#' are preserved.
 #' @examples
 #' \donttest{
-#'   monthly = data = hydro_imgw("monthly", year = 1969, col_names = "polish")
+#'   monthly = data = hydro_imgw("monthly", year = 1969)
 #'   
 #'   if (is.data.frame(monthly)) {
-#'   abbr = hydro_shortening_imgw(data = monthly,
-#'       col_names = "full",
-#'       remove_duplicates = TRUE)
+#'   abbr = hydro_shortening_imgw(data = monthly, remove_duplicates = TRUE)
 #'   head(abbr)
 #'   }
 #' }
 #'
 
 hydro_shortening_imgw = function(data,
-                                 col_names = "short",
                                  remove_duplicates = TRUE) {
 
-  if (col_names != "polish") {
-    abbrev = climate::imgw_hydro_abbrev
-    # additional workarounds for mac os but not only...
-    abbrev$fullname = gsub(x = abbrev$fullname, pattern = "'", replacement = "")
-    abbrev$fullname = gsub(x = abbrev$fullname, pattern = "\\^", replacement = "")
-    # end of workaround
-    orig_columns = trimws(gsub("\\s+", " ", colnames(data))) # remove double spaces
+  data = as.data.frame(data)
+  column_attributes = lapply(data, attributes)
 
-    matches = match(orig_columns, abbrev$fullname)
-    matches = matches[!is.na(matches)]
+  abbrev = climate::imgw_hydro_abbrev
+  # additional workarounds for mac os but not only...
+  abbrev$fullname = gsub(x = abbrev$fullname, pattern = "'", replacement = "")
+  abbrev$fullname = gsub(x = abbrev$fullname, pattern = "\\^", replacement = "")
+  abbrev$fullname = stringi::stri_trans_general(abbrev$fullname, "LATIN-ASCII")
+  # end of workaround
+  orig_columns = trimws(gsub("\\s+", " ", colnames(data))) # remove double spaces
+  orig_columns = gsub(x = orig_columns, pattern = "'", replacement = "")
+  orig_columns = gsub(x = orig_columns, pattern = "\\^", replacement = "")
+  orig_columns = stringi::stri_trans_general(orig_columns, "LATIN-ASCII")
+  # `Data` is created by the package and is not an IMGW parameter.
+  orig_columns[orig_columns == "Data"] = NA_character_
 
-    if (col_names == "short") {
-      # abbrev english
-      colnames(data)[orig_columns %in% abbrev$fullname] = abbrev$abbr_eng[matches]
-    }
-
-    if (col_names == "full") {
-      # full english names:
-      colnames(data)[orig_columns %in% abbrev$fullname] = abbrev$fullname_eng[matches]
-    }
-  }
+  matches = match(orig_columns, abbrev$fullname)
+  matches = matches[!is.na(matches)]
+  colnames(data)[orig_columns %in% abbrev$fullname] = abbrev$abbr_eng[matches]
 
   # removing duplicated column names:  (e.g. station's name)
   if (remove_duplicates == TRUE) {
-    data = data[, !duplicated(colnames(data))]
+    keep = !duplicated(colnames(data))
+    data = data[, keep, drop = FALSE]
+    column_attributes = column_attributes[keep]
   }
   data = unique(data)
+  for (ind in seq_along(column_attributes)) {
+    if (!is.null(column_attributes[[ind]])) {
+      attributes(data[[ind]]) = column_attributes[[ind]]
+    }
+  }
   rownames(data) = NULL
   return(data)
 }
diff --git a/R/meteo_imgw.R b/R/meteo_imgw.R
index f53fe598..9b801e4d 100644
--- a/R/meteo_imgw.R
+++ b/R/meteo_imgw.R
@@ -15,8 +15,6 @@
 #'   (default `FALSE` — status columns are deleted). Not used when `rank = "telemetry"`.
 #' @param coords add coordinates of the station (logical value `TRUE` or `FALSE`).
 #'   Default `FALSE`.
-#' @param col_names column name style: `"short"` (default), `"full"` (English descriptions),
-#'   or `"polish"` (original dataset names). Not used when `rank = "telemetry"`.
 #' @param station name of meteorological station(s).
 #'   For ranks `"synop"`, `"climate"`, `"precip"`: station name(s) in CAPITAL LETTERS.
 #'   Please note that station names may change over time — sometimes two names are required,
@@ -33,6 +31,7 @@
 #' @return A data.frame with meteorological parameters where each row is a measurement.
 #'   For ranks `"synop"`, `"climate"`, `"precip"`: measurements at a given hour, day, or month,
 #'   depending on `interval`. If `coords = TRUE` two additional coordinate columns are appended.
+#'   IMGW parameter columns also carry a `label` attribute with the original Polish metadata label.
 #'   For `rank = "telemetry"`: a data.table with 10-minute interval observations (not
 #'   expert-validated). If `coords = TRUE` columns `name`, `lon`, `lat`, and `alt` are appended.
 #' @examples
@@ -52,7 +51,6 @@ meteo_imgw = function(interval = NULL,
                       status = FALSE,
                       coords = FALSE,
                       station = NULL,
-                      col_names = "short",
                       parameters = NULL,
                       ...) {
   if (rank == "telemetry") {
@@ -70,22 +68,19 @@ meteo_imgw = function(interval = NULL,
                               year      = year,
                               status    = status,
                               coords    = coords,
-                              station   = station,
-                              col_names = col_names, ...)
+                              station   = station, ...)
   } else if (interval == "monthly") {
     result = meteo_imgw_monthly(rank      = rank,
                                 year      = year,
                                 status    = status,
                                 coords    = coords,
-                                station   = station,
-                                col_names = col_names, ...)
+                                station   = station, ...)
   } else if (interval == "hourly") {
     result = meteo_imgw_hourly(rank      = rank,
                                year      = year,
                                status    = status,
                                coords    = coords,
-                               station   = station,
-                               col_names = col_names, ...)
+                               station   = station, ...)
   } else {
     stop("Wrong `interval` value. It should be either 'hourly', 'daily', or 'monthly'.")
   }
diff --git a/R/meteo_imgw_daily.R b/R/meteo_imgw_daily.R
index 0a41cc79..e22b089d 100644
--- a/R/meteo_imgw_daily.R
+++ b/R/meteo_imgw_daily.R
@@ -14,10 +14,6 @@
 #' database and thus providing both names is needed
 #' (e.g. `station = c("POZNAŃ", "POZNAŃ-ŁAWICA", "WARSZAWA", "WARSZAWA-OKĘCIE")`).
 #' Stations' IDs (numeric) are no longer valid
-#' @param col_names three types of column names possible:
-#' "short" - default, values with shorten names,
-#' "full" - full English description,
-#' "polish" - original names in the dataset
 #' @param allow_failure logical - whether to proceed or stop on failure. By default set to TRUE (i.e. don't stop on error). For debugging purposes change to FALSE
 #' @param ... other parameters that may be passed to the 'shortening' function that
 #' shortens column names
@@ -35,7 +31,6 @@ meteo_imgw_daily = function(rank = "synop",
                              status = FALSE,
                              coords = FALSE,
                              station = NULL,
-                             col_names = "short",
                              allow_failure = TRUE,
                              ...) {
   if (allow_failure) {
@@ -45,8 +40,7 @@ meteo_imgw_daily = function(rank = "synop",
         year,
         status,
         coords,
-        station,
-        col_names
+        station
       ),
       error = function(e) {
         message(paste(
@@ -63,7 +57,6 @@ meteo_imgw_daily = function(rank = "synop",
       status,
       coords,
       station,
-      col_names,
       ...
     )
   }
@@ -74,7 +67,6 @@ meteo_imgw_daily_bp = function(rank,
                                 status,
                                 coords,
                                 station,
-                                col_names,
                                 ...) {
   translit = check_locale()
   base_url = "https://danepubliczne.imgw.pl/data/dane_pomiarowo_obserwacyjne/"
@@ -399,7 +391,7 @@ meteo_imgw_daily_bp = function(rank,
   }
 
   all_data = imgw_rename_params_to_labels(all_data, meta)
-  all_data = meteo_shortening_imgw(all_data, col_names = col_names, remove_duplicates = TRUE)
+  all_data = meteo_shortening_imgw(all_data, remove_duplicates = TRUE)
 
   # check if there any messages gathered in env$logs and if it is not empty then print them:
   if (length(env$logs) > 0) {
diff --git a/R/meteo_imgw_hourly.R b/R/meteo_imgw_hourly.R
index 467add7e..b34e7456 100644
--- a/R/meteo_imgw_hourly.R
+++ b/R/meteo_imgw_hourly.R
@@ -9,9 +9,6 @@
 #' (default status = FALSE - i.e. the status columns are deleted)
 #' @param coords add coordinates of the station (logical value TRUE or FALSE)
 #' @param station name of meteorological station(s) (character vector)
-#' @param col_names three types of column names possible: "short" - default,
-#' values with shorten names, "full" - full English description,
-#' "polish" - original names in the dataset
 #' @param allow_failure logical - whether to proceed or stop on failure. By default set to TRUE (i.e. don't stop on error). For debugging purposes change to FALSE
 #' @param ... other parameters that may be passed to the 'shortening'
 #' function that shortens column names
@@ -31,7 +28,6 @@ meteo_imgw_hourly = function(rank = "synop",
                               status = FALSE,
                               coords = FALSE,
                               station = NULL,
-                              col_names = "short",
                               allow_failure = TRUE,
                               ...) {
   if (allow_failure) {
@@ -41,8 +37,7 @@ meteo_imgw_hourly = function(rank = "synop",
         year,
         status,
         coords,
-        station,
-        col_names, ...
+        station, ...
       ),
       error = function(e) {
         message(paste(
@@ -58,8 +53,7 @@ meteo_imgw_hourly = function(rank = "synop",
       year,
       status,
       coords,
-      station,
-      col_names, ...
+      station, ...
     )
   }
 }
@@ -70,8 +64,7 @@ meteo_imgw_hourly_bp = function(rank,
                                  year,
                                  status,
                                  coords,
-                                 station,
-                                 col_names, ...) {
+                                 station, ...) {
   translit = check_locale()
   stopifnot(rank == "synop" | rank == "climate") # for hourly data only synop and climate has data
   base_url = "https://danepubliczne.imgw.pl/data/dane_pomiarowo_obserwacyjne/"
@@ -279,7 +272,7 @@ meteo_imgw_hourly_bp = function(rank,
   }
 
   all_data = imgw_rename_params_to_labels(all_data, meta)
-  all_data = meteo_shortening_imgw(all_data, col_names = col_names, ...)
+  all_data = meteo_shortening_imgw(all_data, ...)
 
   # check if there any messages gathered in env$logs and if it is not empty then print them:
   if (length(env$logs) > 0) {
diff --git a/R/meteo_imgw_monthly.R b/R/meteo_imgw_monthly.R
index f8c3d9db..e40e9394 100644
--- a/R/meteo_imgw_monthly.R
+++ b/R/meteo_imgw_monthly.R
@@ -12,9 +12,6 @@
 #' It accepts names (characters in CAPITAL LETTERS). Stations' IDs (numeric) are no longer supported.
 #' Please note that station names may change over time and thus sometimes 2 names
 #' are required in some cases, e.g. `c("POZNAŃ", "POZNAŃ-ŁAWICA")`.
-#' @param col_names three types of column names possible: "short" - default,
-#' values with shorten names, "full" - full English description,
-#' "polish" - original names in the dataset
 #' @param allow_failure logical - whether to proceed or stop on failure.
 #' By default set to TRUE (i.e. don't stop on error). For debugging purposes change to FALSE
 #' @param ... other parameters that may be passed to the
@@ -29,12 +26,6 @@
 #' monthly = meteo_imgw_monthly(rank = "climate", year = 1969)
 #' head(monthly)
 #'
-#' # a descriptive (long) column names:
-#' monthly2 = meteo_imgw_monthly(
-#'   rank = "synop", year = 2018,
-#'   col_names = "full"
-#' )
-#' head(monthly2)
 #' }
 #'
 meteo_imgw_monthly = function(rank = "synop",
@@ -42,7 +33,6 @@ meteo_imgw_monthly = function(rank = "synop",
                               status = FALSE,
                               coords = FALSE,
                               station = NULL,
-                              col_names = "short",
                               allow_failure = TRUE,
                               ...) {
   if (allow_failure) {
@@ -53,7 +43,6 @@ meteo_imgw_monthly = function(rank = "synop",
         status,
         coords,
         station,
-        col_names,
         ...
       ),
       error = function(e) {
@@ -71,7 +60,6 @@ meteo_imgw_monthly = function(rank = "synop",
       status,
       coords,
       station,
-      col_names,
       ...
     )
   }
@@ -84,7 +72,6 @@ meteo_imgw_monthly_bp = function(rank,
                                   status,
                                   coords,
                                   station,
-                                  col_names,
                                   ...) {
   translit = check_locale()
   base_url = "https://danepubliczne.imgw.pl/data/dane_pomiarowo_obserwacyjne/"
@@ -245,6 +232,6 @@ meteo_imgw_monthly_bp = function(rank,
   }
 
   all_data = imgw_rename_params_to_labels(all_data, meta)
-  all_data = meteo_shortening_imgw(all_data, col_names = col_names, ...)
+  all_data = meteo_shortening_imgw(all_data, ...)
   return(all_data)
 }
diff --git a/R/meteo_shortening_imgw.R b/R/meteo_shortening_imgw.R
index 4c87d894..8f44fa39 100644
--- a/R/meteo_shortening_imgw.R
+++ b/R/meteo_shortening_imgw.R
@@ -3,30 +3,36 @@
 #' Shortening column names of meteorological parameters to improve the readability of downloaded dataset from the danepubliczne.imgw.pl collection and removing duplicated column names 
 #'
 #' @param data downloaded dataset with original column names
-#' @param col_names three types of column names possible: "short" - default, values with shorten names, "full" - full English description, "polish" - original names in the dataset
 #' @param remove_duplicates whether to remove duplicated column names (default TRUE - i.e., columns with duplicated names are deleted)
 #' @export
-#' @returns data.frame with modified names of meteorological parameters
+#' @returns data.frame with short English names of meteorological parameters.
+#' Existing column attributes, including the original IMGW `label` metadata,
+#' are preserved.
 #' 
 #' @examples 
 #' \donttest{
 #'   monthly = meteo_imgw("monthly", rank = "climate", year = 1969)
 #'   
-#'   abbr = meteo_shortening_imgw(data = monthly,
-#'       col_names = "full", 
-#'       remove_duplicates = TRUE)
+#'   abbr = meteo_shortening_imgw(data = monthly, remove_duplicates = TRUE)
 #'   head(abbr)
 #' }
 #'
 
-meteo_shortening_imgw = function(data, col_names = "short", remove_duplicates = TRUE) {
+meteo_shortening_imgw = function(data, remove_duplicates = TRUE) {
 
   data = as.data.frame(data)
+  column_attributes = lapply(data, attributes)
+
   # removing duplicated column names:  (e.g. station's name)
   if (remove_duplicates == TRUE) {
-    data = data[, !duplicated(colnames(data))]
+    keep = !duplicated(colnames(data))
+    data = data[, keep, drop = FALSE]
+    column_attributes = column_attributes[keep]
     # fix for merged station names with suffixes
     if (any(colnames(data) %in% c("Nazwa stacji.x", "Nazwa stacji.y"))) {
+      keep = colnames(data) != "Nazwa stacji.y"
+      data = data[, keep, drop = FALSE]
+      column_attributes = column_attributes[keep]
       data$`Nazwa stacji.y` = NULL
       colnames(data)[colnames(data) == "Nazwa stacji.x"] = "Nazwa stacji"
     }
@@ -34,30 +40,26 @@ meteo_shortening_imgw = function(data, col_names = "short", remove_duplicates =
     # fix for mean air temperature which is stated sometimes in two files as:
     # "Srednia dobowa temperatura[°C]" and "Srednia temperatura dobowa [°C]"
     if (any(grepl(x = colnames(data), "Srednia dobowa temperatura"))) {
-      data[, which(grepl(x = colnames(data), "Srednia dobowa temperatura"))] = NULL
+      keep = !grepl(x = colnames(data), "Srednia dobowa temperatura")
+      data = data[, keep, drop = FALSE]
+      column_attributes = column_attributes[keep]
     }
   }
   
-  if (col_names != "polish") {
-    abbrev = climate::imgw_meteo_abbrev
-    orig_columns = trimws(gsub("\\s+", " ", colnames(data))) # remove double spaces
-    orig_columns = trimws(gsub("\\[.*?]", "", orig_columns)) # remove brackets and content inside
-
-    abbrev$fullname = trimws(gsub("\\[.*?]", "", abbrev$fullname))
-    matches = match(orig_columns, abbrev$fullname)
-    matches = matches[!is.na(matches)]
-
-    if (col_names == "short") {
-      # abbrev english
-      colnames(data)[orig_columns %in% abbrev$fullname] = abbrev$abbr_eng[matches]
-    }
+  abbrev = climate::imgw_meteo_abbrev
+  orig_columns = trimws(gsub("\\s+", " ", colnames(data))) # remove double spaces
+  orig_columns = trimws(gsub("\\[.*?]", "", orig_columns)) # remove brackets and content inside
 
-    if (col_names == "full") {
-      # full english names:
-      colnames(data)[orig_columns %in% abbrev$fullname] = abbrev$fullname_eng[matches]
+  abbrev$fullname = trimws(gsub("\\[.*?]", "", abbrev$fullname))
+  matches = match(orig_columns, abbrev$fullname)
+  matches = matches[!is.na(matches)]
+  colnames(data)[orig_columns %in% abbrev$fullname] = abbrev$abbr_eng[matches]
+  data = unique(data)
+  for (ind in seq_along(column_attributes)) {
+    if (!is.null(column_attributes[[ind]])) {
+      attributes(data[[ind]]) = column_attributes[[ind]]
     }
   }
-  data = unique(data)
   rownames(data) = NULL
   return(data)
 }
diff --git a/R/utils.R b/R/utils.R
index 997b5396..6c19161c 100644
--- a/R/utils.R
+++ b/R/utils.R
@@ -1,4 +1,4 @@
-#' Rename parameter code columns to full Polish label names
+#' Rename parameter code columns to full Polish label names and attach metadata labels
 #'
 #' Translates short internal parameter codes (e.g. NSP, POST, TMAX) to the
 #' original full Polish labels stored in the metadata, so that the result can
@@ -18,7 +18,10 @@ imgw_rename_params_to_labels = function(data, meta) {
   param_map = setNames(meta$label, meta$parameters)
   cols = colnames(data)
   in_map = cols %in% names(param_map)
-  colnames(data)[in_map] = param_map[cols[in_map]]
+  for (ind in which(in_map)) {
+    attr(data[[ind]], "label") = unname(param_map[cols[ind]])
+  }
+  colnames(data)[in_map] = unname(param_map[cols[in_map]])
   data
 }
 
diff --git a/README.md b/README.md
index c370eb49..2311c4ba 100644
--- a/README.md
+++ b/README.md
@@ -61,6 +61,8 @@ Downloading measurements of the vertical profile of atmosphere (aka rawinsonde d
 Downloading hourly, daily, and monthly meteorological data from the Polish met service across 
 all types of stations (i.e. SYNOP/CLIMATE/PRECIP ) available in the danepubliczne.imgw.pl collection. 
 It is a wrapper for `meteo_monthly()`, `meteo_daily()`, `meteo_hourly()` and `meteo_imgw_datastore()` which gives access from montly to even to 10-min dataset.
+IMGW archive outputs use short English parameter names; the original IMGW names are
+available as `label` attributes on the corresponding columns.
 
   
 ### (Polish) Hydrological data
@@ -69,6 +71,8 @@ It is a wrapper for `meteo_monthly()`, `meteo_daily()`, `meteo_hourly()` and `me
 Downloading hourly, daily, and monthly hydrological data from stations available in the
 danepubliczne.imgw.pl collection.
 It is a wrapper for previously developed set of functions such as: `hydro_monthly()`, and `hydro_daily()`
+IMGW archive outputs use short English parameter names; the original IMGW names are
+available as `label` attributes on the corresponding columns.
 
 - 🇵🇱 [**hydro_imgw_datastore**](https://bczernecki.github.io/climate/reference/hydro_imgw_datastore.html) - 
 Downloading hourly and subhourly hydrological data from the IMGW-PIB hydro telemetry stations.
@@ -193,7 +197,7 @@ h = hydro_imgw(interval = "daily", year = 2010:2011)
 head(h)
 ```
 
-| id        | station | riv_or_lake | date       | hyy  | idhyy | dd | H   | Q   | T  | mm | thick |
+| id        | station | riv_or_lake | Data       | hyy  | idhyy | dd | H   | Q   | T  | mm | thick |
 |-----------|---------|-------------|------------|------|-------|----|-----|-----|----|----|-------|
 | 150210180 | ANNOPOL | Wisła (2)   | 2009-11-01 | 2010 | 1     | 1  | 287 | 436 | NA | 11 | NA    |
 | 150210180 | ANNOPOL | Wisła (2)   | 2009-11-02 | 2010 | 1     | 2  | 282 | 412 | NA | 11 | NA    |
@@ -202,6 +206,9 @@ head(h)
 | 150210180 | ANNOPOL | Wisła (2)   | 2009-11-05 | 2010 | 1     | 5  | 264 | 336 | NA | 11 | NA    |
 | 150210180 | ANNOPOL | Wisła (2)   | 2009-11-06 | 2010 | 1     | 6  | 260 | 320 | NA | 11 | NA    |
 
+The original IMGW description of `H` can be inspected with
+`attr(h$H, "label")`.
+
 ## Example 5
 #### Create Walter & Lieth climatic diagram based on downloaded data
 
diff --git a/data-raw/hydro_parametry_skroty.csv b/data-raw/hydro_parametry_skroty.csv
index 25582403..da5bcc39 100644
--- a/data-raw/hydro_parametry_skroty.csv
+++ b/data-raw/hydro_parametry_skroty.csv
@@ -27,3 +27,4 @@ Grubość lodu [cm];thick;Thickness of ice [cm];;
 Kod zjawiska lodowego (słownik poniżej);id_ice;Ice phenomena;;
 "Procent udziału zjawiska lodowego [mnożnik *10; np. 3 oznacza 30% udziału zjawisk lodowych]";p_ice;Percentage of ice phenomena;;
 Kod zarastania (słownik poniżej);id_over_grow;Overgrowing;;
+Kod zarastania (szczegółowe informacje poniżej);id_over_grow;Overgrowing;;
diff --git a/data-raw/parametry_skrot.csv b/data-raw/parametry_skrot.csv
index 425c383b..9af97863 100644
--- a/data-raw/parametry_skrot.csv
+++ b/data-raw/parametry_skrot.csv
@@ -253,3 +253,35 @@ Nazwa stacji.x,station,Station
 Nazwa stacji.y,station,Station
 rank_code,rank,Station rank code
 Zapas wody w sniegu [mm],water_in_snow,Water content in snow [mm]
+Absolutna maksymalna temperatura powietrza [C],tmax_abs,Absolute maximum air temperature [C]
+Absolutna minimalna temperatura powietrza [C],tmin_abs,Absolute minimum air temperature [C]
+Charakterystyka tendencji cisnienia [kod],press_tend,Pressure tendency
+Cisnienie atmosferyczne na poziomie morza [hPa],slp,Sea level pressure [hPa]
+Cisnienie atmosferyczne na poziomie stacji [hPa],press,Station level pressure [hPa]
+Liczba dni ze zmetnieniem opalizujacym,hazyness_days,Days with hazyness
+Maksymalna dobowa suma opadow w miesiacu [mm],rr_max_daily,Maximum daily precipitation [mm]
+Maksymalna dobowa temperatura powietrza [C],tmax_daily,Maximum daily air temperature [C]
+Minimalna dobowa temperatura powietrza [C],tmin_daily,Minimum daily air temperature [C]
+Minimalna dobowa temperatura powietrza przy gruncie [C],t5cm_min,Minimum daily near surface air temperature [C]
+Minimalna temperatura powietrza przy gruncie za 12 godzin [C],tmin_soil_12h,Minimum near surface air temperature for 12 hours [C]
+Minimalna temperatura powietrza za 12 godzin [C],tmin_12h,Minimum air temperature for 12 hours [C]
+Minimalna temperatura powietrza przy gruncie [C],tmin_soil,Minimum near surface air temperature [C]
+Temperatura minimalna powietrza przy gruncie za 12 godzin [C],tmin_soil_12h,Minimum near surface air temperature for 12 hours [C]
+Niedosyt wilgotnosci [hPa],vapor_press_deficit,Vapour pressure deficit [hPa]
+Srednia dobowa temperatura powietrza [C],t2m_mean_daily,Daily mean air temperature [C]
+Srednie dobowe cisnienie na poziomie morza [hPa],slp_mean_daily,Daily mean sea level pressure [hPa]
+Srednie dobowe cisnienie na poziomie stacji [hPa],press_mean_daily,Daily mean pressure at station level [hPa]
+Srednie dobowe cisnienie pary wodnej [hPa],vapor_press_mean_daily,Daily mean vapour pressure [hPa]
+Srednia maksymalna temperatura powietrza [C],tmax_mean,Monthly mean maximum air temperature [C]
+Srednia minimalna temperatura powietrza [C],tmin_mean,Monthly mean minimum air temperature [C]
+Srednia miesieczna temperatura powietrza [C],t2m_mean_mon,Monthly mean air temperature [C]
+Srednie miesieczne cisnienie na poziomie morza [hPa],slp_mean_mon,Monthly mean sea level pressure [hPa]
+Temperatura gruntu na glebokosci 5 cm [C],t5_soil,Soil temperature 5 cm below surface [C]
+Temperatura gruntu na glebokosci 10 cm [C],t10_soil,Soil temperature 10 cm below surface [C]
+Temperatura gruntu na glebokosci 20 cm [C],t20_soil,Soil temperature 20 cm below surface [C]
+Temperatura gruntu na glebokosci 50 cm [C],t50_soil,Soil temperature 50 cm below surface [C]
+Temperatura gruntu na glebokosci 100 cm [C],t100_soil,Soil temperature 100 cm below surface [C]
+Temperatura maksymalna powietrza za 12 godzin [C],tmax_12h,Maximum air temperature for 12 hours [C]
+Temperatura minimalna powietrza za 12 godzin [C],tmin_12h,Minimum air temperature for 12 hours [C]
+Wartosc tendencji cisnienia [wartosc],tendency,Pressure tendency [value]
+Zachmurzenie ogolne [0-10 do dn. 31.12.1988/oktanty od dn. 01.01.1989],cl_tot,"Total cloudiness [0-10 to 1988, octants since 1999]"
diff --git a/data/imgw_hydro_abbrev.rda b/data/imgw_hydro_abbrev.rda
index 75a94a80c3bcbced76cce22a5a232326431ef8f1..1690126b84f1d064f3911b06e7ee597d6f316716 100644
GIT binary patch
literal 1123
zcmV-p1f2UqT4*^jL0KkKS?kY*ivR*vf9?PO|I1SM{#p<3{_DU0-|$v|00;;a4ge4U
z00Gbh9sq;@1EAE<nx~T}^n=nipaGyBrhpn`(<9SRWYf}mhJYFk00V#k0iXZ_BSRnn
z0000P82|ue5lmFm00009fB*p)3537^00dx61AqVlpa26SLm&VE000>o003kOO(2Pr
zKujm7#MH`WhM|xE2AQUlC#W<UWEf<RQIYB-px6~aUdJ0E3{5cKmw)j3d0(T(M03D9
z(!+j<gg!4ES%$pG8~x|=;9xw-yhn%PfXo3}4qTixQQ~y+<j>msB23H+zsnm}LJGKT
zx$P`ZC~`h5HCUdFs&x45b4_piktGS)7`!81vCNBPFs5d1lK>>Z$uN*4B*;iWkc5(8
z7(khl24)ffbj%Fi6mX_tpnz<5Ek~O+V9sC4xUTl}cQ+-b6Yee(g@i8N7!J0J2{aP6
z8_3gyiIN<qe8mk?d!AAn?%yZP3P?XnzCN{7Y&l4&b&9g6ND+<|@hS3ZabS$zmlh{a
zuZ<nSujRVgq)N;A)ezb!$E#Vcy1|b6^W?7QsOXKvxT<>yF`COdfq0RwPy)hDsFH;-
zT*JBG@S#V@L?yEVM&`VR_$k1EAc9m=+EAF{1zuqy-lrMNP^dT}#;thRrsL7Vt6jpw
z77LdCJC{}6b*fZZFgHa>Z<ax6%{kPvGN(`?v^mI$@MHloGf<EU6ozLR^@_~!TcMFj
zLA5Tp=Mt7kNk;9mwc4l~0Z1|$k;45mC2U0?l!J-9WQm?b0dziw_dxFd4b0OQAP(4p
zcAPX|<NmI-2L}*n0-0>g#tOM1Gdn*|7WvpV$+FZGwgwXsSwv~lVA&u#A(6!hAlKn@
zQVN2B0xJ$>xumY6G!IKY!BpBuR`{Ss=-^Gab^<_f6Tntg4fg4W^5hCbMl0G(CRw)0
zik}1Z4#^R^6AeaE!u&C|O@Z^oqQkw95pz}c&T1|f3rHIK*$va8ssmkn0BCHmX)`RR
z$=8gD7%mW0l#+vN{Z|6xZQ9f}Zd7O%<2^18<@6Mk&w3CiY^apC!hncSAD?5VJtkIp
zqk`xw`xV}37l$BHLgYzGfk@L-Vnk^Qk;(!CnD7TjM1i%<(h{KEa5|or+Z67!j5nR5
z7DeND(^zabFl`1fyc)TS#6k^2HsN8q-3&!FlH78Vp?WWlIf@kL;JjO4NVEjTaoEn|
zy#xnKL4rhBf!Em`7yzb6cIE!ID+mRfXPp)Sd^V`LY;5TzzixPj(|Q7kZ&uJZz;K5Z
zB;!B5wA#sD#7J^JhX6@HC?Xo1gLz4g3uf)5bQXr-QK%PG3k*4Qi{KnH8xYrg0Lu0d
zWJXZHfC3NBoiY%FA<!-#;BaCnLO$?-qyTWdf-f;RfoU1k>mdpZu(hIgLQYNF9Ld;w
peC)?N3Aw}ty<Gv<=fWz97uxmp5T&2~?+9P`yOJrwgoj>yELc5g_(T8z

literal 1010
zcmV<O0}cEiiwFP!000002Ax)4ZyQArztnc2Nkjj*DM)}cLP8{!lQuPxP^&)Gq=~>Z
ze~c)y+SBH4ZLjZk*V?<w#rK3F5Fer+gE!PSdf$Z4L7267<I6b|sSmbi_BX$o@yzVp
z9kgCF7hfzoj#F`})di<gL(!@I*x6qH#&N1wP+4%Q&NckK#dtq;2V-AJmmD1_dW@v{
zgy}dp@$;E6T>MV@5XHnBGILUGk(d~{*2)+aC}`!XR`Oqsw$f%aVsC%*ARr+XK2e!W
z$A(w`fn+HGmC=CRAJU9TWnEsMLMG(Y*R4IY9Ob9u6uSIX^M3Tw0NV%B&&Uu`M#Dhn
zuRR8<JFu?V0#fOZ8?cY+$kH3oK{cj)L=}k>MJ7gTwkQWOfJse4*=JrRY`PyET7S^>
z_)By8NvGrox>39VPi=DbpQ3=|e+m`=4Wd|4o^S&_xdejMP0Wom!BA<yLsatDiu!_)
z1~fH9f=FX)#>v{u2{*OIJV41Hx%e+^)U6p|^w%HxsghnsFsyc-F^*w2zeX`6u|#nT
zMf7E(t)yeD>ImD-s4IQfG<#Z#&r9GlYzq_QQJaYbhh^&F<C2jb5aZ)^3C3nz7TTBZ
z2#w(Xheo4!g+_Q5y8cWhM>tq-e}_=|$+Z2Z{Q72OLy|JAszzl>`w~{8{8Xls8cV_Q
zH*{z#_-YRePDKm}W=QhW<Pvt53z@&s?0cIJ9)lQdVC!Wji0)ML!IyK~DRt}(Y+7GT
z(jtoJJHW_)<Ws!jSD0_FwohVePkn!0;@n6}iXFRBxgiNG;R6g9GbxO^uy38VX_Oeb
zG8o%CwXkD_jx}UHJ-I}#F@HuwU?v~Y^VCot@R`gD@_XCY;sNslqg-Xa%RKvMYomg4
z9U*y4mD^V`Ef;sURoI)BQE8F!yjW{X5f6;Ae899ztnv4S8XsP$@zQoE^-65IkzGn#
zZAzj97m6OVu5Jm%d%V+to;F>F{RT)?*nk@CSn`dX^R)qEN|b#D98TO6R^N?I!9n+_
zJ^t$36TQnQ?6xo~CeZh012`pGdf0$AIYFWMp60hp*}-{f8$VA~=;~dTG7(c{mRXc=
zK2h{K(=6+){@n8trKs>I^pxb_pZ8&<0AU4H3^C)j1$XB_Ug+99wNS@9bKopyF&%tb
z79-v%LB)7wS|Qp>M^CSyQ3SFFIIYD7F7Ki~8c-o|Q0$}k(><#6f+BrtDP@zT=Qxyl
z%$ewyUMXsVa2lEMub9IRdr25-!0P&C!=*?zSI6R$vywwEB$0J!4;OJok&T|xE-#v2
g!>l%J4$J!L{eoX-XJ>!o&p-O}FOFgjYaj^#09lCYU;qFB

diff --git a/data/imgw_meteo_abbrev.rda b/data/imgw_meteo_abbrev.rda
index 891716ed54eda98103be99834ef9eb59da47a557..23e8ab6c476ca560eb709a528ebc4a9a9cea8d7d 100644
GIT binary patch
literal 5230
zcmV-!6p`yfT4*^jL0KkKS)g=RR{#_a|H1$N|NsC0|KorE-}=A*|MEZp5C8xH;2A%C
z-BPb&@%6VGi?MBXU7Xt<^m}8ILyEhJ5zkPfgHDEl(G@5{D5U@>7Qn3rw9-Zzo}*Jh
zLroeH<jJIIl*K(TlhpL|WEAxWsM<gS)G<7oCXZ7YK*LQ5pdn0COo@VeYI#7=XwU!v
z4FC-^7>zLi$)FOM5eSV>(JAPSJw|GJnjVD7r=$ZzKr|Wv007VfKmtuk6F^hc^)h-%
zr?mw24K&HL6DC7OMuva@00E!?$OuwBLSkT;CYT7sz!8an0GePX0E|W<h{=EgB#00K
zXcHAr#Pu39r<C@MQ`1w_+M6a#5C8*4gF{URsh~74mUvA-Zs@+-ql{1&{a_0Z=leEX
zK|uWcc=QMy5`YI`KtE?bq>uyvL_kiIB0_G{q!19mAupFX1Li2fQ;kmFtl!Ye^6vL@
z7fdg;hHRI&W{}u=cZm=bYm*Jy$fOngo=$(;&C}3N>O{sUVI-`J2SlV{1fnr`TIuME
z%hucCYU;;lcaGHD`bcgu5KVa3h|PNwy)%>b=Xh|ruU9IHwka>hsUwC5#hDTj98Hbr
z(V$s=C$L)==oL~l3rN4Hj4THgxk1(Yz3+ltM2b~|?@@;I@yl4L?o355qjPv<wOOW2
zF4;3vRSG=XVkde|vRQ?oWfO1iNn!jbNs#5Z7=(8kt!e2~H6+wV6Kddu$hpm-X<|v2
zi%uLEZQ73uclS72oL%x8bciOqw}{PnCwgZm;GN;Z=DaynRxz)Uw30^*4~sG+A~>5H
z(W5}J{7+j{!it%Q;UoT9#z$y->1Ecb(~eq*;GJ?t{TG;O^7Uae%0a-O!FO8k&%b(H
z8AQaH^dd+yb>P<_GT;(|8@S^}%xz2o!bgFeG7A*jwWD3&IS3c*ddfj&X1U7juz)Nm
zB@@Fe%b7qHR4u5I5lRdeC?zSfbBH8NXrlC<is6=5k~q`2uST?KDS-w#I#6EEww10V
z5)2}cl>-}6;SYp92ml{&Kv9B%3<xBE3XFqLKo#`C6#Q{j3qV8<QxJANC$GPXZ~ooc
zctw%_9?&~ppe~xmVBpMkQDz5sb0LOxhLH@^Q_5>pKt=?L5MmI6ATff%1KFJ+eMA*$
z5WlrTC00zHW+Le;>%GuQLP`v9krQ;j2p0&NBBc%D+SX#nHeG^M)vDG}vVk@%qlv01
zMY4ee+_p-^r#sM5s7)nnyUs>UF-kilZv~pDh)`loz_T$@75Q-szkCZU;g#K**BHM;
zZ^@RN7*yAf8Uq4O&7lS}22{;;UBZ!?9hL7&#*Bjo(9HDxxbk%Xw<RZ(iyD!!u#$^5
zAtVo$hki@hg1WH(+RfNpO}iyWkEYa0lDkqz4*h{5$VmIS;s)i2ZVZwTnaMc!btSMz
zYh6-qI70x|o)U5p;6Zc}NLCG_EVgD<^x!0xFcisAP;Y_vk98%r4AN9Fj5L!bMG;)d
zLlt#pGe{<QS}>Yn(XRRitXeVI%7)>zgc|1!(oYy%c?PXE8XF0ggoVJO5(073*-~VY
z45)?O(niD7^JLgIQ+~Ke^<}FL>~g^h*_R@}A#-<zguOvlaB5B!N-3qJp&F`?M1*Z2
zsdkd3I4W5QK@x<4q~P1RvKVTDP5w->4+^kuGr+NydX)wu>Gt8fX4>hwSRze>DPU28
z1qevCsF-?Pq7bsH2@R!|TOc#3v6vj3#a${g<fZub*QQ9oaz*na1<ENRNePi8B@Vr3
zkX<&UIUJ`@STqvgZgIbt&*!c8RMRoktP-H5OKB#FSR41eBi>SiEH=Ob&7!EUGHULH
zgr-7^3yroLqM^4mp%)m;15Uck=q{yIywd=@5(*k9GutN>yH3es+Z(HF;Ut^7@wt8w
zh6XLJTg4DR>|!DBhP_!A6z?Dyb)*mXt1eH7c95$MeVUG5Fj!}ZGpEWIWbO>aL~Z4(
zypm>R%aluz%(-f2Tn6MD6b=FNW5zb(g$5YQc*;S&iAMOE@-O6G2!c@1at*O7o!%-y
z(oh)71}Ek&S;3MVAAvV}ci^ERxZPmfW^*wD(HvIOa4oGYY!Vq0A&Eu7g9#(3YAt?3
zLHhxL0vZvlv4}10?e89DXR0vW!r<OyodE0yX`AeI{yno3PL#N}M9E1dEy$s%&YN=*
z5qd#N3LbpFmq8?osCTcc%jZJLTPqW-&r7g91-pMY&#{9CJ9<2h{!19P^?BXA@fqiA
zRx`^vEN5)%o^FmJ#PYS0Nj(s?#2ppDM0MKaq7#ptb_PyF{w{%&4^L3^C!qs}Us4di
zOGcvK$LdT-9F6%0Pt~=6Vj5!@SPTo`0{jHx3t2R4kKIhGS=heRL_~xrK<lC)qzHtP
zPC&Ke0)D?&S7-{u1yx$)Pm_nEgH=ZQ_Ykpx#4_bxS96w8);jof5~?Xrz;nfTFOH`k
z?99xW1qu%!oM4?mOgN5s#^$`e?#cV;<C_k@j~&j<Ev6O07@Ja-DU4$jv3m9Kcz1%4
z379l@0vrTtg@{8018?5;ZgBeo#cnz7=*%}*Hfe)*wjLiDsj9F(FtNy}s^5SI0N=gJ
zt3h(jglrwcJRF<tK+@70*0Ul=;RIlC=9M0K+3K7|)}Fi7)OAj)#W~^`-<|+bM;nca
zL4e3-ON+OiiDQhLhgd3-XB@FqlQ{vw9B8$zt*jpxhS16fnMwd*JC=juo%*(RFL*{l
zJ+ox2uK4&alCjxvEINeVvPmS8fx*lXGmK)yyxP3Z>sK1}^E7>&96Kpk3tH9dTGpJw
zF|!2v3{8W*fK~6e_aexMiCjAAID*zmB)}mIVKdM(Z)acgpj);xf;hn$=w4r=CfaMQ
zoNl1t$(i?F`#W>+yCCFOC|!R3=F1}z%f+H6p{~5nP{;RpYVmQhz=9Z=C3U}=#fTU8
zFi0B8CP^uk>@HBcuC|cIRu=8thFMn3)SBh@W`XOy$0cNvzJoMw*?2QKQi%!#5e6c}
zWCUcgLYSZ(zgdw;B0-3;837qTw1qK1H1l_{bhp>8bM$p;l~q>7QmU^|s>F7STg*6{
zuH(<}9+T1C^!hx$V+1<7&I;OsF{;~xj9wD7o@WL<JQQtYs=&3l?7ROC9J!<7amO6K
z{6|l)Yu3D?sZ~`~eQb*M){k!E9N}!^?d|QGnmim+w=CI@h3)KH&C<!>(43!tkzs?!
zq=1o7q*X~Ipydgal1U__xOiso1wI?VK>_ZcMqfKVuZmSv(G0~?+4#C&d(M4liq(83
zSmH2u^jcP(+#;G9_N6x_lM2-;YObFdq$u5`t!RRk9Q9H)PiQ<yL&%3v@^=jRP93Fb
zENfk=Ys0pf7^*u-=yA|^oY|1TR_fC+wyv+D_!1OYsd7@3&?F5I)HDo0vRIR<T2W*+
znoSR0oNAcitXq|%CIggxjyj{U+9g$mI*bu^`%y4!m(SAX)XTGq0BVYaqFe^UQ**DR
z5aws6B{ESZO$v@EoU1!5S>-vcj2UHTvVyb>1A&2YNgOC&++9MlC{}j7h9?kcDMVxl
zoMWx<p(rK|R$2uQ6fB;ne#RRA3d+?)6bWR7l(f`Ji!nnZkvLVBYy2Q4HIs29g%5sU
znN&fj1rD1?5e>JaP2pIi7UaG{l|o$8d|J4|`VhTps*zP`iil?~%a0L$B;Y1RE(Et8
zKs&vx;6TS8utW92jLtD~$L!zQOp0nOigcZW$w|!x79=kWgn@+GHbq5}<N445nxl1H
ztS@s8VVLcGP|L9#$OBMQuU6VX)Z(W^Ltr2I7-zFWfaG|OxN9&}9wCqf$nR6Ed`aI7
zFgbD#zAJ6JM-aGYBoaW95eWyN5(O@0AD|cBC+yN7-4TJa7uU$RG)WE7V8Jx(VIrc4
z#Z-ihAtW}ACjWKpL!VfNTUI`jfHoodtDp}dT}g-*MDi3B@&SNb%Xq0-1((j$WGp1A
z1MOBe99KZaV2#1SFaqHZWa|sw3dBk2dtW!NU#;++VF-I=1W<Yg@E2WJgW$Z=kZ?T*
z`9$bWaJ>j%r|FP90}pFeWtWac#<kf`0d5u`%xM~{(z#=aV|XUz4&Tc|VMr{%HHIKI
zF7eho5G;?cu)V`Xa19f4=L}1Qonb;0BUA%Lg2@rE&|1=9-MFuxA9#-|5Vfs34=rM+
zkIl4H?&V`Ftt2ru9S~4~Fi;tcVWk61V#1>e$1LchRv3e?Rm}8pwSdGd_25QAY74hW
z2CZpMT5B;#MKh8GoG)M!;GismqOOy)3?S4(0~uOZ7Ithik>i2|)L+lA3O(a!WSfI{
ztOpmGJ8s~?911rh1A=>U?rC<A)6&U5Qm}4B<UrOO#n^qwV(eE3W>;vyswAT6^oa;Y
z3`P+Jmm4M2x;_TrFeDV#F=|Hm>^8c2+b&6&a*(aELX;%~l$Au)6((>B3<z$jCGfCP
zpyD?nFioX(QSufxn}oQr<Rb`010d<y7-I2*kOnpv5Iz|rpkOb00vL#FH7667x>M%!
zP<nrUgv<A6^g74e)7<Je$pS}DY_J?qY<$DlTLWj~r1Lm90hw_$5TvKIAdav#;)mb3
z{$8D${-x;q4mR#-fQNL6C?N3!D8#$0HFSL?hvGZ7OBsVI1+4^Vgz>N<O0Z(twS~d}
ztR2c43vah{%zlGXPnhp0v}6w@5kiX{762|qIt#wXS_u0hTnh)d^!gJIoI~vOIuFq7
z+7pCE;{_sFfwUi-x%yEIO^;oL)I>SNL17U-++B;f8R8&n0L{nA0@I)4h#n-DB6y~{
zoWdPcVn`<%hl!Lkfp`S0Ko^?<1oHS{+PoO3s>B$P8915m^Wp`eNp-Og7?(TUo1~0^
zAea~&UPvmDM9yD01!Xq`6hR{t78Bgf_5Mzh*yA^dX&0CXdBk}k^u%5P^nTIedA+mc
zj5)fGVT6ez*m4R?kvi24hnR{UWYr9%ya-KEzab_2k52V!+-|7N;j^cmB(zP;2SC_?
z_^dQ)?crNdyfO(u6eS977|2J?Z>^*7n#UMqnGAWGf(3xNz%Q<71_hifE@|tk8}tw`
zik!+yLkcu@pdCfcsR9W_KD-AlL^A?_^g$#vAxNUKBC0BiiW^lxsO~;LCwbrsF<mKC
zn%f1TMKgr}w2^^U8YLuB^RwwL=?n>TBTCf;P}LL^%2b#`AuM26NwgC-VPXpvR27oP
z1DS!qD1ic$x*Uv@lHqg$3V56h3ub#Pak(<G2-az9NJ0g|>SAauAXmk}Q)^&u4hFPA
z$ZQ85QfmP@q`x>Vl9s5ZnIVHj7^755Z2<(7!J&YbDX6w1OiL5U7g`2hG$l~)mwRm;
zvlyDwVTD<WRS|`V6#^3|A2%&+*k{RU7?ov9DPf?+GF<@AEDfe`W+K?8(t<%^D~v!u
z&R`))s)iCog#?6XAkt8gLZRBEZ>5MuXozd&F6IbeaJ%-Vj?%UpoRcvpWl>dBH<Cj)
zpo=tuGgTxd8Pi}eYTFRO*TztW&j5tE9wvN8q`?zUW0_0Z>!ag7HFh2VJt^WGYwoxS
z_sAil0v__8AxY#%xR8*LkdT9@&XM+$3y;=u3;^0+p2fHwr=aSSagw0)iF6>yqp}hR
z8lNX;Ppcq^5G|4>{f96i4P#XhH|;KFIqrJY$c;2oQFCqst~^YK>LD2k3=TuFIG`CP
z%ibG+f)7v})p6e-Z6x_y7C?#;(9&8+VlB5=3oXoKngCtwfa7H@KZx@?$N?uH^Zh>1
zm?N0u)QCA`gtvMn#=whlgaRo0Cqgv>c!z}#Ab0wu<O3`oSYr+$MKowFZ0TSvRS6Vf
z5u_~C0~JD(1{1QDl43ts#bWMe5d^O+nv8YpJ)fftGX&HqFauOX4rm_J2Y6K;ksQOk
z#NIxQzbHV>K+urr$U}ho{(q-t(B;oFq3B@(5J;10=`!9i>P;=$9o?qS!RIT7-mV?q
zx&df8Cs1=DpvB;}ATXc^Oa|d#=6N9!_#t94Zx)kW3|r35|D%Mpgrx?MzQD99AwKQM
z)bx=_s!)<aQ1F<T1wx8SS&*nYji4E-#>#_*@=8e;eHATDLbZoxhN1;gO`#xe?oR@h
zb}3=sV|trFfZeVCIYonqp!^4*?l4mn2!{}W9&aT4Iez!2qxFmCr{Vr@y;}jg$7Wp!
z`he>LfaG<=dpiLGthpqzRd^O@l8P;9D+Dw{VL?1v4}&jQc)Gl!AqYrThC$^~)>0;J
z2ngpIffQCFwjLACT$?hovZo4aNRmn<(S{aUlvpfWm)23>d!gp;e9ewzKCsYFd$Vnj
z;DBxP4x}}b_ywK1s)Z+t!BrVxw5EzgP=V^zsA5B`F2DruAa5dc^|-F`_u71+@0Q%D
zLlwv}9E)IshfTxA=XZ6(K@M<Z5FNN~4`A;pD4_P7PGJn3jtPB4FivP5Q0Z`Aj`$it
zCBPyF1P-gHp!9)}ai=1%<D;O@G1cAgrsNU!3aX5$=)1MhJ7{hqplz%b%jnkt8jtoZ
z#B@MUDaC~9I3W{S>Pj8f^F(%@&n%>M#4=JIsi}&hfkttQ0t|qJU`rT@_2MApKA1Bk
zMGajk_abmQLO}98+<Tv7`2F+Ld>d1U876}Itfr_5X9Lt>alC|8o%NV0ASFtKvup|F
zLx98qYA_BCWdY9UQry6ip{fU>z#a{QR8X@O2?sBGNe-L=*llZ3K@8CPAW`>VA;3Jn
zt>QPc3=qt>nzAH?8PaPz@b6;qC@?3vzzT>0W{H5pfuwJ&V9*%a<^q5e77G@b%R~WW
zEDSIZbIka{L4zF3j3`c=4hq@{6_Ic_Sh2OCx3B@ULlDp~NKugH=Dsoz<`7mcPOT&v
o77TP~+de~C$J-jfD+!N9?1^X;I)7i>|A+q<az!{$kf3x|SLtO{ga7~l

literal 5237
zcmV-*6pHH~iwFP!000001GQZLQyWLJ=WKH}#NRO37!#7U!Lj4B4Z=?;=ZCv`1jWEJ
z0u>UG4=SZvv<uoutJT|GS&^z<)%(jYuV;2=rsvy&Q}?P&!ghb=TTgdSch4?Il~3<x
zKg~`}O<kY5apU^b^=bHZYWleG;m!N0sT<ee;V)A+re@&(b3u49YK4C4M=kg4EcQP`
zqj}nLYU)1x{3{2pZJi}iZ<zL7C-uXDAG_%=cAe1eCgafU0r_<M*BrZQ!3O<6g)VG$
z({Vq73B&==|7*P+ZZh7t+Z{J{yQv?KlXToAo%KDx-~KIdPP>uEVE48=;V>RM+q)8Z
zPLNZ-o4`a={Jt5#Z^Kv3_^N`R+YbDR8#{wII_ttzFIWp_J4pJ05C0sH!62FhA^aLf
z@dQ@hsksdI)QC$`m#!QEmj!U(!c<0)*Y|&uAYX4!+{8)ak=qCPo(<#47-kfClb}Bq
z^TL>ms2BBo0BAV5uZ!6Ah8kvB!h{z+eYiCV5)DW<o}WzG6TQ2Q0_+520!zOL{7`GR
zYG_Bnm(4hi5`D(<!V}jdoD~@}iKD(g)^iCo3fwdvI{)~4X(L=uE|C)hm&zuryoD1N
z1n-mIRNT6X3;cFKB{LzK(+`qvKs2au8?<JlaS|EZJ(b&ml~Z-q=QUlK03&~RO&^>H
zG#>Dz_tYI^?XsQZPRGB9@-}a>HUNu93F!Eh8ZGRAW9s`|t{TX^Z|=|;y7}Xi{88$7
z5XIxsFh|VqO(Jkrq3d`N^>OX+SB-TM#WL%`N0>34L%zYMjI;DH@|-T^c88pcW;ja#
znO9()?Fa428B74AHlfDOroSk0xQR3Nyvzp67T!!^fnKpT%Mydg80;{&%^dQg0D{ez
z$=C$XwuTv^gcw9!kiVPjp=+El*t2<&)QTBvMPgA|!_fcuyF)G0xL^*V25A&CGef?>
zRh3Q5h|ix|ADhuyF_c;y!#>1VKl43&%q5x0ZIA<=RmBnRMQ0Ha=rGKzCK@@XA=lNX
z2#L@2yp=fAzehki(<o`@yEPVG16je%YTze<+h&rckA9o1)H0WsDrA`SVAb$!b}&wr
z2Oot~|IUmXkaP0_OtB4Ttw1jUdYBeCh1&=2ga{N{>63%N_u!aGL5F%5QQ8OTWas+-
zFu!yN>lF@bj2Hfp@*raP!O510nXeCXpXpI=cs0*iV?Of>8ZFc!7##u%66&;F$uBY9
z5SZXDB%V`?!hps`lQDQFV&Em`Wj~lC6Zd=Pv{r!IbY70eD(+d0NE?d^Bm;`(Cr+~o
zl37o<qs@@ak0Q>)TI_pNMGee_BM`FtTP}&y`g(KSVZ3i7p0p)s$tyWZ*M9`MeiU`t
z(vPB1ACwX_xl1FuF+ZNda9!BUS{<fMHYwNPLBT4%QLtFGk?t!fp`OlY0-nqxN7r+P
zE=g$LX{$kNWQkf|SsU7d*q5wDD2%>nxz7xE^O9c~kj2a9qAYYH&rB#|u4G1v*$CuQ
zGg8l~)-#}Q-5GP*@mgsHhD|b@32g}`4fgMxM%~oyr`)Z+V{Q11uRL7&Dot%V^tQ0E
zO^-GVt*)}lxjH;j_n*okn9Mj?KjIM|Fw$WHYaniH=xlv#956?uvT8L;qo0mWwrs59
zql$&Ky}P|@XYJcrRU2!sc4TAi^H%DzwrfpGqsq@8_N=VBZCSuNv9b2*2R7En%10Zk
zeqbHz<R_b)DkrtNh4tY>X~WJcTjcbid1U3-M%vkGSXqbrwp9V^(8fC2vbH*^*jZM&
z@74}&t$<a~%BD2CwWC9eoObJV3uo`7V&&{u2Kso^v}yZ(rDhvxzhc{v>eeTl)T=xD
zHv0yaZR!VyJ64Wu5{I?wCkw0g(YlGX!y}vV0L!*hwOY;QNq}XWVXa=VP2=dW*08eb
zHVFdDrh&k+P5!86mBvxi%BpW!Z3I|0w^iRcvUXUl{A6Pt{A^R*`oW$}Y3dC-$F^Jb
z<Eou?Y}?5CiS=Xxmdz9zKUZz1U}I;aY-MqeyI^41ym4b^)6Oc}SykJy)c_;6vi5A7
z+^8Pdb`n@N`8KMbtgKqY=DHe<lcq%hj_bCw4mh^CH%l9~U2c}jHeGA(SsU$DY=#Uh
zo7Ob<cWkl%md(bShd<cl-aNEgYxD5HHjR_5s%>2-Th<`sWUFaqRevsE*_`1C92!>E
zflYT$s@5>$<Zy@i|Hb++<Y|sxjd`T9hf>L+&T2TryvBf-^N87x+L-mEffAiO$bM7W
zK$*_68G#|jPBt-O$%H6h)<qe2@uFc)^C0CZZStakNVFTq49-$ydUpl%mkne=OM)0g
zqrgvT(q<O-VA$`1AhHV+4tb6tq~=Yp_J+;`vg%1hDo6<Y^>H_uq;4;b$OVDVOC&hO
zx_gWZdwx!REK~JMou(V7B!fp+hNFzcbGK$dD!@IO=;IL=F~c3EQRt?imRw4h9nrN$
zq^1y}_C*DR`DW71iaWiC*vbBSQ*_N(BOlVwUY7PY*AIBw+Z5GQ57J_g70q7BrR8?j
zNhag-n8;Etx~=G~XvK|W$D+8agXpzpMD`;Z1W7?}v@R$0_D%CDW#0vX%I``A?_18x
zt6#Jt(DT#o5W`jl4D@vpISF`@9#Jz^G;4?23Nt7m&A_!Q^vPAY$;+j_bJiQf29AO@
zh-_oM#7jIFqss#!A7Nx>S%|xGf@>AoU7c2|sGOJz8XLrKdXeMx-@YroEp5In{rLX<
zdfwEDJi^#0k$=GW(LA7~IC2-4xZnP6q>gjyOzif%t#;&bJ3es-ZbE7`phc+8`RkPY
z_Q9<*bboKTXNh<Y`sjH&h-25wA7%%!pCm0}`i#=<wYGQhz|)k1+55%ps>Wu)@twpC
z2R*+v@`H;`niH=ly#WK$>EsVzeRc;?+`^;>X=eydDh6AN!&WC6#vFXnb0@w;n(eqj
zPom76M;9_Ac0o20ByJ@#B<V*Z8RGdPrI#R-hrZW|+!*{zl1k(|qoh^Z2os%HR*3gt
zrWZk<><uQ*lgVU@?<uod5UdUeSK{66xRbHg*>%rl%66~k!c^h8MRX^QbMpQ3AU4rw
zs7%uFfNk;((l29P4D)5P(P657+zOH??Z7F=sKdcPr(PQ|J-j*cuo_(p7~^)X17uo!
z2eREqlo`7Dh;WPSJqFON#{jbP7(f?^2zLpg9etLEzD2Y_>X#@3BHU)1uF}df?XK8z
zM13IdIO5DJ8;-jY>W8RzxEA1<#4+7l55&c^VK!AMcf`_x#kM45yEaQYiTY`WH+T^A
zlc?vWLDW|`v(z2o&TpmVcGv<cH)WzgA-OJOyML41)j?lu>Buv*JpVjs2h5rLHO5n%
zhNmAoe@@A7w};Sr<hw*a8ErtikZIGQ?<OVuu*qgk7t&<O)eZ8<^_*?MO&A1O(ro@H
zrj67N;jWFD_#RR{9k~!ew*3Dd2A|=m?z3@-bAc(z0>7JXhC205+|h~mGZ}02M|cM`
z?pJKokZ;*-KBl03#YPk)^VDXjLK1R8Z!=V}5a=$PD~1)dd(qI->9;u=%?7PcP#3hd
zqZX)lLA{qm;Eae4``D+=3NoW)X(hS>naL^q<|Buo3lByS_Z&F|U2h1=zDB7k6mf^a
zDB|iOpN_b($RX&Wf>Y`hjdG$;#BBtlh)W1Y5qA%aBCZ@bg%?Y63c6qrl(M`Cah<>k
zN-ya0fKkMqfuI~J6mdhqDHV;Pt{hzh2#StrsuW%)%oia_{5e5X_Bmlk?gZ8Ga|$Zo
z&oGK22~@0S6j7rW6lH8tk>>=ZUsQk3DX8quD5BafD2ilIZCy~5u}~#Fr=S8lr_{GJ
zN=2h6qod+Dqcl|tDuoM*Vj-yT%_ySY%_ySG%_)0|+(l8F6BOy75;mtCEAm1GYe7-8
z0F|d1MO39Zg<I!Wg0i84HZ@9Fqf`~QgX+tivUi|T4t6w3U8Q_dDWYV|D56rFOH~wz
zXE{L`9F>DPrL1g@DE`iJ!j?u+1dl4bf}*G;s^@Zws&c4$n^REHmQz%xj7qeOA}X_j
z(oiWo`zqyNPo=1C85LV^#jzYqD$ld@d7MH_&V!D))hF4fAl*t!8!d7c^GE`e5fg5u
zo6_^mv@EmB3Y>&1rOgiSk3=R2aEB8RD$CQw^L2=E@?Wo|Ecuxx(I*S;?up1y3L-@8
zom1|kXdGvNxAKx7PXj*(NpKeQf^>}I6ig;XKGrKDcZMO4yROA-33r0F-wOI(&?bdH
z-VDN1Sh-4{qoLE8oA^ov;f`QH3dCf}i8<k}K#_Ea4JJg?!8Hd~h>(bn1o%Agd!Edl
zJ%@m=(;}EJpv?407yY1rp>JB$Zi}p$q>-5L*LDvrmpr9WYJ5_XBjxc1sUxvfXITQ#
zJu*Y{S1Y~)+8cz!&~bwpy+DQQ><CLL<%QS>d!*hJxc#=Te6%uGBOx?FD|)!N3qo_~
zhB{PK&MDD{U+wJF|2;)-ntEQ(x3@O}a1W36qmjM2&c$u#>Qu9wti**rX?a}tU8m=N
z_Iu8Nty}v1q@L|t5$549ka`SbGt;RI<E)41+`)x1KoO(o%f~ocI8gwPF@V%j6#gRS
zpGOx(%}JG|X#Aw-`>6t0mH?osU_C`s0MjX)0(g3bzap)W6-yM*D-41;cuq3#ea~^v
zVAq|}h7%@V)}pMpSW!fo6@V1bQw##1;PA?i0(>Sdix`Ow`L(>oSzSoMT(9u@;Of2$
z;(>z{R+V<Rg8(&b2ZKQJH)|hAIpu2D`*oO}L#i6Gh27Xi<;0yKbcf1`^ZTTn=&-_}
z^zq0ssVN#+3r(^eUI4S8%RcpdZkZ2N9-hlJ9*aE!c_i7QuyrxNF1XaD*s&zY#FCgJ
z7F5`D9ZIZ5@0Z$&4o0WVf=tc)lRhRGjA2X>5IKTs*88jA$^thpEjQC~ZpfS^X*M@x
zAwEf;^`5yIvvUVaK2z|V^KCWPNg65UnwdLa`5O)lv!8G)GEhh!E8ujMltG!`8d2)@
zoV+qB+XkM3f)in!>LIH6F?`un-^#$g@5T-sB<C(1&h}RG^xF-edo@(ZyjhnldPO~;
z^6IQ~P%Tv@pP(;$wI&JVJc=DsyLEn`ejt&c|39AT(e^NoAx?GjTCeOA=4)h7BWtzl
zyJX=ZZJZg~X9$Cxl0&aS&<7b-=%z#OMi6w3NsNmhbb)-f0k9zQq&NESH8}>8?>B{o
zV)|oQ)|oS|i;MPG2EM)_QhZbw24L{*{{)?5+Cq^uD0M{XTKRZmxk6QylvJ^ok-phL
z5W5yT2o?Lpj=)Vp;u3g}!JxMeULhRUV>fjYDL#`UFVW(2L4>wa%?-I|F^&YXFXx$r
zD_YQ1dweP3F3$q~WWSeY49_KeQXDp~O5l|GOB$TBNJ>DzsNdHT_-c9}Z%L6Re8H^e
z>JdcXG{KiG={{zl4=Nt!@1}vbcH8VP2CztOSuh8d>}y3pq$~n#^MapX%N7v2T3iy4
z4o}$d+5cDodMu^Fd^ARgGBzoE5Re4tUF@W{jUY=B!7^v0Y+Vj}aS6<kuOgocxW_{g
z(_!-sk|oJivc2IaqdrN7G9^!VEd^oiqNO)>!Xd#a1q?I*igak{0#3{}k0F(LSYyZ{
z9Eu?qOgMc?)Sj~7MuR;OFj}Y4W8`}RfrEIY$(a$TfV)W5MU?;8qyUJ&VwNNXJf!-h
zuRE(bpqkecNxh+M5h9BziV#`3OW2ZbhN{>$jh2j3FtaPr0gr83tJ_jf3i#A`bPcIO
z?qYREZ%~_pSa8v_`yy@o+Qe2e4kTpN*R!(;>q=ld7)+jpg5BSZQ(rJ?m@X+Ps-*59
z=!1WQa7=dxs=(3mt%6=P96IVjs{2CEp3DAC>QD_&NgG&59-Qi;AoU{YFk3VYDw5>K
za#PLKEg9_?y{w_SOf}-w#nc~$XT%ra3DG)Mkw?M+yz@Lv(ESe1fMwNLHjR)M<8+-h
z%~uAm85LX!0(0zlx)B$kpX=z%ZcNL1W&>-NCNGLtr<NO5jBY>yp=~p-<+ZYHI6y8t
z?a(>%Ng#6(M??O`=9_{xaMGF^y)A5HYgIJV=AUU)wl#W(W4|n6&>Dk{)5D;&DOvia
z?$*^@;|eV;S&F6Qy56#|wqj89g51c4Hm^>1I>c)YBJk-37*fd*1e5eoCK&_wwaAlb
zHge^I=s${$sW3wP;<hMiM5mLut#n?}>&ox-BL^o4QJcKGEz)vxPx(g=zW2SpNnSy{
zLHDvY2%27Yq8SBzC7+H|f?k1JV3&A4wS@b(Vz{ZS`}z{!YAW@9cS##L;=*x+;i|)6
zBc_WO{b2TNS22cT&UzV>WB=1073QBIocZSwamYWufA6@U<{NL`J4rxpQl)PXvMO<L
zHTz1S)rNj>l68;^I1+Lr=<@<-KfC%we}br9bTBa|iT)~*0bkfaA7HZ(|3br;zpVRm
z@pP(|gNIa3+0h~4uX1PzeECZF%Lcn5etye3kEM#~WWhzs1|sW63_>oVs<ENiHM^&t
zJrDM5w;@Z~dj-xZ1lmL1+1xL}gUeNsoB{N+F)PQT@7W}Fbe2GlL>3`F-#jWP!_M!u
zN&Jw%^&)-Ew)_iyCnj4Af}ro&9!oDRkOW?y`og|a0l=Tx9V@xN{K~?D&PJcBifzg3
z1L#s~ahI12TvmWXZ0CbgFCfJYH!E*o+(QlbZ(%O%2M~vo(v_aGE!R!8E9o1ug^HBy
zNxUWP9S|j4c-Js{^faR0{I!$gxjX-nyYnY{Q>Am|(oK1sV=t1BH-8-TBGr>7%{$wW
z_mFEA2xY0<vTw+wKXp6ZcSAp!%3ct<o)d4K4|_erV9(@_;I@38;oW9@N>e&)L)6SV
vWW-s|P2IQWG3io<&&1K_Et|s~___AK@Si{a_=AA{MSlJd^|)_fxK{uG?2suQ

diff --git a/man/hydro_imgw.Rd b/man/hydro_imgw.Rd
index 66904cca..d38b884f 100644
--- a/man/hydro_imgw.Rd
+++ b/man/hydro_imgw.Rd
@@ -25,6 +25,7 @@ function that shortens column names}
 A data.frame with columns describing the hydrological parameters
 (e.g. flow, water level) where each row represent a measurement,
 depending on the interval, at a given hour, month or year.
+IMGW parameter columns carry a \code{label} attribute with the original metadata label.
 If \code{coords = TRUE} additional two columns with geographic coordinates are added.
 }
 \description{
diff --git a/man/hydro_shortening_imgw.Rd b/man/hydro_shortening_imgw.Rd
index a477bbcf..32de1490 100644
--- a/man/hydro_shortening_imgw.Rd
+++ b/man/hydro_shortening_imgw.Rd
@@ -4,21 +4,18 @@
 \alias{hydro_shortening_imgw}
 \title{Shortening column names for hydrological variables}
 \usage{
-hydro_shortening_imgw(data, col_names = "short", remove_duplicates = TRUE)
+hydro_shortening_imgw(data, remove_duplicates = TRUE)
 }
 \arguments{
 \item{data}{downloaded dataset with original column names}
 
-\item{col_names}{three types of column names possible: "short" - default,
-values with shorten names,
-"full" - full English description,
-"polish" - original names in the dataset}
-
 \item{remove_duplicates}{whether to remove duplicated column names
 (default TRUE - i.e., columns with duplicated names are deleted)}
 }
 \value{
-data.frame with shorten names of hydrological parameters
+data.frame with short English names of hydrological parameters.
+Existing column attributes, including the original IMGW \code{label} metadata,
+are preserved.
 }
 \description{
 Shortening column names of hydrological parameters to improve the readability of downloaded dataset from
@@ -26,12 +23,10 @@ the danepubliczne.imgw.pl collection and removing duplicated column names
 }
 \examples{
 \donttest{
-  monthly = data = hydro_imgw("monthly", year = 1969, col_names = "polish")
+  monthly = data = hydro_imgw("monthly", year = 1969)
   
   if (is.data.frame(monthly)) {
-  abbr = hydro_shortening_imgw(data = monthly,
-      col_names = "full",
-      remove_duplicates = TRUE)
+  abbr = hydro_shortening_imgw(data = monthly, remove_duplicates = TRUE)
   head(abbr)
   }
 }
diff --git a/man/meteo_imgw.Rd b/man/meteo_imgw.Rd
index a300d83e..b5bbe120 100644
--- a/man/meteo_imgw.Rd
+++ b/man/meteo_imgw.Rd
@@ -11,7 +11,6 @@ meteo_imgw(
   status = FALSE,
   coords = FALSE,
   station = NULL,
-  col_names = "short",
   parameters = NULL,
   ...
 )
@@ -40,9 +39,6 @@ e.g. \code{c("POZNAŃ", "POZNAŃ-ŁAWICA")}.
 For \code{rank = "telemetry"}: station name(s) as listed by \code{stations_meteo_imgw_telemetry()}.
 \code{NULL} (default) downloads all available stations.}
 
-\item{col_names}{column name style: \code{"short"} (default), \code{"full"} (English descriptions),
-or \code{"polish"} (original dataset names). Not used when \code{rank = "telemetry"}.}
-
 \item{parameters}{character vector of parameter codes to download.
 Only used when \code{rank = "telemetry"}. \code{NULL} (default) downloads all available parameters.
 Accepted values: \code{"wd"}, \code{"t2m"}, \code{"t0m"}, \code{"rr_24h"}, \code{"rr_1h"}, \code{"rr_10min"},
@@ -55,6 +51,7 @@ Accepted values: \code{"wd"}, \code{"t2m"}, \code{"t0m"}, \code{"rr_24h"}, \code
 A data.frame with meteorological parameters where each row is a measurement.
 For ranks \code{"synop"}, \code{"climate"}, \code{"precip"}: measurements at a given hour, day, or month,
 depending on \code{interval}. If \code{coords = TRUE} two additional coordinate columns are appended.
+IMGW parameter columns also carry a \code{label} attribute with the original Polish metadata label.
 For \code{rank = "telemetry"}: a data.table with 10-minute interval observations (not
 expert-validated). If \code{coords = TRUE} columns \code{name}, \code{lon}, \code{lat}, and \code{alt} are appended.
 }
diff --git a/man/meteo_imgw_daily.Rd b/man/meteo_imgw_daily.Rd
index 8d51bfbb..1eaf7f17 100644
--- a/man/meteo_imgw_daily.Rd
+++ b/man/meteo_imgw_daily.Rd
@@ -10,7 +10,6 @@ meteo_imgw_daily(
   status = FALSE,
   coords = FALSE,
   station = NULL,
-  col_names = "short",
   allow_failure = TRUE,
   ...
 )
@@ -32,11 +31,6 @@ database and thus providing both names is needed
 (e.g. \code{station = c("POZNAŃ", "POZNAŃ-ŁAWICA", "WARSZAWA", "WARSZAWA-OKĘCIE")}).
 Stations' IDs (numeric) are no longer valid}
 
-\item{col_names}{three types of column names possible:
-"short" - default, values with shorten names,
-"full" - full English description,
-"polish" - original names in the dataset}
-
 \item{allow_failure}{logical - whether to proceed or stop on failure. By default set to TRUE (i.e. don't stop on error). For debugging purposes change to FALSE}
 
 \item{...}{other parameters that may be passed to the 'shortening' function that
diff --git a/man/meteo_imgw_hourly.Rd b/man/meteo_imgw_hourly.Rd
index 7f116074..630c0cb3 100644
--- a/man/meteo_imgw_hourly.Rd
+++ b/man/meteo_imgw_hourly.Rd
@@ -10,7 +10,6 @@ meteo_imgw_hourly(
   status = FALSE,
   coords = FALSE,
   station = NULL,
-  col_names = "short",
   allow_failure = TRUE,
   ...
 )
@@ -27,10 +26,6 @@ meteo_imgw_hourly(
 
 \item{station}{name of meteorological station(s) (character vector)}
 
-\item{col_names}{three types of column names possible: "short" - default,
-values with shorten names, "full" - full English description,
-"polish" - original names in the dataset}
-
 \item{allow_failure}{logical - whether to proceed or stop on failure. By default set to TRUE (i.e. don't stop on error). For debugging purposes change to FALSE}
 
 \item{...}{other parameters that may be passed to the 'shortening'
diff --git a/man/meteo_imgw_monthly.Rd b/man/meteo_imgw_monthly.Rd
index 3d3f2224..31f03d4e 100644
--- a/man/meteo_imgw_monthly.Rd
+++ b/man/meteo_imgw_monthly.Rd
@@ -10,7 +10,6 @@ meteo_imgw_monthly(
   status = FALSE,
   coords = FALSE,
   station = NULL,
-  col_names = "short",
   allow_failure = TRUE,
   ...
 )
@@ -30,10 +29,6 @@ It accepts names (characters in CAPITAL LETTERS). Stations' IDs (numeric) are no
 Please note that station names may change over time and thus sometimes 2 names
 are required in some cases, e.g. \code{c("POZNAŃ", "POZNAŃ-ŁAWICA")}.}
 
-\item{col_names}{three types of column names possible: "short" - default,
-values with shorten names, "full" - full English description,
-"polish" - original names in the dataset}
-
 \item{allow_failure}{logical - whether to proceed or stop on failure.
 By default set to TRUE (i.e. don't stop on error). For debugging purposes change to FALSE}
 
@@ -52,12 +47,6 @@ SYNOP / CLIMATE / PRECIP stations available in the danepubliczne.imgw.pl collect
 monthly = meteo_imgw_monthly(rank = "climate", year = 1969)
 head(monthly)
 
-# a descriptive (long) column names:
-monthly2 = meteo_imgw_monthly(
-  rank = "synop", year = 2018,
-  col_names = "full"
-)
-head(monthly2)
 }
 
 }
diff --git a/man/meteo_shortening_imgw.Rd b/man/meteo_shortening_imgw.Rd
index df754523..deba7a17 100644
--- a/man/meteo_shortening_imgw.Rd
+++ b/man/meteo_shortening_imgw.Rd
@@ -4,17 +4,17 @@
 \alias{meteo_shortening_imgw}
 \title{Shortening column names for meteorological variables}
 \usage{
-meteo_shortening_imgw(data, col_names = "short", remove_duplicates = TRUE)
+meteo_shortening_imgw(data, remove_duplicates = TRUE)
 }
 \arguments{
 \item{data}{downloaded dataset with original column names}
 
-\item{col_names}{three types of column names possible: "short" - default, values with shorten names, "full" - full English description, "polish" - original names in the dataset}
-
 \item{remove_duplicates}{whether to remove duplicated column names (default TRUE - i.e., columns with duplicated names are deleted)}
 }
 \value{
-data.frame with modified names of meteorological parameters
+data.frame with short English names of meteorological parameters.
+Existing column attributes, including the original IMGW \code{label} metadata,
+are preserved.
 }
 \description{
 Shortening column names of meteorological parameters to improve the readability of downloaded dataset from the danepubliczne.imgw.pl collection and removing duplicated column names
@@ -23,9 +23,7 @@ Shortening column names of meteorological parameters to improve the readability
 \donttest{
   monthly = meteo_imgw("monthly", rank = "climate", year = 1969)
   
-  abbr = meteo_shortening_imgw(data = monthly,
-      col_names = "full", 
-      remove_duplicates = TRUE)
+  abbr = meteo_shortening_imgw(data = monthly, remove_duplicates = TRUE)
   head(abbr)
 }
 
diff --git a/tests/testthat/test-hydro_imgw.R b/tests/testthat/test-hydro_imgw.R
index 67c00e79..9a7c6f8e 100644
--- a/tests/testthat/test-hydro_imgw.R
+++ b/tests/testthat/test-hydro_imgw.R
@@ -6,9 +6,11 @@ test_that("hydro_imgw_not_available", {
                           station = "not available", allow_failure = FALSE)))
   
   # monthly download for a single station:
-  w = hydro_imgw(interval = "monthly", year = 2010, station = "WIGRY", 
+  w = hydro_imgw(interval = "monthly", year = 2010, station = "WIGRY",
                  allow_failure = FALSE)
   expect_true(is.data.frame(w) && nrow(w) == 36)
+  expect_true(all(c("station", "riv_or_lake", "H") %in% names(w)))
+  expect_identical(attr(w$H, "label"), "Stan wody [cm]")
   
   h2022_2023 = hydro_imgw(interval = "monthly", 
                           year = 2022:2023, 
@@ -30,6 +32,7 @@ test_that("hydro_imgw_not_available", {
       testthat::expect_true(is.data.frame(h2022_2023d))
       testthat::expect_true(nrow(h2022_2023d) > 50000)
       testthat::expect_true(class(h2022_2023d$Data) == "Date")
+      testthat::expect_true("H" %in% names(h2022_2023d))
     }
   }
   
diff --git a/tests/testthat/test-hydro_shortening_imgw.R b/tests/testthat/test-hydro_shortening_imgw.R
index 1d734a98..e00c835b 100644
--- a/tests/testthat/test-hydro_shortening_imgw.R
+++ b/tests/testthat/test-hydro_shortening_imgw.R
@@ -10,10 +10,38 @@ test_that("hydro_shortening_imgw", {
                  `Przepływ [m3/s]` = c(288, 413.266), `Temperatura wody [st. C]` = c(NA_real_, NA_real_), 
                  `Miesiąc kalendarzowy` = c(11L, 11L)), row.names = 7165:7166, class = "data.frame")
   
-  eng_full = hydro_shortening_imgw(data = df, col_names = "full")
-  eng_short = hydro_shortening_imgw(data = df, col_names = "short")
+  eng_short = hydro_shortening_imgw(data = df)
   
-  expect_true("River_or_Lake" %in% colnames(eng_full))
   expect_true("riv_or_lake" %in% colnames(eng_short))
-  
+  expect_true("Q" %in% colnames(eng_short))
+
+})
+
+test_that("hydro IMGW column labels are retained as attributes", {
+  data = structure(
+    list(
+      `Stan wody [cm]` = c(258L, 287L),
+      `Przepływ [m3/s]` = c(288, 413.266)
+    ),
+    row.names = c(1L, 2L),
+    class = "data.frame"
+  )
+  attr(data[[1]], "label") = "Stan wody [cm]"
+  attr(data[[2]], "label") = "Przepływ [m3/s]"
+
+  result = hydro_shortening_imgw(data)
+
+  expect_equal(
+    unname(vapply(result, attr, character(1), which = "label")),
+    c("Stan wody [cm]", "Przepływ [m3/s]")
+  )
+})
+
+test_that("generated hydro date column is not renamed", {
+  data = data.frame(Data = as.Date("2020-01-01") + 0:1, check.names = FALSE)
+
+  result = hydro_shortening_imgw(data)
+
+  expect_identical(names(result), "Data")
+  expect_s3_class(result$Data, "Date")
 })
\ No newline at end of file
diff --git a/tests/testthat/test-imgw_meteo_abbrev.R b/tests/testthat/test-imgw_meteo_abbrev.R
index 211a4840..43e2d922 100644
--- a/tests/testthat/test-imgw_meteo_abbrev.R
+++ b/tests/testthat/test-imgw_meteo_abbrev.R
@@ -1,5 +1,5 @@
 context("meteo-abbrev")
 
 test_that("meteo-metadata works!", {
-  expect_equal(dim(imgw_meteo_abbrev), c(254, 3))
+  expect_equal(dim(imgw_meteo_abbrev), c(286, 3))
 })
diff --git a/tests/testthat/test-meteo_imgw.R b/tests/testthat/test-meteo_imgw.R
index ffb0fb0d..7c3a8698 100644
--- a/tests/testthat/test-meteo_imgw.R
+++ b/tests/testthat/test-meteo_imgw.R
@@ -18,8 +18,8 @@ test_that("meteo_imgw works!", {
     x = meteo_imgw("monthly", "precip", year = y)
     x = meteo_imgw("monthly", "synop", year = y, status = TRUE)
     x = meteo_imgw("monthly", "synop", year = y, coords = TRUE)
-    x = meteo_imgw("monthly", "synop", year = y, col_names = "full")
-    x = meteo_imgw("monthly", "synop", year = y, coords = TRUE, col_names = "polish")
+    x = meteo_imgw("monthly", "synop", year = y)
+    x = meteo_imgw("monthly", "synop", year = y, coords = TRUE)
     expect_message(suppressWarnings(meteo_imgw_daily(rank = "synop", year = 2001, station = "blabla")))
     leszno = meteo_imgw(interval = "monthly", rank = "synop", year = 2020:2021, station = "LESZNO")
     testthat::expect_equal(nrow(leszno), 24)
diff --git a/tests/testthat/test-meteo_imgw_daily.R b/tests/testthat/test-meteo_imgw_daily.R
index cf277f60..2ae5b317 100644
--- a/tests/testthat/test-meteo_imgw_daily.R
+++ b/tests/testthat/test-meteo_imgw_daily.R
@@ -80,3 +80,43 @@ test_that("check_encoding_in_non_synop", {
     expect_identical(nchar(non_synop$NSP), nchar(trimws(non_synop$NSP)))
   }
 })
+
+test_that("current IMGW daily labels use short names", {
+  data = data.frame(
+    "Maksymalna dobowa temperatura powietrza [°C]" = 1,
+    "Minimalna dobowa temperatura powietrza [°C]" = 2,
+    "Minimalna dobowa temperatura powietrza przy gruncie [°C]" = 3,
+    check.names = FALSE
+  )
+
+  expect_named(
+    meteo_shortening_imgw(data),
+    c("tmax_daily", "tmin_daily", "t5cm_min")
+  )
+})
+
+test_that("IMGW column labels are retained as attributes", {
+  data = data.frame(
+    TMAX = 1:2,
+    TMIN = 3:4,
+    TMNG = 5:6,
+    check.names = FALSE
+  )
+  meta = data.frame(
+    parameters = c("TMAX", "TMIN", "TMNG"),
+    label = c(
+      "Maksymalna dobowa temperatura powietrza [°C]",
+      "Minimalna dobowa temperatura powietrza [°C]",
+      "Minimalna dobowa temperatura powietrza przy gruncie [°C]"
+    ),
+    stringsAsFactors = FALSE
+  )
+
+  data = imgw_rename_params_to_labels(data, meta)
+  data = meteo_shortening_imgw(data)
+
+  expect_equal(
+    unname(vapply(data, attr, character(1), which = "label")),
+    meta$label
+  )
+})
diff --git a/vignettes/articles/pl.Rmd b/vignettes/articles/pl.Rmd
index 121543e1..87067c35 100644
--- a/vignettes/articles/pl.Rmd
+++ b/vignettes/articles/pl.Rmd
@@ -86,31 +86,31 @@ h = hydro_imgw(interval = "monthly", year = 2001:2010)
 head(h)
 ```
 
-Zmienna `MCWSKEX` reprezentuje etykietę ekstremum, gdzie "1" oznacza
+Zmienna `idex` reprezentuje etykietę ekstremum, gdzie "1" oznacza
 minimum, "2" oznacza średnią, a "3" maksimum. [^1]
 
 [^1]: Więcej informacji na ten temat można znaleźć w zestawie danych
-    `hydro_abbrev`.
+    `imgw_hydro_abbrev`.
 
 Analizy hydrologiczne często koncentrują się na jednej grupy zjawisk,
 np. związanych z przepływami maksymalnymi. W tym celu pozostaną w ramce
 danych tylko wartości przepływów maksymalnych oraz kolumny zawierające
 interesujące nas informacje, tj. identyfikator stacji - `id`, rok
-hydrologiczny (`MCROKH`), szerokość geograficzną `X` i długość
-geograficzną`Y`.
+hydrologiczny (`hyy`), nazwę stacji (`station`), szerokość geograficzną
+`X` i długość geograficzną `Y`.
 
 Następnie obliczymy średnią maksymalną wartość przepływu na stacjach w
 każdym roku za pomocą `dplyr::summarise()`, oraz rozdzielimy dane według
 roku używając `spread ()` aby uzyskać roczne średnie maksymalne
-przepływy (`MCPRZP`) w kolejnych kolumnach.
+przepływy (`Q`) w kolejnych kolumnach.
 
 ```{r filtering, eval=TRUE, echo=TRUE, message=FALSE, warning=FALSE, include=TRUE, paged.print=TRUE}
 h2 = h %>%
-  dplyr::filter(MCWSKEX == 3) %>%
-  dplyr::select(id, PSNZWP, X, Y, MCROKH, MCPRZP) %>%
-  dplyr::group_by(MCROKH, id, PSNZWP, X, Y) %>%
-  dplyr::summarise(srednie_roczne_Q = round(mean(MCPRZP, na.rm = TRUE), 1)) %>% 
-  spread(MCROKH, srednie_roczne_Q)
+  dplyr::filter(idex == 3) %>%
+  dplyr::select(id, station, X, Y, hyy, Q) %>%
+  dplyr::group_by(hyy, id, station, X, Y) %>%
+  dplyr::summarise(srednie_roczne_Q = round(mean(Q, na.rm = TRUE), 1)) %>%
+  spread(hyy, srednie_roczne_Q)
 ```
 
 ```{r filtering2, eval=TRUE, echo=TRUE, message=FALSE, warning=FALSE, include=TRUE, paged.print=TRUE}
diff --git a/vignettes/getstarted.Rmd b/vignettes/getstarted.Rmd
index e420fb69..ba5d7df0 100644
--- a/vignettes/getstarted.Rmd
+++ b/vignettes/getstarted.Rmd
@@ -42,6 +42,9 @@ default for `interval = "hourly"`) and **HTML scraping** (`source = "html"`, def
 SYNOP/CLIMATE/PRECIP stations available in the danepubliczne.imgw.pl collection.
 It is a wrapper for `meteo_monthly()`, `meteo_daily()`, `meteo_hourly()` and
 `meteo_imgw_datastore()` which gives access from monthly to even 10-min datasets.
+IMGW archive outputs use short English parameter names. The original IMGW names are
+available as `label` attributes on the corresponding columns, for example
+`attr(df$t2m_mean_mon, "label")`.
 
 - **meteo_noaa_hourly()** - Downloading hourly NCEI/NOAA Integrated Surface Hourly (ISH)
 meteorological data - some stations have > 100 years of observations.
@@ -59,6 +62,9 @@ data frames. For a full walkthrough see the
 
 - **hydro_imgw()** - Downloading daily and monthly hydrological data from stations available in the
 danepubliczne.imgw.pl collection. It is a wrapper for `hydro_monthly()` and `hydro_daily()`.
+IMGW archive outputs use short English parameter names. The original IMGW names are
+available as `label` attributes on the corresponding columns, for example
+`attr(h$Q, "label")`.
 
 - **hydro_imgw_datastore()** - Downloading hourly and sub-hourly hydrological data from the
 IMGW-PIB hydro telemetry stations.