Potential Layer Selection Error in Storm Motion Calculation (`wrf_relhl.f90`)

Issue Description:​

In wrf_relhl.f90, the storm motion vector is estimated as 75% of the mean wind speed and 30° to the right of the mean wind direction (lines 103-113). This method is similar to Maddox et al. (1976), which computes the mean wind using surface, 850-, 700-, 500-, 300-, and 200-hPa levels. However, a potential layer selection issue exists in the mean wind calculation between 3000–10,000 m AGL:

1. The code identifies k3(layer at/above 3000 m AGL) and k10(layer at/above 10,000 m AGL) (lines 70–81)
2. The current loop DO k = k3, k10, -1(lines 87–92) likely includes data abovethe target heights, not strictly between3000–10,000 m AGL.

Proposed Fix:​​

Change the loop to DO k = k3+1, k10+1, -1 to ensure averaging only layers within 3000–10,000 m AGL

​Questions:​​

1. Why is the 3000–10,000 m AGL layer used? Is the current inclusion of heights outside this range intentional?
2. Is the proposed fix appropriate to align with the physical intent?

Thanks!

Inspect wrf_relhl.f90 around lines 70–113, starting with how k3 and k10 are selected and how the DO k = k3, k10, -1 loop averages wind layers. Check the Maddox et al. reference and the physical intent of the 3000–10,000 m AGL range. Done means determining whether the indexing is intentional and documenting or implementing the correct bounds.