Added post processing files

1 files changed, 42 insertions, 0 deletions

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')
+