“…The set of the couples (U(s,h), V(s,h)) will be denoted by H*; (2) A 9 × 9 calculation grid G was defined over the spatial domain, and a Kriging technique of vector interpolation (boundary-constrained) was applied to H* to estimate the wind velocity components u(i,j,h) and v(i,j,h) at the nodes (i,j) of G for each h; (3) These estimations were used to calculate the wind velocity components (u, v), the wind divergence γ, and the wind vorticity ω at the cells of the 8 × 8 lattice, L, associated to the calculation grid G. At each cell (p,q) of L, the estimation of the parameters u(p,q), v(p,q), γ(p,q) and ω(p,q) was carried out using the values of (u,v) at the four nodes located at the cell vertexes. Here, the 2D numerical definitions [25,28] of divergence and vorticity were used. The set {(u(p,q,h), v(p,q,h), γ(p,q,h), ω(p,q,h)) | (p,q) ∈ L } will be denoted by W(h).…”