helenjanisz / atacr Goto Github PK

When I use Bell et al. (2015) methods to calculate the tilt angle, I calculate the coherence between H component and Z component.
i.e., abs(specprop.rotation.chz_stack).^2./(specprop.rotation.chh_stack.*specprop.power.czz_stack) in ATaCR
Some days the coherence between H and Z will be greater than 1 (Figure 1), which does not conform to the actual situation.

Figure 1. Tilt angle (bottom) and corresponding coherence (top) during ocean-bottom station deployment. The red line in the top subfigure indicates coherence = 1.

Thank you for the code provided. When I calculated the transfer function, I found that there is an error in function noisecal_dailystaspectra.m, which will cause all transfer functions related to the H component to be incorrect. This problem is caused by the wrong number of windows used for averaging.
I solved this problem as follows:
The code in line 140

I found an error in function findTFs_correctevent that caused the previous day's transfer function file to be used when the current day's transfer function file exists.
I solved this problem as follows:

change to iday = find(dday >= 0);
I hope this is useful for the ATaCR users.
Thanks.