pyepwmorph 3.0.0__tar.gz → 3.2.0__tar.gz

{pyepwmorph-3.0.0 → pyepwmorph-3.2.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyepwmorph
-Version: 3.0.0
+Version: 3.2.0
 Summary: A python package to enable simple and easy gathering of climate model data and morphing of EPW files
 Project-URL: Homepage, https://github.com/justinfmccarty/pyepwmorph
 Project-URL: Issues, https://github.com/justinfmccarty/pyepwmorph/issues
@@ -51,6 +51,7 @@ A Python package for morphing EnergyPlus Weather (EPW) files with climate model
 
 - **CMIP6 data** fetched automatically from Google Cloud (Pangeo)
 - **Custom CSV data** from any climate model, including historical reconstructions
+- **CH2025 station scenarios** for Switzerland, shipped with the package and indexed by global warming level
 
 ## Installation
 
@@ -128,6 +129,29 @@ Custom CSVs should have a `date` column (parseable by pandas) and a column named
 
 The reference and target scenarios must cover **different years**, the same way the CMIP6 `historical` and `sspXXX` experiments do. The two series are concatenated before the baseline and target periods are sliced out, so overlapping years get averaged together and weaken the climate signal. Keep `baseline_range` inside the years the reference scenario covers.
 
+### Switzerland (CH2025 warming levels)
+
+For an EPW inside Switzerland, `data_source="ch2025"` morphs from MeteoSwiss CH2025 station scenarios. There is no target year: each pathway is a global warming level relative to 1991-2020, and the run is offline.
+
+```python
+results = workflow.morphing_workflow(
+    project_name="Zurich_GWL2",
+    epw_file="zurich.epw",
+    user_variables=["Temperature", "Humidity", "Wind", "Radiation", "Dew Point"],
+    user_pathways=["GWL 2.0"],
+    percentiles=[50],
+    output_directory="output/",
+    data_source="ch2025",
+)
+morphed = results["gwl2.0"]["50"]
+```
+
+The site is matched to the nearest station that has the requested variables (or, with `ch2025_full_coverage=True`, the nearest station that has every CH2025 variable), with an elevation penalty so a much higher station is not chosen just because it is close on the map. Supported variables are temperature, humidity, dew point, wind, and radiation (global, diffuse, and direct). Pressure and cloud cover are not in CH2025; requesting them raises an error. Sky cover is left at the EPW's baseline values when radiation is morphed, so longwave sky temperature in EnergyPlus does not follow the shortwave change.
+
+The signal is the change from the 1991-2020 reference climate to the chosen warming level. Each model chain's change is computed first, over the chains that reach that warming level, and the percentiles are taken of those changes. The EPW's period is read from the `Period of Record` in its header. If that is not exactly 1991-2020, a `UserWarning` recommends a TMY built from 1991-2020 data. If the header states no period, a second warning says the period could not be detected, and the mismatch warning still fires if the data rows' years fall outside 1991-2020. The signal is applied to the file as it stands either way. `tools.io.write_period_of_record` adds the statement to files built without one. The same caveats are listed in `MorphConfig.ch2025_notes` and the matched station in `MorphConfig.ch2025_station`.
+
+CH2025 data: MeteoSwiss & ETH Zurich (2025), Climate CH2025 - Daily Datasets, CC-BY 4.0, https://doi.org/10.18751/climate/scenarios/ch2025/data/1.0/
+
 ## Climate scenarios
 
 | Scenario              | SSP    | Description                             | Expected warming |
@@ -143,15 +167,16 @@ The reference and target scenarios must cover **different years**, the same way
 - **Humidity** -- relative humidity, stretched in specific humidity space
 - **Pressure** -- atmospheric pressure (shift)
 - **Wind** -- wind speed (stretch)
-- **Clouds and Radiation** -- global/diffuse/direct radiation and sky cover
+- **Radiation** -- global/diffuse/direct radiation, without changing sky cover (all sources; the only radiation option for CH2025)
+- **Clouds and Radiation** -- global/diffuse/direct radiation and sky cover (CMIP6 and custom data)
 - **Dew Point** -- recalculated from morphed temperature and humidity
 
 ### Variable dependencies
 
 Some variables cannot be morphed on their own:
 
-- **Humidity** requires Temperature and Pressure
-- **Dew Point** requires Temperature, Humidity, and Pressure
+- **Humidity** requires Temperature and Pressure, except with CH2025, where relative humidity is stretched directly
+- **Dew Point** requires Temperature, Humidity, and Pressure (Temperature and Humidity only with CH2025)
 
 Dependencies are added automatically and are **written to the output file**. Asking for `Dew Point` alone therefore returns an EPW with morphed pressure, temperature, relative humidity, and dew point, which keeps the file internally consistent. `MorphConfig.resolved_variables` shows exactly what will be written, in the order it is computed.
 
@@ -220,7 +245,7 @@ See [CHANGELOG.md](CHANGELOG.md) for the full history. The most recent release c
 
 - Python >= 3.9
 - pandas >= 2.2
-- Internet connection (for CMIP6 data download; the custom CSV workflow runs offline)
+- Internet connection (for CMIP6 data download; the custom CSV and CH2025 workflows run offline)
 
 ## License
 

{pyepwmorph-3.0.0 → pyepwmorph-3.2.0}/README.md

@@ -12,6 +12,7 @@ A Python package for morphing EnergyPlus Weather (EPW) files with climate model
 
 - **CMIP6 data** fetched automatically from Google Cloud (Pangeo)
 - **Custom CSV data** from any climate model, including historical reconstructions
+- **CH2025 station scenarios** for Switzerland, shipped with the package and indexed by global warming level
 
 ## Installation
 
@@ -89,6 +90,29 @@ Custom CSVs should have a `date` column (parseable by pandas) and a column named
 
 The reference and target scenarios must cover **different years**, the same way the CMIP6 `historical` and `sspXXX` experiments do. The two series are concatenated before the baseline and target periods are sliced out, so overlapping years get averaged together and weaken the climate signal. Keep `baseline_range` inside the years the reference scenario covers.
 
+### Switzerland (CH2025 warming levels)
+
+For an EPW inside Switzerland, `data_source="ch2025"` morphs from MeteoSwiss CH2025 station scenarios. There is no target year: each pathway is a global warming level relative to 1991-2020, and the run is offline.
+
+```python
+results = workflow.morphing_workflow(
+    project_name="Zurich_GWL2",
+    epw_file="zurich.epw",
+    user_variables=["Temperature", "Humidity", "Wind", "Radiation", "Dew Point"],
+    user_pathways=["GWL 2.0"],
+    percentiles=[50],
+    output_directory="output/",
+    data_source="ch2025",
+)
+morphed = results["gwl2.0"]["50"]
+```
+
+The site is matched to the nearest station that has the requested variables (or, with `ch2025_full_coverage=True`, the nearest station that has every CH2025 variable), with an elevation penalty so a much higher station is not chosen just because it is close on the map. Supported variables are temperature, humidity, dew point, wind, and radiation (global, diffuse, and direct). Pressure and cloud cover are not in CH2025; requesting them raises an error. Sky cover is left at the EPW's baseline values when radiation is morphed, so longwave sky temperature in EnergyPlus does not follow the shortwave change.
+
+The signal is the change from the 1991-2020 reference climate to the chosen warming level. Each model chain's change is computed first, over the chains that reach that warming level, and the percentiles are taken of those changes. The EPW's period is read from the `Period of Record` in its header. If that is not exactly 1991-2020, a `UserWarning` recommends a TMY built from 1991-2020 data. If the header states no period, a second warning says the period could not be detected, and the mismatch warning still fires if the data rows' years fall outside 1991-2020. The signal is applied to the file as it stands either way. `tools.io.write_period_of_record` adds the statement to files built without one. The same caveats are listed in `MorphConfig.ch2025_notes` and the matched station in `MorphConfig.ch2025_station`.
+
+CH2025 data: MeteoSwiss & ETH Zurich (2025), Climate CH2025 - Daily Datasets, CC-BY 4.0, https://doi.org/10.18751/climate/scenarios/ch2025/data/1.0/
+
 ## Climate scenarios
 
 | Scenario              | SSP    | Description                             | Expected warming |
@@ -104,15 +128,16 @@ The reference and target scenarios must cover **different years**, the same way
 - **Humidity** -- relative humidity, stretched in specific humidity space
 - **Pressure** -- atmospheric pressure (shift)
 - **Wind** -- wind speed (stretch)
-- **Clouds and Radiation** -- global/diffuse/direct radiation and sky cover
+- **Radiation** -- global/diffuse/direct radiation, without changing sky cover (all sources; the only radiation option for CH2025)
+- **Clouds and Radiation** -- global/diffuse/direct radiation and sky cover (CMIP6 and custom data)
 - **Dew Point** -- recalculated from morphed temperature and humidity
 
 ### Variable dependencies
 
 Some variables cannot be morphed on their own:
 
-- **Humidity** requires Temperature and Pressure
-- **Dew Point** requires Temperature, Humidity, and Pressure
+- **Humidity** requires Temperature and Pressure, except with CH2025, where relative humidity is stretched directly
+- **Dew Point** requires Temperature, Humidity, and Pressure (Temperature and Humidity only with CH2025)
 
 Dependencies are added automatically and are **written to the output file**. Asking for `Dew Point` alone therefore returns an EPW with morphed pressure, temperature, relative humidity, and dew point, which keeps the file internally consistent. `MorphConfig.resolved_variables` shows exactly what will be written, in the order it is computed.
 
@@ -181,7 +206,7 @@ See [CHANGELOG.md](CHANGELOG.md) for the full history. The most recent release c
 
 - Python >= 3.9
 - pandas >= 2.2
-- Internet connection (for CMIP6 data download; the custom CSV workflow runs offline)
+- Internet connection (for CMIP6 data download; the custom CSV and CH2025 workflows run offline)
 
 ## License
 

pyepwmorph-3.2.0/pyepwmorph/data/__init__.py

@@ -0,0 +1 @@
+"""Shipped climate datasets."""

{pyepwmorph-3.0.0 → pyepwmorph-3.2.0}/pyepwmorph/models/assemble.py

@@ -109,3 +109,27 @@ def calc_model_climatologies(baseline_range, future_range, baseline_data, future
 
 
     return baseline_means, future_means
+
+
+def calc_gwl_climatologies(baseline_data, future_data, variable):
+    """Return CH2025 monthly climatologies unchanged.
+
+    A warming-level state is already a stationary 30-year sample, so there
+    is no year range to slice the way ``calc_model_climatologies`` does for
+    a transient CMIP6 run.
+
+    Parameters
+    ----------
+    baseline_data : pd.Series
+        Twelve monthly values for the reference state.
+    future_data : pd.Series
+        Twelve monthly values for the warming-level state.
+    variable : str
+        Name used to rename both series.
+
+    Returns
+    -------
+    tuple
+        ``(baseline_data, future_data)``, each renamed to *variable*.
+    """
+    return baseline_data.rename(variable), future_data.rename(variable)

pyepwmorph-3.2.0/pyepwmorph/models/ch2025.py

@@ -0,0 +1,301 @@
+"""Offline access to MeteoSwiss CH2025 station scenarios.
+
+The monthly climatologies in ``pyepwmorph/data`` are reduced from the
+DAILY-LOCAL product. Each row is one model chain's mean for a calendar
+month, over the 30 synthetic climate years of a reference state or a
+global warming level.
+
+CH2025 data © MeteoSwiss & ETH Zurich, licensed CC-BY 4.0.
+https://doi.org/10.18751/climate/scenarios/ch2025/data/1.0/
+"""
+
+import logging
+from functools import lru_cache
+from importlib.resources import files
+
+import numpy as np
+import pandas as pd
+
+logger = logging.getLogger(__name__)
+
+#: (west, south, east, north) of the CH2025 gridded domain, in degrees.
+CH2025_BBOX = (5.96, 45.82, 10.49, 47.81)
+
+CH2025_REFERENCE = "ref91-20"
+CH2025_BASELINE_RANGE = (1991, 2020)
+CH2025_STATES = ("gwl1.5", "gwl2.0", "gwl2.5", "gwl3.0")
+CH2025_VARIABLES = frozenset({"tas", "tasmax", "tasmin", "hurs", "rsds", "sfcWind"})
+
+#: Variables whose change is a difference (future - reference). The others
+#: (``hurs``, ``rsds``, ``sfcWind``) change as a ratio, matching how the
+#: morphing procedures apply them.
+CH2025_ADDITIVE_VARIABLES = frozenset({"tas", "tasmax", "tasmin"})
+
+#: Short attribution written into morphed EPW comments. No commas: EPW
+#: header fields are comma delimited.
+CH2025_CITATION = (
+    "CH2025 (c) MeteoSwiss and ETH Zurich CC-BY 4.0 "
+    "doi:10.18751/climate/scenarios/ch2025/data/1.0"
+)
+
+#: Kilometres of horizontal distance treated as equal to 100 m of elevation,
+#: so a station far below or above the site loses to a nearer-altitude one.
+DEFAULT_ELEVATION_WEIGHT_KM_PER_100M = 10.0
+WARN_DISTANCE_KM = 25.0
+WARN_ELEVATION_M = 300.0
+
+
+def _read_parquet(name: str) -> pd.DataFrame:
+    resource = files("pyepwmorph.data").joinpath(name)
+    with resource.open("rb") as handle:
+        return pd.read_parquet(handle)
+
+
+@lru_cache(maxsize=1)
+def load_table() -> pd.DataFrame:
+    """Return the shipped monthly climatology table.
+
+    Columns are ``station_id``, ``variable``, ``state``, ``chain``,
+    ``month``, and ``value``. The frame is cached; do not mutate it.
+    """
+    return _read_parquet("ch2025_monthly.parquet")
+
+
+@lru_cache(maxsize=1)
+def _load_stations() -> pd.DataFrame:
+    return _read_parquet("ch2025_stations.parquet")
+
+
+def available_stations(variables=None) -> pd.DataFrame:
+    """Return stations, optionally limited to those carrying every variable.
+
+    Parameters
+    ----------
+    variables : iterable of str or None
+        Model variable ids that a station must provide. ``None`` returns
+        every station in the shipped table.
+    """
+    stations = _load_stations().copy()
+    if variables is None:
+        return stations.reset_index(drop=True)
+    needed = set(variables)
+    mask = stations["variables"].map(lambda present: needed <= set(present))
+    return stations.loc[mask].reset_index(drop=True)
+
+
+def reference_key(state: str) -> str:
+    """Key under which the reference paired with *state* is stored.
+
+    Each warming level is compared with the reference over the model chains
+    the two share, so the reference climatology differs per warming level.
+    """
+    return f"{CH2025_REFERENCE}|{state}"
+
+
+def baseline_matches(baseline_range) -> bool:
+    """Return whether an EPW baseline is exactly the CH2025 reference period."""
+    return tuple(int(year) for year in baseline_range) == CH2025_BASELINE_RANGE
+
+
+def in_switzerland(latitude: float, longitude: float) -> bool:
+    """Return whether a point lies inside the CH2025 domain."""
+    west, south, east, north = CH2025_BBOX
+    return south <= latitude <= north and west <= longitude <= east
+
+
+def _haversine_km(latitude, longitude, latitudes, longitudes) -> np.ndarray:
+    radius = 6371.0
+    lat1 = np.radians(latitude)
+    lat2 = np.radians(np.asarray(latitudes, dtype=float))
+    dlat = lat2 - lat1
+    dlon = np.radians(np.asarray(longitudes, dtype=float) - longitude)
+    a = np.sin(dlat / 2.0) ** 2 + np.cos(lat1) * np.cos(lat2) * np.sin(dlon / 2.0) ** 2
+    return 2.0 * radius * np.arcsin(np.sqrt(np.clip(a, 0.0, 1.0)))
+
+
+def nearest_station(
+    latitude: float,
+    longitude: float,
+    elevation: float,
+    variables=None,
+    elevation_weight_km_per_100m: float = DEFAULT_ELEVATION_WEIGHT_KM_PER_100M,
+) -> pd.Series:
+    """Pick the station that best matches a site.
+
+    The score is horizontal distance plus an elevation penalty, so a
+    station hundreds of metres higher loses to one farther away at a
+    similar altitude. Only stations that carry every requested variable
+    are eligible.
+
+    Returns
+    -------
+    pd.Series
+        The station row, with ``distance_km`` and ``elevation_difference_m``.
+    """
+    candidates = available_stations(variables)
+    if candidates.empty:
+        needed = ", ".join(sorted(variables or []))
+        raise ValueError(f"No CH2025 station provides {needed or 'any variables'}")
+
+    distance = _haversine_km(
+        latitude, longitude, candidates["latitude"], candidates["longitude"],
+    )
+    elevation_difference = candidates["elevation"].to_numpy(dtype=float) - float(elevation)
+    penalty = np.abs(elevation_difference) / 100.0 * elevation_weight_km_per_100m
+    position = int(np.argmin(distance + penalty))
+    chosen = candidates.iloc[position].copy()
+    chosen["distance_km"] = float(distance[position])
+    chosen["elevation_difference_m"] = float(elevation_difference[position])
+    chosen["weak_match"] = bool(
+        chosen["distance_km"] > WARN_DISTANCE_KM
+        or abs(chosen["elevation_difference_m"]) > WARN_ELEVATION_M
+    )
+
+    logger.info(
+        "CH2025 station %s (%s): %.1f km away, %+.0f m elevation",
+        chosen["station_id"],
+        chosen["name"],
+        chosen["distance_km"],
+        chosen["elevation_difference_m"],
+    )
+    if chosen["weak_match"]:
+        logger.warning(
+            "CH2025 station %s is a weak match (%.1f km, %+.0f m). "
+            "The morphed file uses that station's climate signal.",
+            chosen["station_id"],
+            chosen["distance_km"],
+            chosen["elevation_difference_m"],
+        )
+    return chosen
+
+
+def _check_variable(variable):
+    if variable not in CH2025_VARIABLES:
+        raise ValueError(
+            f"Variable '{variable}' is not in the CH2025 station table. "
+            f"Available: {sorted(CH2025_VARIABLES)}"
+        )
+
+
+def _chains_by_month(station_id, variable, state) -> pd.DataFrame:
+    """Return a month x chain table for one station, variable, and state."""
+    table = load_table()
+    subset = table[
+        (table["station_id"] == station_id)
+        & (table["variable"] == variable)
+        & (table["state"] == state)
+    ]
+    if subset.empty:
+        raise ValueError(
+            f"No CH2025 data for station '{station_id}', variable '{variable}', state '{state}'"
+        )
+    return subset.pivot(index="month", columns="chain", values="value").sort_index()
+
+
+def build_ch2025_change_ensemble(percentiles, variable, station_id, state):
+    """Build reference and warming-level climatologies from paired chain changes.
+
+    Each model chain's change from the reference to *state* is computed
+    first, over the chains present in both, and the percentiles are taken
+    of those changes. Taking percentiles of each state separately and
+    differencing them would mix chains, and at the tails can even flip the
+    sign of the change.
+
+    The reference returned is the median of the paired chains, the same for
+    every percentile column. The warming-level frame is that reference plus
+    the percentile change (temperatures) or times the percentile ratio
+    (humidity, radiation, wind), so the morphing procedures recover exactly
+    the percentile change.
+
+    Parameters
+    ----------
+    percentiles : list
+        Ensemble percentiles. Column labels match these values.
+    variable : str
+        One of ``CH2025_VARIABLES``.
+    station_id : str
+        Lower-case station abbreviation, for example ``"sma"``.
+    state : str
+        A warming level such as ``"gwl2.0"``.
+
+    Returns
+    -------
+    tuple[pd.DataFrame, pd.DataFrame]
+        ``(reference, future)``, each with twelve rows (months 1-12) and one
+        column per percentile. ``attrs["n_chains"]`` on both records how many
+        paired chains contributed.
+    """
+    _check_variable(variable)
+    reference = _chains_by_month(station_id, variable, CH2025_REFERENCE)
+    future = _chains_by_month(station_id, variable, state)
+    common = reference.columns.intersection(future.columns)
+    if len(common) == 0:
+        raise ValueError(
+            f"No CH2025 model chain provides both {CH2025_REFERENCE} and {state} "
+            f"for station '{station_id}', variable '{variable}'"
+        )
+    reference = reference[common]
+    future = future[common]
+    baseline = reference.median(axis=1)
+
+    if variable in CH2025_ADDITIVE_VARIABLES:
+        change = future - reference
+    else:
+        ref_values = reference.to_numpy(dtype=float)
+        change = pd.DataFrame(
+            np.divide(
+                future.to_numpy(dtype=float), ref_values,
+                out=np.ones_like(ref_values), where=ref_values != 0,
+            ),
+            index=reference.index, columns=reference.columns,
+        )
+
+    reference_frame = pd.DataFrame({percentile: baseline for percentile in percentiles})
+    future_columns = {}
+    for percentile in percentiles:
+        quantile = change.quantile(int(percentile) / 100.0, axis=1)
+        if variable in CH2025_ADDITIVE_VARIABLES:
+            future_columns[percentile] = baseline + quantile
+        else:
+            future_columns[percentile] = baseline * quantile
+    future_frame = pd.DataFrame(future_columns)
+
+    for frame in (reference_frame, future_frame):
+        frame.index.name = "month"
+        frame.attrs["n_chains"] = int(len(common))
+    return reference_frame, future_frame
+
+
+def build_ch2025_ensemble(percentiles, variable, station_id, state) -> pd.DataFrame:
+    """Reduce one station, variable, and state to percentile climatologies.
+
+    Parameters
+    ----------
+    percentiles : list
+        Ensemble percentiles. Column labels match these values, as in
+        ``assemble.build_cmip6_ensemble``.
+    variable : str
+        One of ``CH2025_VARIABLES``.
+    station_id : str
+        Lower-case station abbreviation, for example ``"sma"``.
+    state : str
+        ``"ref91-20"`` or a warming level such as ``"gwl2.0"``.
+
+    Returns
+    -------
+    pd.DataFrame
+        Twelve rows (months 1-12) and one column per percentile.
+        ``attrs["n_chains"]`` records how many model chains contributed.
+        Membership differs by variable and warming level.
+
+    Morphing does not use this: differencing two separately reduced states
+    mixes model chains. See ``build_ch2025_change_ensemble``.
+    """
+    _check_variable(variable)
+    wide = _chains_by_month(station_id, variable, state)
+    result = pd.DataFrame(
+        {percentile: wide.quantile(int(percentile) / 100.0, axis=1) for percentile in percentiles}
+    )
+    result.index.name = "month"
+    result.attrs["n_chains"] = int(wide.shape[1])
+    return result

{pyepwmorph-3.0.0 → pyepwmorph-3.2.0}/pyepwmorph/morph/procedures.py

@@ -7,6 +7,8 @@ Container for all of the individual morphing calculations which can be traced ba
         55 514–524 ISSN 0960-1481 URL https://www.sciencedirect.com/science/article/pii/S0960148113000232
 
 """
+import calendar
+
 import numpy as np
 import pandas as pd
 
@@ -33,10 +35,17 @@ def _index_of(*candidates):
 
 
 def _year_of(index):
-    """Return the year of a DatetimeIndex, or the default solar year."""
+    """Return a non-leap year to build solar geometry for the data in *index*.
+
+    EPW files always hold 8760 hours with no 29 February, so their calendar
+    is a non-leap year's even when the first data row says 2016. Building
+    daily series for a leap year would add an empty 29 February and shift
+    every later day by one, so a leap year is replaced by the year before.
+    """
     if index is None or len(index) == 0:
         return morph_solar_utils.DEFAULT_SOLAR_YEAR
-    return int(index[0].year)
+    year = int(index[0].year)
+    return year - 1 if calendar.isleap(year) else year
 
 
 def _as_series(values, index, name, decimals=2):
@@ -221,6 +230,34 @@ def morph_relhum(present_relhum, present_psl, present_dbt, future_psl, future_db
     return _as_series(np.clip(morphed_relhum, 1, 100), index, "relhum_percent")
 
 
+def morph_relhum_direct(present_relhum, future_hurs, baseline_hurs):
+    """Stretch relative humidity by the modelled relative-humidity ratio.
+
+    CH2025 provides ``hurs`` directly, so the specific-humidity round trip
+    used by ``morph_relhum`` is not needed.
+
+    Parameters
+    ----------
+    present_relhum : pd.Series
+        Hourly present-day relative humidity in percent.
+    future_hurs : pd.Series
+        Twelve monthly future relative-humidity values.
+    baseline_hurs : pd.Series
+        Twelve monthly baseline relative-humidity values.
+
+    Returns
+    -------
+    pd.Series
+        Morphed relative humidity, clipped to ``[1, 100]``.
+    """
+    index = present_relhum.index
+    ratio = morph_utils.relative_delta(
+        morph_utils.as_array(future_hurs, 12), morph_utils.as_array(baseline_hurs, 12),
+    )
+    morphed = stretch(present_relhum.to_numpy(dtype=float), morph_utils.month_factors(index, ratio))
+    return _as_series(np.clip(morphed, 1, 100), index, "relhum_percent")
+
+
 def morph_psl(present_psl, future_psl, baseline_psl):
     """
     Pressure at sea level morph requires a shift
@@ -327,6 +364,31 @@ def morph_wspd(present_wspd, future_vas, baseline_vas, future_uas, baseline_uas)
     return _as_series(np.clip(morphed_wspd, 0, None), index, "windspd_ms")
 
 
+def morph_wspd_direct(present_wspd, future_sfcwind, baseline_sfcwind):
+    """Stretch wind speed by a scalar surface-wind ratio.
+
+    Parameters
+    ----------
+    present_wspd : pd.Series
+        Hourly present-day wind speed in m/s.
+    future_sfcwind : pd.Series
+        Twelve monthly future surface wind speeds.
+    baseline_sfcwind : pd.Series
+        Twelve monthly baseline surface wind speeds.
+
+    Returns
+    -------
+    pd.Series
+        Morphed wind speed, clipped at zero.
+    """
+    index = present_wspd.index
+    ratio = morph_utils.relative_delta(
+        morph_utils.as_array(future_sfcwind, 12), morph_utils.as_array(baseline_sfcwind, 12),
+    )
+    morphed = stretch(present_wspd.to_numpy(dtype=float), morph_utils.month_factors(index, ratio))
+    return _as_series(np.clip(morphed, 0, None), index, "windspd_ms")
+
+
 def morph_glohor(present_glohor, future_glohor, baseline_glohor):
     """
     Global horizontal radiation morph requires a stretch