%% Some experiment constants.
dia_thickfibers = 511e-6; % In meters.
dia_thinfibers = 127e-6; % In meters.
measured_vol_bigdrop = 2*1.5; % microliters.
air_speed_20Hz = 2.429; % m/s, from Martin's PIV.
air_speed_22Hz = 0.0995*22 + 0.4749; % m/s, interpolated from Martin's PIV.
air_speed_24Hz = 0.0995*24 + 0.4749; % m/s, interpolated from Martin's PIV.
air_speed_18Hz = 0.0995*18 + 0.4749; % m/s, interpolated from Martin's PIV.
framespersec = 30; % More or less, really.
r_smalldrop = 1.4/2000; % in meters.  % Very sneaky error! 1.4 mm is a DIA.
pixels_per_cm = 791.6882;
pixels_per_m = pixels_per_cm*100;

%% Some material constants.
oil_viscosity = 0.0046; % Viscosity of oil in Pa.s. (Motile drop.)
air_viscosity = 18.0e-6; % Viscosity of air in Pa.s.
air_density = 1.225; % Density of air in kg/m^3.


%% Include the data from 5May2021:
% 20 Hz controller data.  Small drops.
speed_5May = 14.9086;  % pixels per frame
stdev_5May = 4.7291; % pixels per frame
pixels_per_cm = 791.6882;

speed_SI_5May = (framespersec*speed_5May)/(pixels_per_cm*100)
stdev_SI_5May = (framespersec*stdev_5May)/(pixels_per_m)
reynolds_5May = (air_speed_20Hz*2*r_smalldrop*air_density)/air_viscosity

%% Include data from 28-29APR2021
% -- New zither.
% -- Thick fibers.
% -- 20, 22, and 24 Hz.

slopes_24Hz = [];
cd('ThickFibers/24 Hz')
folder = dir('Images_DSC_*');
for f = 1:length(folder)
    cd(folder(f).name)
    load('slope_rsquared.mat')
    slopes_24Hz = [slopes_24Hz; slopes];
    cd ..
end
cd ..
cd ..
meanall_24Hz = mean(abs(slopes_24Hz));
stdall_24Hz = std(abs(slopes_24Hz));
% This data is in pixels per frame.

speed_SI_24Hz = (framespersec*meanall_24Hz)/(pixels_per_m)
stdev_SI_24Hz = (framespersec*stdall_24Hz)/(pixels_per_m)
reynolds_24Hz = (air_speed_24Hz*(2*r_bigdrop)*air_density)/air_viscosity

clear slopes
slopes_22Hz = [];
cd('ThickFibers/22 Hz')
folder = dir('Images_DSC_*');
for f = 1:length(folder)
    cd(folder(f).name)
    load('slope_rsquared.mat')
    slopes_22Hz = [slopes_22Hz; slopes];
    cd ..
end
cd ..
cd ..
meanall_22Hz = mean(abs(slopes_22Hz));
stdall_22Hz = std(abs(slopes_22Hz));
% This data is in pixels per frame.

speed_SI_22Hz = (framespersec*meanall_22Hz)/(pixels_per_m)
stdev_SI_22Hz = (framespersec*stdall_22Hz)/(pixels_per_m)
reynolds_22Hz = (air_speed_22Hz*(2*r_bigdrop)*air_density)/air_viscosity

clear slopes
slopes_20Hz = [];
cd('ThickFibers/20 Hz')
folder = dir('Images_*');
for f = 1:length(folder)
    cd(folder(f).name)
    load('slope_rsquared.mat')
    slopes_20Hz = [slopes_20Hz; slopes];
    cd ..
end
cd ..
cd ..
meanall_20Hz = mean(abs(slopes_20Hz));
stdall_20Hz = std(abs(slopes_20Hz));
% This data is in pixels per frame.

speed_SI_20Hz = (framespersec*meanall_20Hz)/(pixels_per_m)
stdev_SI_20Hz = (framespersec*stdall_20Hz)/(pixels_per_m)
reynolds_20Hz = (air_speed_20Hz*((2*r_bigdrop))*air_density)/air_viscosity


clear slopes
slopes_18Hz = [];
cd('ThickFibers/18 Hz')
folder = dir('Images_*');
for f = 1:length(folder)
    cd(folder(f).name)
    load('slope_rsquared.mat')
    slopes_18Hz = [slopes_18Hz; slopes];
    cd ..
end
cd ..

meanall_18Hz = mean(abs(slopes_18Hz));
stdall_18Hz = std(abs(slopes_18Hz));
% This data is in pixels per frame.

speed_SI_18Hz = (framespersec*meanall_18Hz)/(pixels_per_m)
stdev_SI_18Hz = (framespersec*stdall_18Hz)/(pixels_per_m)
reynolds_18Hz = (air_speed_18Hz*(2*r_bigdrop)*air_density)/air_viscosity