fails when leaf_wetness = 1. The .diurnal_rh() returns a single value best I can tell from the original meteor package, but this appears to expect a vector?
If weather data for an entire year are supplied to SEIR() and a season is simulated starting later than the beginning of the time-series data, the values returned are only for the earliest dates of weather, not for the desired season.