%%%%%%%%            Data for Disks settling in Turbulence JFM
%%%%%%%%            Luis Blay Esteban

%%%%%%%%            The X,Y,Z coordinates are given for all figures
%%%%%%%%            presented in the paper. The framerate was always fixed
%%%%%%%%            at 60fps. Each section has explanation if required.
clear all
load('DataSet.mat')

%% All figures showing particle trajectory can be created as

figure
plot3(xfig1,yfig1,zfig1,'b.')
axis equal
%% Figure 8 shows trajectory but also angle, the angle is saved under the variable alpha
figure
plot3(xfig_ang,yfig_ang,zfig_ang,'b.')
axis equal
figure
plot(real(alpha))

%% All figures showing frequency content can be created as

PSDX_d100=[PSDX_02f_d100t2; PSDX_02f_d100t10; PSDX_02f_d100t20];
Trap_d100=trapz(f(1:50),mean(PSDX_d100(:,1:50),1));

figure
plot(f(1:50),mean(PSDX_d100(:,1:50),1)./Trap_d100,'r:','LineWidth',2)


%% For Figure 10
%variable mean_vz and std_mean: %first row quiescent, second dt=2,
%third dt=10, forth dt=20; first columns: disk1, second disk2, third disk3;

%variable Disk_variables: first row Archimedes, second Inertia, third
%turbulence level; first columns: disk1, second disk2, third disk3;

figure
semilogx(Disk_variables(3,:),(mean_vz(2:4,1)./mean_vz(1,1)),'ko:','LineWidth',2)%./Disk_variables(2,1))%Disk1
hold on
semilogx(Disk_variables(3,:),(mean_vz(2:4,2)./mean_vz(1,2)),'ko-.','LineWidth',2)%./Disk_variables(2,2))%Disk2
semilogx(Disk_variables(3,:),(mean_vz(2:4,3)./mean_vz(1,3)),'ko-','LineWidth',2)%./Disk_variables(2,3))%Disk3
ylabel('$\langle V_{z_t} \rangle / \langle V_{z} \rangle $','Interpreter','Latex','FontSize',14)
xlabel('$q^2[m^2s^2]$','Interpreter','Latex','FontSize',14)

%% Figure 14 can be created as

figure
p1=plot(fftmax_d100_quies./mean(fftmax_d100_quies), (V_z_d100_quies./mean(V_z_d100_quies)),'ko','MarkerFaceColor','k')
axis([0 1.4 0.6 2.6])
xlabel('$f/f_q$','Interpreter','Latex')
ylabel('$V_z/V_{z_q}$','Interpreter','Latex')