added post processing script and sample case parameters

1 files changed, 82 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);