pressureRange: [0, max(Pressure)] (default)
profile: [ ] (all profiles default).
direction: up, down or both.
windowLength: 21 (default)
This function estimates the optimal conductivity time shift relative to temperature in order to minimise salinity spiking. The algorithm works by computing the salinity for conductivity lags of -20 to 20 samples at 1 sample increments. The optimal lag is determined by constructing a high-pass filtered version of the salinity time series for every lag, and then computing the standard deviation of each. The optimal lag is the one that yields the smallest standard deviation.
pressureRange argument allows for the possibility of estimating the lag on specific sections of the profile. This can be useful when unreliable measurements near the surface are found to impact the optimal lag, or when the profiling speed is highly variable.
The lag output by
RSKcalculateCTlag is compatible with the
lag input argument for
lagunits is specified as
rsk = RSKopen('file.rsk'); rsk = RSKreadprofiles(rsk, 'profile', 1:10, 'direction', 'down'); % 1. All downcast profiles with default smoothing lag = RSKcalculateCTlag(rsk); % 2. Specified profiles (first 4), reference salinity found with 13 pt boxcar. lag = RSKcalculateCTlag(rsk, 'profile',1:4, 'windowLength',13);