import xskillscore as xs
weights = np.cos(np.deg2rad(obs2.lat))
_, weights = xr.broadcast(obs2, weights)
weights = weights.isel(time=0)