I have started a python implementation of the AlGaAs noise budget. All parameters, functions, etc. are defined in a single notebook, and this same notebook generates the plot. The python uncertainties package facilitates estimation of uncertainties in material parameters, optical parameters, etc.
Currently, the coating thermo-optic trace is not an actual calculation; it is just a flat line culled from figure 5.9 of Tara's thesis.
The PDH shot noise trace is shown assuming an incident power of 1 mW on each cavity, a PDH modulation index of 0.2 rad, and a cavity visibility of 0.92.