I am interested in the interpolation method used in the meta.top function. In the Read et al., (2011) paper it states that the method is designed to; "interpolate between i and i+1 to yield the approximate depth of the base of the mixed layer". That is, the method interpolates between the depth where the threshold is first met and the consecutive depth. However, I find in the code the interpolation is actually between the depth before the threshold was met and the depth of the depth of the thermocline. Is this possibly a mistake?
Hi,
I am interested in the interpolation method used in the meta.top function. In the Read et al., (2011) paper it states that the method is designed to; "interpolate between i and i+1 to yield the approximate depth of the base of the mixed layer". That is, the method interpolates between the depth where the threshold is first met and the consecutive depth. However, I find in the code the interpolation is actually between the depth before the threshold was met and the depth of the depth of the thermocline. Is this possibly a mistake?
This is the specific line of code:
metaTop_depth = stats::approx(drho_dz[i:thermo_index], sortDepth[i:thermo_index], slope)
which could be changed to:
metaTop_depth = stats::approx(drho_dz[i:i+1], sortDepth[i:i+1], slope)
This is the original full function:
Thanks, Harriet