I ran the code for aLIGO coating structure and tried to reproduce fig5 and fig7 of Hong et al. paper. It turned out that the derivatives of the complex reflectivity were not matching with the paper. I rewrote the code, with a fresh mind without looking at previous code and voila, after increasing computation time slightly due to more brute force calculations, I was able to reproduce the figures. These figures are attached. Since I do not have access to the data of the figures in the paper (I tried to email the authors but got no replies), I could only try to plot it on the same scale and limits and check the values by eyes. The values seem to match. So I am more confident now to declare that this code completely follows the paper's calculations.

However, this does not change the coating Brownian noise. I have updated the noise budgets at the Daily and Latest plots.