Issue by combining ZongGaussianDeficit with EqGridRotorAvg

Hello,

By calling the code in copy, I raise an error and I can’t understand why I get this. Issue_RotorModel.py

I work on a wind farm with a Weibull Distribution per wind direction. I would like, at the end, to measure the losses due to wake losses. I use some gaussian wake deficit model (ZongGaussianDeficit, BastankhahGaussianDeficit or NiayifarGaussianDeficit) and the rotor average model EqGridRotorAvg(9).

Nevertheless, when simulating the situation (line 66), an error occurs : ValueError: operands could not be broadcast together with remapped shapes [original->remapped]: (1,3,360,23,1) and requested shape (1,3,360,1,69). I imagine this error comes from the combination of Gaussian Wake deficit model with EqGridRotorAvg(9) because the error disapears when : -I change the rotor average to RotorCenter() -I change the wake deficit model to NOJDeficit() -I call only 1 wind speed in the simulation (for instance ws = 10 instead of ws = None)

Do you see any mistakes on my code or is it an error you have already highlighted? Thank you for your help and for the development of PyWake