Feature/wind profiles

docs/source/api/abiotic/wind.md

@@ -0,0 +1,24 @@
+---
+jupytext:
+  cell_metadata_filter: -all
+  formats: md:myst
+  main_language: python
+  text_representation:
+    extension: .md
+    format_name: myst
+    format_version: 0.13
+    jupytext_version: 1.13.8
+kernelspec:
+  display_name: vr_python3
+  language: python
+  name: vr_python3
+---
+
+#  API for the {mod}`~virtual_rainforest.models.abiotic.wind` module
+
+```{eval-rst}
+.. automodule:: virtual_rainforest.models.abiotic.wind
+    :autosummary:
+    :members:
+    :special-members: __init__
+```

docs/source/index.md

@@ -108,6 +108,7 @@ team.
   Abiotic Mechanistic Model <api/abiotic/abiotic_model.md>
   Abiotic Mechanistic Constants <api/abiotic/abiotic_constants.md>
   Abiotic Mechanistic Tools <api/abiotic/abiotic_tools.md>
+  Abiotic Mechanistic Wind <api/abiotic/wind.md>
   Hydrology Overview <api/hydrology.md>
   Hydrology Model <api/hydrology/hydrology_model.md>
   Hydrology Above-ground <api/hydrology/above_ground.md>

docs/source/refs.bib

@@ -1,4 +1,14 @@
 
+@article{yasuda_turbulent_1988,
+	title = {Turbulent diffusivity and diurnal variations in the atmospheric boundary layer.},
+	volume = {43},
+	doi = {https://doi.org/10.1007/BF00128403},
+	journal = {Boundary-layer meteorology},
+	author = {Yasuda, N.},
+	year = {1988},
+	pages = {209--221},
+}
+
 @article{henderson-sellers_new_1984,
 	title = {A new formula for latent heat of vaporization of water as a function of temperature},
 	volume = {110},

tests/conftest.py

@@ -319,6 +319,9 @@ def dummy_climate_data(layer_roles_fixture):
         np.full((3, 3), 30),
         dims=["cell_id", "time_index"],
     )
+    data["wind_speed_ref"] = DataArray(
+        [[0, 5, 10], [0, 5, 10], [0, 5, 10]], dims=["time_index", "cell_id"]
+    )
     data["mean_annual_temperature"] = DataArray(
         np.full((3), 20),
         dims=["cell_id"],
@@ -339,6 +342,18 @@ def dummy_climate_data(layer_roles_fixture):
         np.full((3, 3), 400),
         dims=["cell_id", "time_index"],
     )
+    pressure = np.repeat(a=[96.0, np.nan], repeats=[13, 2])
+    data["atmospheric_pressure"] = DataArray(
+        np.broadcast_to(pressure, (3, 15)).T,
+        dims=["layers", "cell_id"],
+        coords={
+            "layers": np.arange(15),
+            "layer_roles": ("layers", layer_roles_fixture),
+            "cell_id": data.grid.cell_id,
+        },
+        name="atmospheric_pressure",
+    )
+
     evapotranspiration = np.repeat(a=[np.nan, 20.0, np.nan], repeats=[1, 3, 11])
     data["evapotranspiration"] = DataArray(
         np.broadcast_to(evapotranspiration, (3, 15)).T,
@@ -472,5 +487,8 @@ def dummy_climate_data(layer_roles_fixture):
         np.full((2, 3), 450),
         dims=("groundwater_layers", "cell_id"),
     )
+    data["canopy_height"] = DataArray([32, 32, 32], dims=["cell_id"])
+    data["sensible_heat_flux_topofcanopy"] = DataArray([100, 50, 10], dims=["cell_id"])
+    data["friction_velocity"] = DataArray([12, 5, 2], dims=["cell_id"])
 
     return data