function master_mblt(theYears)
%
%   master_mblt(theYears)
%
%   Main control function for generating the MBLT report.  Only creates new
%   plots for theYears.

%Soupy Alexander, 10/27/03

dataDirectory = 'C:\SCHTUFF\MASS_BAY\MBLT_REPORT\DATAFILES';
plotDirectory = 'C:\SCHTUFF\MASS_BAY\MBLT_REPORT\PLOTS';
outDirectory = 'C:\SCHTUFF\MASS_BAY\MBLT_REPORT\OTHER';

%Create an archive of the data files (or load it in)
%archive_mblt(dataDirectory)
load(fullfile(dataDirectory, 'mbltArchive.mat'))

%%DEEPEST ADCP BIN%%
adcp_deep(dataDirectory, outDirectory, siteA, siteB)

%%DATA AVAILABLE%%
data_available(theYears, dataDirectory, plotDirectory, siteA, siteB)    

%%MONTHLY STATS%%
monthly_stats(dataDirectory, plotDirectory)
monthly_stats2(dataDirectory, plotDirectory, siteA, siteB)

%%CURRENTS BY SPEED AND DIRECTION%%
sort_currents(theYears, dataDirectory, plotDirectory, siteA, siteB)

%%TIDAL ANALYSIS%%
tides(theYears, dataDirectory, plotDirectory, SiteA, siteB)

%%TIME SERIES PLOTS%%
%Plots of meteorological observations
weath_plot(theYears, dataDirectory, plotDirectory)

%%Mixed Plot
mixed_plot(theYears, dataDirectory, plotDirectory)
mixed_plot_b(theYears, dataDirectory, plotDirectory, siteB)

%Plots of streamflow, salinity, and temperature (Sites A and B)
river_ts(theYears, dataDirectory, plotDirectory)

%Plot of wind stress and point current meter whisker plots from Site A
pcm_whisk(theYears, dataDirectory, plotDirectory)

%Plot of wind stress and profiling current meter whisker plots
adcp_whisk(theYears, dataDirectory, plotDirectory, siteA, 'a')
adcp_whisk(theYears, dataDirectory, plotDirectory, siteB, 'b')

%Plot of the pressure
pres_plot(theYears, dataDirectory, plotDirectory, siteA, siteB)

%Beam attenuation
attn(theYears, dataDirectory, plotDirectory)

%ADCP contour plots
adcp_contour(theYears, dataDirectory, plotDirectory, siteA, 'A')
adcp_contour(theYears, dataDirectory, plotDirectory, siteB, 'B')

%Velocity profile viewer
adcp_vpv(theYears, dataDirectory, plotDirectory, siteA, 'A')
adcp_vpv(theYears, dataDirectory, plotDirectory, siteB, 'B')

%%PRINCIPAL AXES%%

%Point current meter scatter plots
pcm_scat_plots(theYears, dataDirectory, plotDirectory)
%ADCP scatter plots
adcp_scat_plots(theYears, dataDirectory, plotDirectory, siteA, 'A')
adcp_scat_plots(theYears, dataDirectory, plotDirectory, siteB, 'B')

%Ellipse plots (PCM)
pcm_ell_plots(theYears, dataDirectory, plotDirectory)
adcp_ell_plots(theYears, dataDirectory, plotDirectory, siteA, 'A')
adcp_ell_plots(theYears, dataDirectory, plotDirectory, siteB, 'B')

%Principal axes
pcm_ell_tables(theYears, dataDirectory, plotDirectory)
adcp_ell_tables(theYears, dataDirectory, plotDirectory, siteA, 'A')
adcp_ell_tables(theYears, dataDirectory, plotDirectory, siteB, 'B')

%Site Comparison
compare_sites(theYears, dataDirectory, plotDirectory, siteA, siteB)

%Sediment trap plots
sed_trap(theYears, dataDirectory, plotDirectory)
sed_trap_2(max(theYears), dataDirectory, plotDirectory)
sed_trap_month(max(theYears), dataDirectory, plotDirectory)