Sharp H1-norm error estimate of a cosine pseudo-spectral scheme for 2D reaction-subdiffusion equations