Construction of recurrent bivariate fractal interpolation surfaces and computation of their box-counting dimension