dxs = xcoord[1:n, 0:m-1]-xcoord[0:n-1, 0:m-1]