A toolkit for processing Seaglider base station NetCDF files: despiking, smoothing, outlier detection, backscatter, fluorescence quenching, calibration, gridding, interpolation.
Describe the bug
Mixed layer depth computation currently only works if the stratification is stable and the variable used for computation is monotonically increasing with depth.
Pass temperature as input variable, use thresh=0.3, ref_depth=1 (problem is present for any parameters though...)
See error, the mixed layer depth is very noise and not at the same depth as the thermocline for a temperature threshold.
Expected behavior
Mixed layer depth computation that takes the shallowest observation below the reference depth that exceeds the given threshold. Currently it takes the value with is the closes to (variable at reference depth - threshold variable), which can be at any depth.
Screenshots
Additional context
Please note that I already implemented a corrected algorithm. I create this Bug report for documentation and motivation for my pull request to come soon. :)
Describe the bug Mixed layer depth computation currently only works if the stratification is stable and the variable used for computation is monotonically increasing with depth.
To Reproduce Steps to reproduce the behavior:
Expected behavior Mixed layer depth computation that takes the shallowest observation below the reference depth that exceeds the given threshold. Currently it takes the value with is the closes to (variable at reference depth - threshold variable), which can be at any depth.
Screenshots
Additional context Please note that I already implemented a corrected algorithm. I create this Bug report for documentation and motivation for my pull request to come soon. :)