Xskillscore

import xskillscore as xs

weights = np.cos(np.deg2rad(obs2.lat))

_, weights = xr.broadcast(obs2, weights)

weights = weights.isel(time=0)