Function es::fit_lewkowicz_speed(const Eigen::ArrayXd&, const Eigen::ArrayXd&, bool)

• Defined in File fit.h

Function Documentation

Array5d es::fit_lewkowicz_speed(const Eigen::ArrayXd &z, const Eigen::ArrayXd &u, bool print_report = true)

Determine best fit of Lewkowicz relation to input speed profile.

Robustly fits Lewkowicz speed relation \( u = f(z, \Pi, \kappa, U_{inf}, S, \delta_{c}) \) to input data \( u, z \), using ‘sensible’ initial guesses (for Atmospheric Boundary Layer in Northern European conditions) and fixed parameters to constrain the fitting process:

- von Karman constant \form#3 is fixed according to the value in definitions.h.
- Boundary layer thickness \form#24 is fixed to a default value of ``1000 m``.

Return

Array5d containing fitted values [pi_coles, kappa, u_inf, shear_ratio, delta_c]

Parameters

• z: Vertical locations in m (or normalised, if z_ref = 1.0)
• u: Observed speeds at corresponding vertical locations, in m/s
• print_report: bool, default true. If true, the solver prints iterations summary and full final convergence and solution report (for debugging and validation purposes).