How to specify a state-dependent gustiness parameterization?

[glwagner]

The COARE 3.0 algorithm for air-sea fluxes uses a state-dependent parameterization for gustiness (Fairall et al 2003, equation 8):

Here the magnitude of gusts are estimated using the boundary layer depth z_i, combined with the virtual temperature flux at the surface (which is part of the iterative solution to the bulk system of equations). There is a free parameter beta.

In some objects, the gustiness can be given as an input, eg

https://github.com/CliMA/SurfaceFluxes.jl/blob/b60ffa4711a0e4b5b90ac08800a5a305c2e62981/src/SurfaceFluxes.jl#L211

But this does not allow a state-dependent gusintess parameterization.