The area of the ellipse is $\pi a b$, so given that you know the period you know $\frac {dA}{dt}= \frac {\pi a b}{\text{1 year}}$ Then $\frac {d\theta}{dt}=\frac 2{r^2}\frac {dA}{dt}=\frac {2 \pi a b}{r^2\text{1 year}}$ Now use the orbit equation to replace $r$ with a function of $\theta$ and you have an equation you can integrate, at least numerically.