MainDCScript_DELTA.m

% Script for Generating Diurnal Cycle Plots for JGR Manuscript
% AUTHOR: Jacob Wenegrat
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%
close all
clc

% Define the time range
trades = 1:2061;
calm = 2062:5373;
coldtongue = 5374:5757;
tiwwarm = 430:813; %Oct. 31 - Nov. 16th (16 days)
tiwcold = 46:429;
detail1 = tiwcold((13:(13+24+24*7))+(24*7));
detail2 = calm((13:(13+24+24*7))+(24*21));


noisebase = 0.1693;

% Define Variables
tempbins = -.5:.0125:.5; % Used for shear temperature anamoly colorscale
sst = dataset.density.fulldepth.temps(1,:)';

% Define Av variables
shu = (dataset.deepadcp.fullu(:,1) - dataset.deepadcp.fullu(:,6))./3.75;
shv = (dataset.deepadcp.fullv(:,1) - dataset.deepadcp.fullv(:,6))./3.75;
shmag = sqrt(shu.^2 + shv.^2);
avmag = dataset.measures.tauh./(1023*shmag);
% Calculate percentile
avcut = bootstrappercentile(avmag, .01, 1000);
av = avmag;
av(av>avcut) = NaN;

disp(['Trades dates: ', datestr(dataset.descriptors.datesen(trades(1))),'  -  ', datestr(dataset.descriptors.datesen(trades(end)))]);
disp(['Variable dates: ', datestr(dataset.descriptors.datesen(calm(1))),'  -  ', datestr(dataset.descriptors.datesen(calm(end)))]);
disp(['TIW Cold dates: ', datestr(dataset.descriptors.datesen(tiwcold(1))),'  -  ', datestr(dataset.descriptors.datesen(tiwcold(end)))]);
disp(['TIW Warm dates: ', datestr(dataset.descriptors.datesen(tiwwarm(1))),'  -  ', datestr(dataset.descriptors.datesen(tiwwarm(end)))]);

tri=triang(48-1)./24; %Following Kessler PDF

cmodDC = complex_demodulate_tri(sst, 1/24, tri);
% trilong = triang(48*6-1)./(24*6);
% cmodDC = complex_demodulate_tri(dataset.density.fulldepth.tempsraw(:,1), 1/(24*6), trilong);

cmodDC.semimaj = interp1(dataset.descriptors.datesend, deltaDaily(sst,dataset.descriptors.datesen)/2, dataset.descriptors.datesen);

disp(['Std Error on Monthly SSTA=0.28:   ', num2str(bootstrapErrors(noisebase, 1000, 0.28,30))]);

disp(['Average SSTA Trades: ', num2str(nanmean(cmodDC.semimaj(trades)))]);
disp(['Average SSTA Variable: ', num2str(nanmean(cmodDC.semimaj(calm)))]);

disp(['Average Wind Speed Trades: ', num2str(nanmean(dataset.measures.wspdh(trades)))]);
disp(['Average Wind Speed Variable: ', num2str(nanmean(dataset.measures.wspdh(calm)))]);

%Display instrumental errors
InstrumentalErrorCalcs;

%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Background Work
% Stuff that is not included in figures.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%==========================================================================
% Phase Info
% Look at the histograms/distributions of phase info for peak SST signal.
%==========================================================================
calculatePhaseInfo;
% shearratiocompare;


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Section: Introduction/Methods
% 
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%==========================================================================
% Figure 1 
% Map of PIRATA array and location of mooring
%==========================================================================
piratachart_wSST; %Still needs finalizing.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Section: Seasonal Cycle
%
% Long time series
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%==========================================================================
% Figure 2
% Long Time series of dSST and atmos variables
%==========================================================================
plot23Long(data23);


%==========================================================================
% Figure 3
% Climatology: Wind, SWR, dSST
%==========================================================================
plot23ClimSWR_wError(data23)


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Section: Enhanced Monitoring Period
%
% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
loadBuoyancyFlux;


%==========================================================================
% Figure 4
% Detailed overview panel.
%==========================================================================
Plot23wShortOver;

%==========================================================================
% Figure 5
% Detail plot of diurnal cycle
%==========================================================================
plotDetails(dataset, detail1);

%==========================================================================
% Figure 6
% Shear Plot: Trades
%==========================================================================
trange = trades;
makeShearPlot(dataset.variables.Emmpersec20mh, dataset.variables.Nmmpersec20mh, ...
    dataset.measures.uwndh, dataset.measures.vwndh, dataset.flux.netheat, ...
    dataset.density.fulldepth.fullmldh, av, dataset.deepadcp.riforwardreduce_sort, ...
    dataset, trange, tempbins);


%==========================================================================
% Figure 7
% Ri Plot: Trades
%==========================================================================
trange = trades;
limits = [0 .5];
figure
GradRI2_Stats;
set(gcf, 'Position', [0 0 656 352]);

%==========================================================================
% Figure 8
% Shear Plot: Variable
%==========================================================================

trange = calm;
makeShearPlot(dataset.variables.Emmpersec20mh, dataset.variables.Nmmpersec20mh, ...
    dataset.measures.uwndh, dataset.measures.vwndh, dataset.flux.netheat, ...
    dataset.density.fulldepth.fullmldh, av, dataset.deepadcp.riforwardreduce_sort, ...
    dataset, trange, tempbins);


%==========================================================================
% Figure 9
% Ri Plot: Variable
%==========================================================================
trange = calm;
limits = [0 .5];
figure
GradRI2_Stats;
set(gcf, 'Position', [0 0 656 352]);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Section: EMP - TIW
%
% Brief focus on TIW
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%==========================================================================
% Figure 10
% Overview of TIW: V, Sh^2, N^2, Rsh^2, with density contours
%==========================================================================
plotTIWOver;

%==========================================================================
% Figure 11
% Ri plots (2x): N-S phases.
%==========================================================================

limits = [0 .7];

figure
gap = [.05 .1]; margh= .1; margw = .2;
subtightplot(2,1,1, gap, margh, margw)
trange = tiwcold;
GradRI2_Stats;
set(gca, 'XTickLabel','');
set(get(gca, 'Xlabel'), 'String', '');
set(cb, 'YTick', [0:.1:limits(end)])
axes_label('a)', 10, 60);
subtightplot(2,1,2, gap, margh, margw)
trange = tiwwarm;
GradRI2_Stats;
set(get(gca, 'Title'),'String', '');
axes_label('b)', 10, 60);
set(gcf, 'Position', [0 0 665 740]);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Section: Marginal Instability
%
% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%==========================================================================
% Figure 12
% MI Plot
%==========================================================================
plotMIwithRI;


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Supplementary Information
%
% Error analysis
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Saving Figures
%
% Instructions for manual changes to figures prior to saving.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Fig 1
%%%%%%%%
% Move colorbars up and align to grid.

%Fig 2
%%%%%%%%
% Adjust labels (use align grid)

%Fig 4
%%%%%%%%
% Lots of adjusting, can try loading .fig file
% Resize to:
% set(gcf, 'Position', [-185 2000 1209 1028]);
% Adjust all labels throughout
%


%Fig 5
%%%%%%%%
% Move to big screen, adjust size
% set(gcf, 'Position', [-185 2000 1270 717]);
% Move labels

%Fig 6
%%%%%%%%
% Move to big screen, adjust size
% set(gcf, 'Position', [-185 2000 1174 1128]);
% Adjust vector legend text location

%Fig 8
%%%%%%%%
% Move to big screen, adjust size
% set(gcf, 'Position', [-185 2000 1174 1128]);
% Adjust vector legend text location

%Fig 10
%%%%%%%%
% Move to big screen, adjust size 
% set(gcf, 'Color', 'w', 'Position', [-185 2000 1340 1000]);
% Move labels
% Shift X label to left.

%Fig 11
%%%%%%%%
% Position (a) and (b) labels.

%Fig 12
%%%%%%%%
% set(gcf, 'Color', 'w', 'Position', [0 2000 1400 900])
% Move labels.
% Move 'Sh^2_red label up a bit.


% End Manuscript