dys = ycoord[1:n, 0:m-1]-ycoord[0:n-1, 0:m-1]