MAST-467/667: Introduction to Arctic Oceanography (2021)

HW 10 - Tmax below Mixed Layer

% The next step will be to use our bin average profile data of potential temperature, absolute salinity,

% and potential density to estimate the potential temperature maximum below the surface mixed layer (SML)

% as defined by the bin above the first stability maximum N2max

load /Users/fbenz/Desktop/PhD/courses/1sem/arctic_seminar/9HW/bin_av_data_profile31_binsize1_zSML.mat

binDepth = binDepth';

clf,

plot(binSA,binDepth,'linewidth',2)

ylabel('Depth [m]')

xlabel('Absolute Salinity [ppt]')

set(gca,'ydir','reverse','linewidth',2,'fontsize',15)

clf,

plot(binCT,binDepth,'linewidth',2)

ylabel('Depth [m]')

xlabel('Conservative Temperature [deg. C]')

set(gca,'ydir','reverse','linewidth',2,'fontsize',15)

clf,

plot(binP,binDepth,'linewidth',2)

ylabel('Depth [m]')

xlabel('Pressure [dbar]')

set(gca,'ydir','reverse','linewidth',2,'fontsize',15)

% Buoyancy frequency

[N2, p_mid] = gsw_Nsquared(binSA,binCT,binP);

clf,

plot(N2(1,2:end),p_mid(1,2:end),'linewidth',2)

b = hline(zSML,'linewidth',2);

legend([b],'SML')

ylabel('Depth [m]')

xlabel('N^2 [1/s^2]')

set(gca,'ydir','reverse','linewidth',2,'fontsize',15)

% What is the depth and salinity of the potential temperature maximum Tmax, below the depth of the first

% stability maximum?

% fist stability maximum

maxN2 = max(N2);

ind_maxN2 = find(N2 == maxN2); % index of the first stability maximum

% Pot. temperature maximum (Tmax), below the first stability maximum

Tmax = max(binPT(ind_maxN2+1:end)) % +1 bin

Tmax = 4.7180

ind_Tmax = find(binPT == Tmax);

% Depth and salinity at Tmax

depth_Tmax = binDepth(ind_Tmax)

depth_Tmax = 26.0453

SA_Tmax = binSA(ind_Tmax)

SA_Tmax = 33.7685