Possible issue with meshing blade in c2d branch
I found this issue on 5/9/2022 and I think I have almost forget...
As I remember, there are two different methods (orders) to generate the blade mesh. And the default one is wrong while the non-default one is correct.
I guess, in: https://gitlab.windenergy.dtu.dk/frza/PGL/-/blob/c2d/PGL/components/loftedblade.py#L400
twist_corr = self._twist_correction(R, Vtheta_unit[i])
should be changed to:
twist_corr = np.arctan2(R[0, 1], R[0, 0])
I should check this later and well document things down.