Nlevels_with_doppler_with_z_4wm

%% field propagation
z_3=z_3*100; % z in cm
z_1=z_1*100; % z in cm
t_3=t_3*1e6; % time now measured in uS
t_1=t_1*1e6; % time now measured in uS

[Nv, Nt, Nz]=size(I2_out_v_3);

Npl=Nv/2+1;
Iv_out_v = reshape(I2_out_v_3(1:Npl,:,Nz),Npl,Nt);

for i=1:Npl
	l_str{i}=num2str(v_3(i));
end


figure(1); hold off;
plot(t_3, Iv_out_v'); hold on

plot(t_1, I2_out_1(:,end), '-b', 'LineWidth',2);
l_str{Npl+1}=num2str('Doppler averaged out');

plot(t_1, I2_out_1(:,1), '-r', 'LineWidth',2);
l_str{Npl+2}=num2str('Doppler averaged in');

legend(l_str);

xlabel('time (uS)');

figure(2); hold off;
imagesc(z_1, t_1, I2_out_1)
ylabel('time (uS)');
xlabel('z (cm)');