clear all;close all; clc;
% jfm 2018 
% fig. 3
% grid size 1024 x 1024 x 500 
% note : every 2nd point was written
%        out to data_fig3_flow-statistics.mat

load data_fig3_flow-statistics.mat

% data = 
%   struct with fields:
%
%       urms: [250×1 double]
%       wrms: [250×1 double]
%          z: [250×1 double]
%        L11: [250×1 double]
%     stdL11: [250×1 double]
%         Rt: [250×1 double]
%      stdRt: [250×1 double]
     
figure;
plot(data.urms, data.z,'k',...
    data.wrms, data.z,'--k');
xlim([0 0.4])
ylim([0 5])
xlabel('u_{rms}/U, w_{rms}/U')
legend('u_{rms}/U','w_{rms}/U')
set(legend,'position',[0.63 0.7 0.001 0.01],'box','off')
ylabel('z/L');

figure;
plot(data.L11, data.z,'k') ;
hold on;
plot(data.L11+data.stdL11, data.z,':k') ;
plot(data.L11-data.stdL11, data.z,':k') ;
xlabel('$$\mathrm{\overline{L_{11}}/L}$$','interpreter','latex')
ylabel('z/L')
xlim([0 8])

figure;
plot(data.Rt, data.z,'k');
hold on
plot(data.Rt + data.stdRt, data.z,':k') ;
plot(data.Rt - data.stdRt, data.z,':k') ;
hold off
xlabel('$$\mathrm{R_T}$$','interpreter','latex')
ylabel('z/L')
xlim([0 3000])

