View source for Macro Deodatis 1996 Wind Simulation Points

{{Macro
|Name=Macro Deodatis 1996 Wind Simulation Points
|Icon=Macro_Deodatis_1996_Wind_Simulation_Points.png
|Description=This macro create simulation point according to G. Deodatis 1996.
|Author=Koffi Daniel
|Version=1.0
|Date=2024-04-15
|RPSVersion=All
|Download=[https://wiki.freecad.org/images/f/f5/Macro_3D_Parametric_Curve.png ToolBar Icon]
|SeeAlso=[[Macro_Deodatis_1996_Mean_Wind_Profile|Macro Deodatis 1996 Mean Wind Profile]]
}}
This macro create simulation point according to G. Deodatis 1996. There are three simulation points created with WindLabFeatureLocationDistribution feature.
[[File:Deodatis1996WindSimulationPoints.png|600px]]
{{Caption|Simulation Point Distribution}}

[[File:Macro_Deodatis_1996_Wind_Simulation_Points_Heighs.png|600px|Points Heights‎|center]]

==Script== 

You can find the source code on the following Github repository:
[https://github.com/LabRPS/LabRPS-macros/blob/master/WindLab/Deodatis1996/Deodatis1996_LocationDistribution.RPSMacro Get the code here!]

'''Macro Deodatis 1996 Wind Simulation Points.FCMacro'''

{{MacroCode|code=
 
# -*- coding: utf-8 -*-
# (c) Koffi Daniel 2024

import LabRPS
from LabRPS import Vector as vec

import WindLab
import WindLabObjects

#install the WindLab plugin with all its plugged features
installResuslt = WindLab.installPlugin("WindLabPlugin")

doc = LabRPS.newDocument()

# create new simulation with default name "Simulation"
sim = WindLabObjects.makeSimulation(doc)
    
# check if the simulation is successfully created
if not sim:
    LabRPS.Console.PrintError("The simulation does not exist.\n")

# create new location distribution in the simulation called "Simulation"
loc = WindLabObjects.makeFeature("SimulationPoints", "Simulation", "Location Distribution",  "General Distribution")
    
# check if the created location distribution feature is null
if not loc:
    LabRPS.Console.PrintError("Error on creating the location distribution.\n")

v1 = vec(0, 0, 35)
v2 = vec(0, 0, 40)
v3 = vec(0, 0, 140)
loc.Locations = [v1, v2, v3]

# run the method "computeLocationCoordinateMatrixP3" of the feature and store the result.
# Note that here, this method will be called for the active location distribution feature in case they are more than one
locRes = sim.computeLocationCoordinateMatrixP3()

if LabRPS.GuiUp:
    import WindLabGui
    #show the location coordinate in Alphaplot
    WindLabGui.WindLabPythonPluginUtilities.showArray(sim.getSimulationData().numberOfSpatialPosition, 4, locRes)
    loc.ViewObject.Document.activeView().viewAxonometric()
    loc.ViewObject.Document.activeView().fitAll()

    import numpy
    import matplotlib.pyplot as plt
    arr = numpy.asarray(locRes)
    fig, ax = plt.subplots()
    fig.suptitle('Location Distribution')
    ax.plot(arr[:,0], arr[:,3], 'o-')
    ax.set_xlabel('Simulation Point Number')
    ax.set_ylabel('Height (m)')
    ax.grid(True)
    plt.show()
  
doc.recompute()

}}

{{clear}}