y = np.linspace(-csz/2., (nrows-1)*csz - csz/2., nrows)