From c077271257a08d607aff2dd53a1522f78557304d Mon Sep 17 00:00:00 2001
From: Brian Groenke <brian.groenke@awi.de>
Date: Wed, 1 Jan 2025 19:20:14 +0100
Subject: [PATCH] Add long overdue param annotations to SEB

---
 src/Physics/Surface/SEB/seb.jl | 15 +++++++++++++--
 1 file changed, 13 insertions(+), 2 deletions(-)

diff --git a/src/Physics/Surface/SEB/seb.jl b/src/Physics/Surface/SEB/seb.jl
index b5245cdb..0db86369 100644
--- a/src/Physics/Surface/SEB/seb.jl
+++ b/src/Physics/Surface/SEB/seb.jl
@@ -48,8 +48,19 @@ Utils.@properties SEBParams(
     # surface propertiess; includes soil and snow by default;
     # this can be easily extended by the user to include other layers by adding more
     # fields to `SEBParams` and then implementing `surfaceproperties` for those layers.
-    soil = SurfaceProperties(),
-    snow = SurfaceProperties(α=0.8, ϵ=0.99, z₀=5e-4),
+    soil = SurfaceProperties(
+        α = param(0.2, desc="Initial snow-free surface albedo"),
+        ϵ = param(0.97, desc="Initial snow-free surface emissivity"),
+        zâ‚€ = param(1e-3, units=u"m", desc="Initial snow-free surface roughness length"),
+        râ‚› = param(50.0, units=u"s/m", desc="Initial snow-free surface resistance against evapotranspiration and sublimation"),
+    ),
+    
+    snow = SurfaceProperties(
+        α = param(0.8, desc="Initial snow-covered surface albedo"),
+        ϵ = param(0.99, desc="Initial snow-covered surface emissivity"),
+        zâ‚€ = param(5e-4, desc="Initial snow-covered surface roughness length"),
+        râ‚› = param(50.0, units=u"s/m", desc="Initial snow-covered surface resistance against evapotranspiration and sublimation"),
+    ),
 
     # "natural" constant
     σ = 5.6704e-8u"J/(s*m^2*K^4)",   # Stefan-Boltzmann constant
-- 
GitLab