gael/profilesMatlabProcessing/find_depth.m

function DEPTH=depth(P,LAT);
% DEPTH   Computes depth given the pressure at some latitude
%         D=DEPTH(P,LAT) gives the depth D (m) for a pressure P (dbars)
%         at some latitude LAT (degrees).
%
%         This probably works best in mid-latiude oceans, if anywhere!
%
%         Ref: Saunders, Fofonoff, Deep Sea Res., 23 (1976), 109-111
% 

%Notes: RP (WHOI) 2/Dec/91
%         I copied this directly from the UNESCO algorithms

% CHECKVALUE: DEPTH = 9712.653 M FOR P=10000 DECIBARS, LATITUDE=30 DEG
%     ABOVE FOR STANDARD OCEAN: T=0 DEG. CELSUIS ; S=35 (IPSS-78)
      X = sin(LAT/57.29578);
%**************************
      X = X.*X;
% GR= GRAVITY VARIATION WITH LATITUDE: ANON (1970) BULLETIN GEODESIQUE
      GR = 9.780318*(1.0+(5.2788E-3+2.36E-5*X).*X) + 1.092E-6.*P;
      DEPTH = (((-1.82E-15*P+2.279E-10).*P-2.2512E-5).*P+9.72659).*P;
      DEPTH=DEPTH./GR;


1	function DEPTH=depth(P,LAT);
2	% DEPTH Computes depth given the pressure at some latitude
3	% D=DEPTH(P,LAT) gives the depth D (m) for a pressure P (dbars)
4	% at some latitude LAT (degrees).
5	%
6	% This probably works best in mid-latiude oceans, if anywhere!
7	%
8	% Ref: Saunders, Fofonoff, Deep Sea Res., 23 (1976), 109-111
9	%
10
11	%Notes: RP (WHOI) 2/Dec/91
12	% I copied this directly from the UNESCO algorithms
13
14	% CHECKVALUE: DEPTH = 9712.653 M FOR P=10000 DECIBARS, LATITUDE=30 DEG
15	% ABOVE FOR STANDARD OCEAN: T=0 DEG. CELSUIS ; S=35 (IPSS-78)
16	X = sin(LAT/57.29578);
17	%**************************
18	X = X.*X;
19	% GR= GRAVITY VARIATION WITH LATITUDE: ANON (1970) BULLETIN GEODESIQUE
20	GR = 9.780318(1.0+(5.2788E-3+2.36E-5X).X) + 1.092E-6.P;
21	DEPTH = (((-1.82E-15P+2.279E-10).P-2.2512E-5).P+9.72659).P;
22	DEPTH=DEPTH./GR;
23
24