  In this paper, we present a general framework to construct recurrent fractal
interpolation surfaces (RFISs) on rectangular grids. Then we introduce bilinear
RFISs, which are easy to be generated while there are no restrictions on
interpolation points and vertical scaling factors. We also obtain the box
dimension of bilinear RFISs under certain constraints, where the main
assumption is that vertical scaling factors have uniform sums under a
compatible partition.
