Parameterization of subgrid stirring in eddy resolving ocean models. Part 1: Theory and diagnostics