A two-equation turbulence model for canopy flows