Figure reescaling
Hello everyone, I wrote this script to simulate Blockage and Wake coupling:
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import time
from py_wake.examples.data.hornsrev1 import Hornsrev1Site, HornsrevV80
from py_wake.deficit_models import *
from py_wake.wind_farm_models import *
from py_wake.rotor_avg_models import *
from py_wake.superposition_models import *
from py_wake.deflection_models import *
from py_wake.turbulence_models import *
from py_wake.ground_models import *
from py_wake.deficit_models.utils import *
from py_wake.examples.data.hornsrev1 import Hornsrev1Site, V80, wt_x, wt_y, wt16_x, wt16_y
from py_wake.examples.data.iea37 import IEA37Site, IEA37_WindTurbines
from py_wake.wind_turbines import WindTurbine, WindTurbines
from py_wake.examples.data.dtu10mw import DTU10MW
dtu10mw ='Mywindturbine'
site=IEA37Site()
windTurbines = DTU10MW()
x,y = ([565761.6,
566664.0,568428.4,570192.7,571095.2,571957.0,572859.5,573721.4,574623.8,575526.3,577290.6,565517.9,566420.4,578111.9,580819.3,565477.3,567282.3,571754.0,573625.3,575363.9,582583.6,583486.1,568169.1,569046.6,569949.1,574420.8,579835.6,581640.5,583445.5,584347.9,585250.4,586152.9,566298.6,569908.5,572772.3,575434.2,576185.1,579795.0,580697.5,565355.5,582502.4,583404.9,584307.3,585209.8,586112.3,588819.7,568965.4,565314.9,573437.1,577047.0,588779.1,589681.5,590584.0,567119.9,572534.7,575201.5,576103.9,577908.9,580616.3,582421.2,583323.7,585128.6,586056.1,587836.0,588785.9,589640.9,591445.9,567981.7,592348.3,570689.1,577868.3,578770.7,585990.5,586892.9,589600.3,590502.8,591405.3,577827.7,578730.1,579632.6,581437.5,582340.0,583242.5,584144.9,585949.9,586852.3,587754.8,588657.3,589559.7,590462.2,577787.1,578689.5,579592.0,585909.3,577746.4,578648.9,579551.4,580453.8,581356.3,582258.8,583161.3,584063.7,584966.2,585868.7,586771.1,588623.6,589478.5,586730.5,588535.5,582177.6,583080.0,583982.5,584885.0,587592.4],[6264687.0,6264165.9,6264616.7,6265067.4,6264546.3,6265518.1,6264997.1,6265968.8,6265447.8,6264926.8,6265377.5,6273643.8,6273122.8,6267842.1,6266278.9,6275136.6,6274094.5,6272982.2,6272014.8,6270898.0,6266729.7,6266208.6,6275109.6,6274545.3,6274024.2,6272911.8,6269785.6,6268743.5,6267701.4,6267180.4,
6266659.4,6266138.3,6277601.2,6275517.0,6275668.9,6274054.4,6273362.6,6271278.4,6270757.4,6279615.0,6269715.3,6269194.2,6268673.2,6268152.2,6267631.1,6266068.0,6277530.9,6281107.9,6276418.5,6274334.3,6267560.8,6267039.8,6266518.7,6280065.8,6276939.5,6276869.2,6276348.2,6275306.1,6273743.0,6272700.9,
6272179.9,6271137.8,6270573.4,6269574.7,6269055.0,6268532.6,6267490.5,6281037.5,6266969.5,6279474.4,6276798.9,6276277.9,6272109.5,6271588.5,6270025.4,
6269504.3,6268983.3,6278291.7,6277770.7,6277249.6,6276207.6,6275686.5,6275165.5,6274644.4,6273602.4,6273081.3,6272560.3,6272039.2,6271518.2,6270997.2,
6279784.5,6279263.5,6278742.4,6275095.2,6281277.3,6280756.3,6280235.2,6279714.2,6279193.2,6278672.1,6278151.1,6277630.0,6277109.0,6276588.0,6276066.9,
6275026.2,6274503.8,6277559.7,6276517.7,6281657.7,6281136.7,6280615.7,6280094.6,6278531.5])
t = time.time()
wfm = All2AllIterative(
site,
windTurbines,
wake_deficitModel=NOJDeficit(
k=0.1,
rotorAvgModel=None,
groundModel=None),
superpositionModel=LinearSum(),
blockage_deficitModel=VortexDipole(
limiter=1e-10,
sct = 1,
exclude_wake=True,
superpositionModel=None,
rotorAvgModel=None,
groundModel=None,
upstream_only=False),
deflectionModel=None,
turbulenceModel=STF2017TurbulenceModel(),
convergence_tolerance=1e-06,
rotorAvgModel=None)
sim_res = wfm(x, y)
fig, ax = plt.subplots(figsize=(12, 8))
wake_map = sim_res.flow_map(wd=270).plot_wake_map(cmap='viridis')
plt.title('AEP: %.2fGWh'%(sim_res.aep().sum()))
plt.show()
The issue here is the velocity scale is from 1m/S ton 10m/S so I can't visualize the Blockage effects well, can you please help me to do the reescaling from 0.5m/s to 1m/s !
Thank you
Best regards
Edited by Mads M. Pedersen