On formulating two‐dimensional vertically integrated hydrodynamic numerical models with an enhanced representation of bed stress