-
Notifications
You must be signed in to change notification settings - Fork 19
Description
Hi there,
I'm trying to follow along with the advanced vignette and am not getting the same result. I've been trying to track down the issue. Not sure I found my main problem but did find a bug.
in calcH the calculation of u200.v is looking for values that do not exist in the provided WeatherStation data. Therefore, it is reverting to the default value of 4 m/s for u200.v. Namely WeatherStation$wind and WeatherStation$height. I've found what I think are the right values WeatherStation$daily$wind_mean and WeatherStation$location$height and can succesfully calculate a u200.v value with them:
Z.om.ws <- 0.03 #from function defaults
u.ws <- WeatherStation$daily$wind_mean * 0.41/log(WeatherStation$location$height/Z.om.ws)
u200.v <- u.ws/0.41 * log(200/Z.om.ws) #returns 6.29
Thanks for all your work on this!