[TI.by_speed] is not working with Pandas v.1.3.3 and above

[BiancaMorandi]

bw.TI.by_speed(DATA[['Spd80mN']], DATA[['Spd80mNStd']]) is raising an error with Pandas v.1.3.3 and above. This error is due to an update of pandas.DataFrame.agg from Pandas v.1.3.3. see https://pandas.pydata.org/pandas-docs/stable/whatsnew/v1.3.3.html

This update is generating an error when deriving pandas.DataFrame.agg(lambda x: np.percentile(x, q=percentile)). The error is generated for dist(var_series=ti['Turbulence_Intensity'], var_to_bin_against=ti['wspd'], bins=speed_bin_array, bin_labels=speed_bin_labels, aggregation_method=lambda x: np.percentile(x, q=percentile), return_data=True)[-1].rename("Rep_TI")

when called inside TI.by_speed

[BiancaMorandi]

additional error raised with numpy==1.23.1 and pandas==0.25.3 (or 1.4.3 ??)

[BiancaMorandi]

additional error raised for

[BiancaMorandi]

while looking at this issue a bug was found for the TI.by_speed and TI.by_sector functions when bin_labels are given as input.

Figure below shows that if speed_bin_labels is given as input then TI.by_speed is not showing the Mean_TI and Rep_TI curves in the plot.

for TI.by_sector if direction_bin_labels is given as input then an error is raised because plot cannot be generated.

[BiancaMorandi]

As part of the debugging of this issue it was observed that:

• turbulence intensity is derived for both TI.by_speed and TI.by_sector for speed > and not >= of 3 m/s. This bug is now fixed.
• the min_speed input given to TI.by_sector was actually not applied to the TI calculation. This bug is now fixed

[BiancaMorandi]

the TI.by_speed function can accept IEC_class inputs as DataFrame but if this doesn't contain three 4 columns (windspeed+ three classes) then TI.by_speed function raises and error. This bug is now fixed.