DBYr = np.array((y + w * np.sin(d + math.pi)))