pressure_to_depth_sea_bisset.m

function D = pressure_to_depth_sea_bisset(P,L)
% Inputs
%   P: pressure \ kPa
%   L: latitude \ degree \ -90 < L < 90
% Outputs
%   D: depth \ m

    P = P.*0.0102;
    D = 9.7512*P./(1+(5.3e-3)*sind(L).^2) - (2.07e-4)*P.^2;
end