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

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.

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