Add 3 field light propagation and output to input comparison in one plot

1 files changed, 31 insertions, 52 deletions

diff --git a/xmds2/Nlevels_no_dopler_with_z/pp.m b/xmds2/Nlevels_no_dopler_with_z/pp.m
index 09253ac..bdfff57 100644
--- a/xmds2/Nlevels_no_dopler_with_z/pp.m
+++ b/xmds2/Nlevels_no_dopler_with_z/pp.m
@@ -1,80 +1,59 @@
 Nlevels_no_dopler_with_z
 
+%% field propagation
 figure(1)
-imagesc(z_1, t_1, I1_out_1); colorbar
+subplot(1,3,1); imagesc(z_1, t_1, I1_out_1); colorbar
 xlabel('z')
 ylabel('t')
 zlabel('I_1')
 title('I_1')
-
-
-
-%% populations plot
-figure(20)
-subplot(2,2,1); imagesc (z_2, t_2, r11_out_2);  caxis([0,1]); colorbar
-xlabel('z')
-ylabel('t')
-zlabel('rho_{11}')
-title('rho_{11}')
-subplot(2,2,2); imagesc (z_2, t_2, r22_out_2);  caxis([0,1]); colorbar
-xlabel('z')
-ylabel('t')
-zlabel('rho_{22}')
-title('rho_{22}')
-subplot(2,2,3); imagesc (z_2, t_2, r33_out_2);  caxis([0,1]); colorbar
+subplot(1,3,2); imagesc(z_1, t_1, I2_out_1); colorbar
 xlabel('z')
 ylabel('t')
-zlabel('rho_{33}')
-title('rho_{33}')
-subplot(2,2,4); imagesc (z_2, t_2, r44_out_2); caxis([0,1]);  colorbar
+zlabel('I_2')
+title('I_2')
+subplot(1,3,3); imagesc(z_1, t_1, I3_out_1); colorbar
 xlabel('z')
 ylabel('t')
-zlabel('rho_{44}')
-title('rho_{44}')
+zlabel('I_3')
+title('I_3')
 
-%% coherences plot
-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)
+
+%% fields before and after the cell
+figure(2)
+subplot(1,3,1);
 plot( ...
 	t_1,I1_out_1(:,1),'-;before;',  ...
 	t_1,I1_out_1(:,end), '-;after;' ...
 	)
 xlabel('t')
 ylabel('I_1')
-title('Pulse propagation')
-
-figure(5)
-imagesc(z_2, t_2, r14_re_out_2); colorbar
-xlabel('z')
-ylabel('t')
-zlabel('Real(rho_{14})')
-title('Real(rho_{14})')
+title('I_1 propagation')
+subplot(1,3,2);
+plot( ...
+	t_1,I2_out_1(:,1),'-;before;',  ...
+	t_1,I2_out_1(:,end), '-;after;' ...
+	)
+xlabel('t')
+ylabel('I_2')
+title('I_2 propagation')
+subplot(1,3,3);
+plot( ...
+	t_1,I3_out_1(:,1),'-;before;',  ...
+	t_1,I3_out_1(:,end), '-;after;' ...
+	)
+xlabel('t')
+ylabel('I_3')
+title('I_3 propagation')
 
-figure(6)
-imagesc(z_2, t_2, r14_im_out_2); colorbar
-xlabel('z')
-ylabel('t')
-zlabel('Imag(rho_{14})')
-title('Imag(rho_{14})')
 
-%% all in one
+%% all density matrix elements in one plot
 % diagonal populations, 
 % upper triangle real part of coherences, 
 % lower diagonal imaginary part of coherences
-figure(22)
+figure(3)
 subplot(4,4,1); imagesc (z_2, t_2, r11_out_2);  caxis([0,1]); colorbar
 xlabel('z')
 ylabel('t')