A parameterization with wave saturation adjustment of subgrid-scale average wave stress over three-dimensional topography