A computational method for model reduction in index-2 dynamical systems for Stokes equations. (1st October 2021)