Nlevels_no_dopler_with_z

figure(1)
imagesc(z_1, t_1, I1_out_1); colorbar
xlabel('z')
ylabel('t')
zlabel('I_1')
title('I_1')

figure(2)
imagesc(z_2, t_2, r13_re_out_2); colorbar
xlabel('z')
ylabel('t')
zlabel('Real(rho_{13})')
title('Real(rho_{13})')

figure(3)
imagesc(z_2, t_2, r13_im_out_2); colorbar
xlabel('z')
ylabel('t')
zlabel('Imag(rho_{13})')
title('Imag(rho_{13})')

figure(4)
plot( ...
	t_1,I1_out_1(:,1),'-;before;',  ...
	t_1,I1_out_1(:,end), '-;after;' ...
	)
xlabel('t')
ylabel('I_1')