points[:, 1] = particles['y'] - 0*dem['originalHeader']['yllcenter']