File Exchange

image thumbnail

Analysis of turbulent wind velocities in the neutral and unstable surface-layer

version 1.1 (4.7 MB) by

A study of surface layer fluxes and wind velocity spectra (simulated data)



View License

A meteorological tower with 3 sonic anemometers at 3 different heights is "simulated" and 4 samples corresponding to both neutral and unstable atmospheric stratifications are provided for each height.
The goal is to use the eddy-covariance method to retrieve the different turbulence statistics and the associated flux of momentum and heat in the surface layer.
The submission contains

• 4 .mat files with the different samples (simulated time series)

• 2 examples files: One to illustrate the spectral model from Højstrup [1], the other one to study the momentum and heat flux and the wind velocity spectra.

• The functions “verticalSpectrum_Hojstrup” and “longitudinalSpectrum_Hojstrup” to compute the spectral model from Højstrup [1].

• The function “getStability” to estimate the Obukhov length, the heat flux, the scaling temperature and friction velocity.

• The function “obukhovLength” to compute the Obukhov length.

• The function “similarityFun” to get the non-dimensional vertical profile of heat and momentum.

• The functions “plotSimilarityFun” and “plotVelocitySpectra” to plot more clearly the velocity spectra and the non-dimensional profiles

There are many typographic mistakes in the present version of the submission, which will be progressively cleared during the next few weeks. This submission was done during my free time as an exercise to better learn and understand the Monin-Obukhov similarity theory (MOST). Therefore, any comments or suggestion to improve this submission is warmly welcomed.

[1] Højstrup, J. (1981). A simple model for the adjustment of velocity spectra in unstable conditions downstream of an abrupt change in roughness and heat flux. Boundary-Layer Meteorology, 21(3), 341-356.

Comments and Ratings (0)



Correct version loaded


The possibility to include the across-wind component "v" is added
Some typo are corrected + the description is improved




Description added


Title modified

MATLAB Release
MATLAB 8.4 (R2014b)

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video