3 files changed, 148 insertions, 0 deletions
diff --git a/xmds2/Genas_system/pp.m b/xmds2/Genas_system/pp.m
new file mode 100644
index 0000000..3762fcb
--- /dev/null
+++ b/xmds2/Genas_system/pp.m
@@ -0,0 +1,14 @@
+Genas_system
+
+%% field propagation
+pp_fields
+
+
+
+
+
+%% all density matrix elements in one plot
+% diagonal populations, 
+% upper triangle real part of coherences, 
+% lower diagonal imaginary part of coherences
+pp_rho
diff --git a/xmds2/Genas_system/pp_fields.m b/xmds2/Genas_system/pp_fields.m
new file mode 100644
index 0000000..c87f73c
--- /dev/null
+++ b/xmds2/Genas_system/pp_fields.m
@@ -0,0 +1,42 @@
+Genas_system
+
+Ex=Ex_re_out_1 + 1i*Ex_im_out_1;
+Ey=Ey_re_out_1 + 1i*Ey_im_out_1;
+
+Ecw= Ex+1i*Ey;
+Eccw=Ex-1i*Ey;
+
+%% field propagation
+t_1=t_1*1e9; % time now measured in uS
+figure(1)
+subplot(1,2,1); imagesc(z_1, t_1, abs(Ecw)); colorbar
+xlabel('z')
+ylabel('t (nS)')
+zlabel('|E_{cw}|')
+title('|E_{cw}|')
+subplot(1,2,2); imagesc(z_1, t_1, abs(Eccw)); colorbar
+xlabel('z')
+ylabel('t (nS)')
+zlabel('|E_{ccw}|')
+title('|E_{ccw}|')
+
+print('fields_map.pdf')
+
+% fields before and after
+[Nr,Nc] = size( Ecw );
+Ecw_b=Ecw(:,1);
+Ecw_a=Ecw(:,Nc);
+Eccw_b=Eccw(:,1);
+Eccw_a=Eccw(:,Nc);
+figure(2)
+subplot(2,1,1);
+plot(t_1, abs(Ecw_b), '-;E_{cw_{before}};', t_1, abs(Ecw_a), '-;E_{cw_{after}};'); 
+xlabel('t (nS)')
+ylabel('Rabi frequency (1/s)')
+subplot(2,1,2);
+plot(t_1, abs(Eccw_b), '-;E_{ccw_{before}};', t_1, abs(Eccw_a), '-;E_{ccw_{after}};')
+xlabel('t (nS)')
+ylabel('Rabi frequency (1/s)')
+
+print('fields_before_after.pdf')
+
diff --git a/xmds2/Genas_system/pp_rho.m b/xmds2/Genas_system/pp_rho.m
new file mode 100644
index 0000000..515f746
--- /dev/null
+++ b/xmds2/Genas_system/pp_rho.m
@@ -0,0 +1,92 @@
+Genas_system
+
+%% all density matrix elements in one plot
+% diagonal populations, 
+% upper triangle real part of coherences, 
+% lower diagonal imaginary part of coherences
+figure(3)
+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})')
+
+print('rho_map.pdf')
+