Variable eddy viscosity Ekman layer in the ABL (1D)

Ekman's equations in the atmospheric boundary layer are solved for a horizontally homogeneous flow and a height-dependant eddy viscosity

Ekman1D

Variable eddy viscosity Ekman layer in the ABL (1D)

Summary

Matlab implementation of the solution to the Ekman equations in the atmospheric boundary layer. The flow is assumed horizontal and homogeneous. however, a height-dependant eddy viscosity can be modelled. The solutions are provided in one-dimension.

Content

The submission includes

• The function EkmanAnalytic that provides analytics solution of Ekman's equations for a constant eddy viscosity in the atmospheric boundary layer.
• The function solveEkman that numerically solves Ekman's equations with an explicit finite difference scheme and allows the use of height-dependant eddy viscosity. The numerical implementation is partly inspired by [1].
• An example file Example0.mlx and reproduces some of the figures displayed in ref [2]

Any question, suggestion or comment is welcomed.

References:

[1] Berger, B. W., & Grisogono, B. (1998). The baroclinic, variable eddy viscosity Ekman layer. Boundary-layer meteorology, 87(3), 363-380.

[2] https://www.whoi.edu/science/PO/people/jprice/website/education_scripts.html