diff options
| -rw-r--r-- | xmds2/realistic_Rb_and_fields/pp_I2.m | 82 | ||||
| -rw-r--r-- | xmds2/realistic_Rb_and_fields/two_field_circular_eit_cold.params | 10 |
2 files changed, 92 insertions, 0 deletions
diff --git a/xmds2/realistic_Rb_and_fields/pp_I2.m b/xmds2/realistic_Rb_and_fields/pp_I2.m new file mode 100644 index 0000000..c3c6932 --- /dev/null +++ b/xmds2/realistic_Rb_and_fields/pp_I2.m @@ -0,0 +1,82 @@ +% suppress output of plots on display +set (0, "defaultfigurevisible", "off") + +realistic_Rb_and_fields + +%% field I2 propagation +z_1=z_1*100; % z in cm +t_1=t_1*1e6; % time now measured in uS +figure(1) +%set(gca,'fontsize',20); +imagesc(z_1, t_1, Ip2_out_1); colorbar +tmin=-0.4; +tmax= 0.4; +ylim([tmin,tmax],'manual'); +xlabel('z (cm)') +ylabel('t (uS)') +zlabel('I_p2') +title('I_p2') +print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_Ip2.eps') + +figure(1) +%set(gca,'fontsize',20); +imagesc(z_1, t_1, Im2_out_1); colorbar +tmin=-0.4; +tmax= 0.4; +ylim([tmin,tmax],'manual'); +xlabel('z (cm)') +ylabel('t (uS)') +zlabel('I_m2') +title('I_m2') +print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_Im2.eps') + + +%% field I4 propagation +figure(3) + +imagesc(z_1, t_1, Ip4_out_1); colorbar +tmin=-0.4; +tmax= 0.4; +ylim([tmin,tmax],'manual'); +xlabel('z (cm)') +ylabel('t (uS)') +zlabel('I_p4') +title('I_p4') +print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_Ip4.eps') + +figure(4) + +imagesc(z_1, t_1, Im4_out_1); colorbar +tmin=-0.4; +tmax= 0.4; +ylim([tmin,tmax],'manual'); +xlabel('z (cm)') +ylabel('t (uS)') +zlabel('I_m4') +title('I_m4') +print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_Im4.eps') + + + + + +%% fields before and after the cell +figure(5) +%set(gca,'fontsize',30); +plot( ... + t_1,Ip2_out_1(:,1),'.-;I_p2 before;', "linewidth", 4 ... + ,t_1,Im2_out_1(:,1),'.-;I_m2 before;', "linewidth", 4 ... + ,t_1,Ip2_out_1(:,end), '-;I_p2 after;', "linewidth", 4 ... + ,t_1,Im2_out_1(:,end), '-;I_m2 after;', "linewidth", 4 ... + ,t_1,Ip4_out_1(:,end), '-;I_p4 after;', "linewidth", 4 ... + ,t_1,Im4_out_1(:,end), '-;I_m after;', "linewidth", 4 ... + ) +xlabel('t (uS)') +ylabel('I (1/s)^2') +title('Fields before and after cell') +legend('location', 'northeast'); +print('-color','-depsc2', 'fields_before_after_cell.eps') + + +fid_stamp=fopen ("pp_I2.stamp", "a+"); +fclose (fid_stamp); diff --git a/xmds2/realistic_Rb_and_fields/two_field_circular_eit_cold.params b/xmds2/realistic_Rb_and_fields/two_field_circular_eit_cold.params new file mode 100644 index 0000000..325be02 --- /dev/null +++ b/xmds2/realistic_Rb_and_fields/two_field_circular_eit_cold.params @@ -0,0 +1,10 @@ + +PARAMS = \ + --Ndens=1e15 \ + --Lcell=10.0e-2 \ + --Temperature=1e-3 \ + --WLx=0 --WLy=0 --WLz=0 \ + --Pwidth=0.4e-6 \ + --delta1=0 --delta2=0 --delta3=0 \ + --Ep1o=5e7 --Ep2o=1e2 --Ep3o=0 --Ep4o=0 --Em1o=0 --Em2o=0 --Em3o=0 --Em4o=0 + |