hotfix: near_sfc_height = 10 for wind power law features.

NREL · Sep 14, 2023 · 611a3b3 · 611a3b3
1 parent 90f7921
commit 611a3b3
Showing 1 changed file with 6 additions and 4 deletions.
diff --git a/sup3r/preprocessing/feature_handling.py b/sup3r/preprocessing/feature_handling.py
@@ -736,9 +736,11 @@ class UWindPowerLaw(DerivedFeature):
     method. Uses power law extrapolation to get values above surface
 
     https://csl.noaa.gov/projects/lamar/windshearformula.html
+    https://www.tandfonline.com/doi/epdf/10.1080/00022470.1977.10470503
     """
 
     ALPHA = 0.2
+    NEAR_SFC_HEIGHT = 10
 
     @classmethod
     def inputs(cls, feature):
@@ -774,18 +776,19 @@ def compute(cls, data, height):
             Derived feature array
 
         """
-        near_surface_height = 1
-        return data['uas'] * (float(height) / near_surface_height)**cls.ALPHA
+        return data['uas'] * (float(height) / cls.NEAR_SFC_HEIGHT)**cls.ALPHA
 
 
 class VWindPowerLaw(DerivedFeature):
     """V wind component feature class with needed inputs method and compute
     method. Uses power law extrapolation to get values above surface
 
     https://csl.noaa.gov/projects/lamar/windshearformula.html
+    https://www.tandfonline.com/doi/epdf/10.1080/00022470.1977.10470503
     """
 
     ALPHA = 0.2
+    NEAR_SFC_HEIGHT = 10
 
     @classmethod
     def inputs(cls, feature):
@@ -821,8 +824,7 @@ def compute(cls, data, height):
             Derived feature array
 
         """
-        near_surface_height = 1
-        return data['vas'] * (float(height) / near_surface_height)**cls.ALPHA
+        return data['vas'] * (float(height) / cls.NEAR_SFC_HEIGHT)**cls.ALPHA
 
 
 class UWind(DerivedFeature):