pydune.physics.dune.courrechdupont2014.flux_at_crest#

flux_at_crest(alpha, theta, Q0, gamma=1.6)[source]#

Compute the sand flux at the dune crest:

\[Q_{\rm crest} = Q_{0}\left[1 + \gamma\sin\vert\theta-\alpha\vert\right]\]
Parameters:

  • alpha (scalar, numpy array) – dune orientation \(\alpha\).

  • theta (scalar, numpy array) – flux orientation \(\theta\) in degrees.

  • Q0 (scalar, numpy array) – flux at the bottom of the dune \(Q_{0}\).

  • gamma (scalar, numpy array) – flux-up ratio \(\gamma\) (the default is 1.6).

Returns:

flux at the dune crest \(Q_{\rm crest}\)

Return type:

scalar, numpy array

Examples

>>> import numpy as np
>>> alpha = 10
>>> theta = np.random.random((1000,))*360
>>> Q0 = np.random.random((1000,))*50
>>> Qcrest = flux_at_crest(alpha, theta, Q0)