test_steady_state_M98.test_steady_state_M98

test_steady_state_M98.test_steady_state_M98(x, y)

Performs the non-stationary ratio test described by Mahrt (1998) after Vitale et al. 2020, Biogeosciences, “Robust data cleaning procedure for eddy covariance flux measurements”” and associated inst_prob_test routine found on Gitlab

Parameters:

• x (raw high-frequency time series of vertical wind component (i.e. W))

• y (raw high-frequency time series of scalar atmospheric variable (e.g. CO2))

Returns:

• S_M98 (Returns the non-stationary ratio test statistic.)

• comments

• --------

• - computation of covariances might be more efficient using rolling and cov functions – on df (e.g. https://www.geeksforgeeks.org/how-to-calculate-rolling-correlation-in-python/), or the np.roll function

• Written by B. Heinesch, 29 October, 2022.

• University of Liege, Gembloux Agro-Bio Tech.