python_codes.meteo_analysis.quadratic_transport_law#

quadratic_transport_law(theta, theta_d, omega)[source]#

Quadratic transport law \(q_{\rm sat}/Q = \Omega \sqrt{\theta_{\rm th}}(\theta - \theta_{\rm th})\), from Duràn et al. 2011.

Parameters
  • theta (scalar, numpy array) – Shield number.

  • theta_d (scalar, numpy array) – Threshold Shield number.

  • omega (scalar, numpy array) – Prefactor of the transport law.

Returns

Sediment fluxes calculated elementwise using the quadratic transport law.

Return type

scalar, numpy array

Examples

>>> import numpy as np
>>> theta = np.random.random((2000, ))
>>> theta_d, omega = 0.0053, 7.8
>>> qsat = quadratic_transport_law(theta, theta_d, omega)

References

[1] Durán, O., Claudin, P., & Andreotti, B. (2011). On aeolian transport: Grain-scale interactions, dynamical mechanisms and scaling laws. Aeolian Research, 3(3), 243-270.

Examples using python_codes.meteo_analysis.quadratic_transport_law#