1 files changed, 20 insertions, 0 deletions
diff --git a/xmds2/Nlevels_with_doppler_with_z_4wm/pp_doppler.m b/xmds2/Nlevels_with_doppler_with_z_4wm/pp_doppler.m
new file mode 100644
index 0000000..b2b6360
--- /dev/null
+++ b/xmds2/Nlevels_with_doppler_with_z_4wm/pp_doppler.m
@@ -0,0 +1,20 @@
+Nlevels_with_doppler_with_z_4wm
+
+%% field propagation
+z_1=z_1*100; % z in cm
+t_1=t_1*1e6; % time now measured in uS
+
+[Nv, Nt, Nz]=size(I2_out_1);
+
+Iv_out = reshape(I2_out_1(1:Nv/2,:,Nz),Nv/2,Nt);
+
+for i=1:Nv/2
+	l_str{i}=num2str(v_1(i));
+end
+l_str{Nv/2+1}=num2str('Doppler averaged');
+
+
+figure(1); hold off;
+plot(Iv_out'); hold on
+plot(I2_out_avgd_3(:,end), 'LineWidth',2);
+legend(l_str);