PLOT_BINS

Torsten Linders, torsten@oceansurface.se, January 4, 2012

  1. Plots dissipation from matices representing horizontal bins during one day.
  2. Plots dissipation from matices of S2 and N2 representing all observations during one leg, day or all three days.

Contents

Cleaning up

clear
close all

Settings

par='eNSav';                % eps, eNSav
par_vector=0;               % 0

Loading prepared data

load mat/timeposlegs
load mat/fjordsection
load mat/mycolormap

if strcmp(par,'eps')
    load mat/sigma sigma_v
    load mat/epsab
elseif strcmp(par,'eNSav')
    load mat/eNSav
end

Plotting

for avday=1:4
    fig=figure(100+avday);
    set(fig,'Position',[1 1 800 900]);
    set(fig,'Colormap',mycolormap);

    if avday<4
        day_title=['Day ' num2str(avday)];
    elseif avday==4
        day_title='All legs';
    end

    if strcmp(par,'eps')
        imagesc(epsab(:,:,avday))
        title([day_title ', average dissipation (log_{10}[W/kg])'])
    elseif strcmp(par,'eNSav')
        subplot(2,1,1)
        imagesc(eNSav(:,:,avday))
        title([day_title ', average dissipation (log_{10}(W/kg)) at ' ...
            num2str(dmin) ' - ' num2str(dmax) ...
            ' m, averaging of un-logarithmized values'])
        caxis([-9 -5])
        plot_extras
        subplot(2,1,2)
        imagesc(log10(eNSavi(:,:,avday)))
        title([day_title ', number of observations (log_{10})'])
        caxis([0 5])
    end

Adjustments and saving

    if strcmp(par,'eps')
        plot_extras_bins
    elseif strcmp(par,'eNSav')
        plot_extras
    end
    printplot_av
    close all
end