Add plot of all density matrix elements in one plot

1 files changed, 90 insertions, 0 deletions

diff --git a/xmds2/Nlevels_no_dopler_with_z/pp.m b/xmds2/Nlevels_no_dopler_with_z/pp.m
index d33d0e1..09253ac 100644
--- a/xmds2/Nlevels_no_dopler_with_z/pp.m
+++ b/xmds2/Nlevels_no_dopler_with_z/pp.m
@@ -54,6 +54,7 @@ plot( ...
 	)
 xlabel('t')
 ylabel('I_1')
+title('Pulse propagation')
 
 figure(5)
 imagesc(z_2, t_2, r14_re_out_2); colorbar
@@ -69,3 +70,92 @@ ylabel('t')
 zlabel('Imag(rho_{14})')
 title('Imag(rho_{14})')
 
+%% all in one
+% diagonal populations, 
+% upper triangle real part of coherences, 
+% lower diagonal imaginary part of coherences
+figure(22)
+subplot(4,4,1); imagesc (z_2, t_2, r11_out_2);  caxis([0,1]); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('rho_{11}')
+title('rho_{11}')
+subplot(4,4,6); imagesc (z_2, t_2, r22_out_2);  caxis([0,1]); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('rho_{22}')
+title('rho_{22}')
+subplot(4,4,11); imagesc (z_2, t_2, r33_out_2);  caxis([0,1]); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('rho_{33}')
+title('rho_{33}')
+subplot(4,4,16); imagesc (z_2, t_2, r44_out_2); caxis([0,1]);  colorbar
+xlabel('z')
+ylabel('t')
+zlabel('rho_{44}')
+title('rho_{44}')
+% real parts of coherences
+subplot(4,4,2); imagesc(z_2, t_2, r12_re_out_2); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('Real(rho_{12})')
+title('Real(rho_{12})')
+subplot(4,4,3); imagesc(z_2, t_2, r13_re_out_2); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('Real(rho_{13})')
+title('Real(rho_{13})')
+subplot(4,4,4); imagesc(z_2, t_2, r14_re_out_2); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('Real(rho_{14})')
+title('Real(rho_{14})')
+subplot(4,4,7); imagesc(z_2, t_2, r23_re_out_2); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('Real(rho_{23})')
+title('Real(rho_{23})')
+subplot(4,4,8); imagesc(z_2, t_2, r24_re_out_2); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('Real(rho_{24})')
+title('Real(rho_{24})')
+subplot(4,4,12); imagesc(z_2, t_2, r34_re_out_2); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('Real(rho_{34})')
+title('Real(rho_{34})')
+% imaginary parts of coherences
+subplot(4,4,5); imagesc(z_2, t_2, r12_im_out_2); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('Imag(rho_{12})')
+title('Imag(rho_{12})')
+subplot(4,4,9); imagesc(z_2, t_2, r13_im_out_2); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('Imag(rho_{13})')
+title('Imag(rho_{13})')
+subplot(4,4,10); imagesc(z_2, t_2, r23_im_out_2); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('Imag(rho_{23})')
+title('Imag(rho_{23})')
+subplot(4,4,13); imagesc(z_2, t_2, r14_im_out_2); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('Imag(rho_{14})')
+title('Imag(rho_{14})')
+subplot(4,4,14); imagesc(z_2, t_2, r24_im_out_2); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('Imag(rho_{24})')
+title('Imag(rho_{24})')
+subplot(4,4,15); imagesc(z_2, t_2, r34_im_out_2); colorbar
+xlabel('z')
+ylabel('t')
+zlabel('Imag(rho_{34})')
+title('Imag(rho_{34})')
+
+