Add planetary evolution functionality: past HZ occupancy and UV flux

[matiscke]

This adds a method evolve() to the Table class, which adds a new attribute evolution to the Table. evolution includes for each planet the time evolution of the habitable zone and its received NUV flux.

The evolution attribute contains for each planet:

• 'time' : time grid in Gyr
• 'lum' : luminosity evolution in L_sun
• 'nuv' : NUV flux evolution in erg/s/cm^2
• 'in_hz' : boolean array indicating whether the planet is in the HZ at each time step

There is further a new function past_hz_uv that determines which planets were both in the HZ and had NUV fluxes above a threshold value for more than a specified minimum time range.

[matiscke]

The current way to simulate measurement errors, drawing from a log-normal distribution (to reflect the scatter of the NUV data in Richey-Yowell et al. (2023) at each point in time leads to extreme outliers.

This is problematic, as we may neet to estimate the maximum past NUV flux, which can then take on extremely high values (especially for large sigma_nuv_dex, see right panel (log scale!)).

I think it would be better to draw a bias for each planet that is applied evenly to the whole time evolution.

[matiscke]

better.