Wake GIF

[rafmudaf]

Here's a simple script to create a gif from FLORIS calculations. I used the gif matplotlib animation library. This is similar to the script used to produce the image on the FLORIS docs page.

import gif
import matplotlib as mpl
import matplotlib.pyplot as plt
import numpy as np

import floris.tools.visualization as wakeviz
from floris.tools import FlorisInterface


mpl.rcParams['figure.figsize'] = [20, 10]


fi = FlorisInterface("gif.yaml")
fi.calculate_wake()
iea15mw_hub_height = fi.floris.farm.hub_heights[0,0,0]


@gif.frame
def plot_horizontal(positive, negative, diff, yaw_angles):
    axes = plt.subplot(3,1,1)
    wakeviz.visualize_cut_plane(positive, ax=axes, color_bar=True)

    axes = plt.subplot(3,1,2)
    wakeviz.visualize_cut_plane(negative, ax=axes, color_bar=True)

    axes = plt.subplot(3,1,3)
    abs_diff = np.abs(diff.df.u)
    min_speed = -1 * np.max(abs_diff)
    max_speed = np.max(abs_diff)
    wakeviz.visualize_cut_plane(diff, ax=axes, color_bar=True, cmap="PuOr", minSpeed=min_speed, maxSpeed=max_speed)


@gif.frame
def plot_yaw(cut_plane, fi):
    axes = plt.subplot(1,1,1)
    wakeviz.visualize_cut_plane(cut_plane, ax=axes, cmap="twilight")
    wakeviz.plot_turbines_with_fi(fi, axes)


fi.reinitialize(layout_x=[0.0, 750.0, 1500.0], layout_y=[0.0, 150.0, -150.0])

frames = []

N = 50
yaw = np.arange(-25, 25, (25 + 25) / N )
yaw = np.append(yaw, np.arange(25, -25, -1 * (25 + 25) / N ) )

for i in range(2 * N):

    yaw_angles = np.zeros((1,1,3))
    yaw_angles[:,:,0] = yaw[i]
    fi.floris.farm.yaw_angles = yaw_angles
    negative = fi.calculate_horizontal_plane(
        height=iea15mw_hub_height,
        wd=[270.0],
        ws=[10.0],
        yaw_angles=yaw_angles
    )

    frame = plot_yaw(negative, fi)
    frames.append(frame)

# Specify the duration between frames (milliseconds) and save to file:
gif.save(frames, 'example.gif', duration=100)

[hikilinkuja]

Hello @rafmudaf ! Tried running your gif but with no success. I'm kind of a noobie. Would you give me a hand?

Here's the error:

C:\Users\guilherme.faria\AppData\Local\anaconda3\envs\py310\python.exe "C:\Users\guilherme.faria\Desktop\docs\floris-3.4.1\examples\floris gif.py" 
Traceback (most recent call last):
  File "C:\Users\guilherme.faria\Desktop\docs\floris-3.4.1\examples\floris gif.py", line 60, in <module>
    frame = plot_yaw(negative, fi)
  File "C:\Users\guilherme.faria\AppData\Local\anaconda3\envs\py310\lib\site-packages\gif.py", line 55, in inner
    plot(*args, **kwargs)
  File "C:\Users\guilherme.faria\Desktop\docs\floris-3.4.1\examples\floris gif.py", line 37, in plot_yaw
    wakeviz.plot_turbines_with_fi(axes, fi)
  File "C:\Users\guilherme.faria\AppData\Local\anaconda3\envs\py310\lib\site-packages\floris\tools\visualization.py", line 97, in plot_turbines_with_fi
    yaw_angles = fi.floris.farm.yaw_angles
AttributeError: 'Axes' object has no attribute 'floris'

After installing gif on the correct environment all i did was change fi = FlorisInterface("gif.yaml") to fi = FlorisInterface("inputs/gch.yaml")

Thanks!

[rafmudaf]

That's my bad @hikilinkuja, I had an error in this line: wakeviz.plot_turbines_with_fi(fi, axes). Previously, the arguments were swapped. It's fixed above, so give it another try.

[donghyeoncheon]

@rafmudaf Hello.. I tried to follow your code.. do i need to change Floris Interface address? like ("Inputs/gch.yaml")..? gif.yaml.. how can i define that file..
Thank you

[rafmudaf]

@donghyeoncheon yes, you need to provide a path to your input file. If you're new to FLORIS, see Quick Start, Background and Concepts, and the input file reference.