pydune.physics.turbulent_flow.andreotti2011_unbounded.calculate_solution

calculate_solution(eta_0, eta_H, alpha, max_z=None, Kappa=0.4, atol=1e-10, rtol=1e-10, method='DOP853', **kwargs)