Output: 4D matrix u = u[3,Nx,Ny,Nz] with
* size(u,1) = turbulent components u, v and w;
* size(u,2) = number of points Nx in the longitudinal direction, x ;
* size(u,3) = number of points Ny in the transversal direction, y ;
* size(u,4) = number of points Nz in the vertical direction, z ;
When defining the grid points Nx,Ny and Nz, the user has to input noteven Ny and Nz in order to count already the hub point.
INPUT
%Box size
Nx= Number of grid points on xaxis;
Ny= Number of grid points on yaxis;
Nz= Number of grid points on zaxis;
%Spatial grid resolution
delta;
%WT / turbulence condition
UHub= Mean Wind Speed at the Hub Point;
HubHt= Hub Height;
z0= Terrain Roughness;
IECturbC= Wind Turbine Turbulence Class;
