We need to do a new huddle test of the Guralps for the Wiener filtering paper. The last test had miserable results.
I tried to use recent data to do this, but it looks like we forgot to turn the Guralp box back on after the power outage or that they're far off center.
So instead I got data from after the previous power outage recovery.
I tried to use our usual Wiener filter method to subtract Guralp1-Z from Guralp2-Z, but that didn't work so well. It was very sensitive to the pre-weighting.
Instead I used the new .m file that Dmass wrote for subtracting the phase noise from his doubling noise MZ. That worked very well. It does all of the subtraction in the frequency domain and so doesn't have to worry about making a stable or causal filter. As you can see, it beats our weighted Wiener filter at all frequencies.
The attached plot shows the Guralp spectra (red & green), the residual using time-domain Wiener filtering (black) and the Dmass f-domain code (yellow).
As soon as Jenne brings in her beer cooler, we're ready to redo the Huddle Test.