Merge branch 'master' of qo.physics.wm.edu:Nresonances

16 files changed, 175 insertions, 101 deletions

diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm/Makefile
index 578f2b0..f8b3559 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm/Makefile
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/Makefile
@@ -15,7 +15,7 @@ M_FILES = $(patsubst %.xmds,%.m,$(XMDS_FILES))
 XMDS = xmds2
 XSIL2GRAPHICS = xsil2graphics
 
-all: $(M_FILES) 
+all: $(RUN_FILES)
 
 %.run: %.xmds
 	$(XMDS) $<
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds b/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds
index 44b13f7..bb0feb2 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.xmds
@@ -36,7 +36,7 @@
 
 		VERY IMPORTANT: all Rabi frequency should be given in [1/s], if you want to
 		normalize it to something else look drho/dt equation.
-		No need to renormalizes eta as long as its express through i
+		No need to renormalizes eta as long as its express through
 		the upper level decay rate in the same units as Rabi frequency.
 	</description>
 
@@ -49,7 +49,8 @@
 				const double lambda=794.7e-9; //wavelength in m
 				const double Kvec = 2*M_PI/lambda;  // k-vector
 				const double Gamma_super=6*(2*M_PI*1e6);  // characteristic decay rate of upper level used for eta  calculations expressed in [1/s]
-				const double eta = 3*lambda*lambda*Ndens*Gamma_super/8.0/M_PI;  // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)]
+				// eta will be calculated in the <arguments> section
+				double eta = 0;  // eta constant in the wave equation for Rabi frequency. Units are [1/(m s)]
 				
 				//  ---------  Atom and cell properties -------------------------
 				// range of Maxwell distribution atomic velocities
@@ -57,9 +58,10 @@
 				// above mass expression is written as (expression is isotopic_mass * atomic_mass_unit)
 
 				// Average sqrt(v^2) in Maxwell distribution for one dimension
-				const double v_thermal_averaged=sqrt(k_boltzmann*Temperature/mass); 
+				// Maxwell related parameters will be calculated in <arguments> section
+				double v_thermal_averaged=0;
 				// Maxwell distribution velocities range to take in account in [m/s]
-				const double V_maxwell_min = -4*v_thermal_averaged, V_maxwell_max = -V_maxwell_min; // there is almost zero probability for higher velocity p(4*v_av) = 3.3e-04 * p(0)
+				double V_maxwell_min = 0, V_maxwell_max = 0;
 
 				// repopulation rate (atoms flying in/out the laser beam)  in [1/s]
 				const double gt=0.01 *(2*M_PI*1e6);
@@ -101,6 +103,18 @@
 					 I am guessing detunings are too large and thus it became a stiff equation-->
 			<!--! make sure it is not equal to zero!-->
 			<argument name="Temperature"  type="real" default_value="5" />
+			<!-- This will be executed after arguments/parameters are parsed -->
+			<!-- Read the code Luke: took me a while of reading the xmds2 sources to find it -->
+			<![CDATA[
+				// Average sqrt(v^2) in Maxwell distribution for one dimension
+				v_thermal_averaged=sqrt(k_boltzmann*Temperature/mass); 
+				// Maxwell distribution velocities range to take in account in [m/s]
+				// there is almost zero probability for higher velocity p(4*v_av) = 3.3e-04 * p(0)
+				V_maxwell_min = -4*v_thermal_averaged; V_maxwell_max = -V_maxwell_min; 
+
+				// eta constant in the wave equation for Rabi frequency. Units are [1/(m s)]
+				eta = 3*lambda*lambda*Ndens*Gamma_super/8.0/M_PI;
+			]]>
     </arguments>
 		<bing />
 		<fftw plan="patient" />
@@ -112,7 +126,7 @@
 	<geometry>
 		<propagation_dimension> z </propagation_dimension>
 		<transverse_dimensions>
-			<dimension name="t"   lattice="1000"   domain="(-14.0e-7, 4.0e-7)" />
+			<dimension name="t"   lattice="1000"   domain="(-0.2e-6, 0.4e-6)" />
 		</transverse_dimensions>
 	</geometry>
 
@@ -310,7 +324,7 @@
 		</group>
 
 		<group>
-      <sampling basis="t(100)" initial_sample="yes">
+      <sampling basis="t(1000)" initial_sample="yes">
 				<dependencies>density_matrix</dependencies>
 				<moments>
 					r11_out r22_out r33_out r44_out 
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_light/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_light/plot_fields_propagation_I4.gp
new file mode 120000
index 0000000..134a380
--- /dev/null
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_light/plot_fields_propagation_I4.gp
@@ -0,0 +1 @@
+../plot_fields_propagation_I4.gp
\ No newline at end of file
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile
index 697fd27..fcc817c 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/Makefile
@@ -9,11 +9,17 @@ GNUPLOT_FILES = $(wildcard *.gp)
 XSIL2GRAPHICS = xsil2graphics
 
 # fast light
-# PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4
+PARAMS =  \
+	  --Ndens=1e15 \
+	  --Lcell=10.0e-2 \
+	  --Temperature=1e-9 \
+	  --Pwidth=0.1e-6 \
+	  --delta1=0 --delta2=0 --delta3=0 \
+	  --E1o=1.9e7 --E2o=1e2 --E3o=3.8e7 --E4o=0
+
+
 # slow light EIT
 #PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0
-#Fast light to Slow light switch
-PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=2e7 --E2o=3e3 --E3o=6e6 --E4o=3e3
 
 all: $(XSIL_FILES) Nlevels_no_dopler_with_z_4wm.xsil $(M_FILES) plot png
 
@@ -23,9 +29,12 @@ Nlevels_no_dopler_with_z_4wm.xsil: ../Nlevels_no_dopler_with_z_4wm.run
 %.m: %.xsil
 	$(XSIL2GRAPHICS) $<
 
-plot: $(M_FILES) $(GNUPLOT_FILES) 
-	octave pp_I2.m
+pretty_plots: plot $(GNUPLOT_FILES)
 	gnuplot plot_fields_propagation_I2.gp
+	gnuplot plot_fields_propagation_I4.gp
+
+plot: $(M_FILES)
+	octave pp_I2.m
 
 clean:
 	rm -f $(CC_FILES) $(RUN_FILES) $(M_FILES) $(XSIL_FILES) *.wisdom.fftw3 *.dat octave-core *.wisdom *.pdf
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/README b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/README
new file mode 100644
index 0000000..ec7d200
--- /dev/null
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/README
@@ -0,0 +1,18 @@
+fast to slow transition
+PARAMS =  \
+	  --Ndens=1e15 \
+	  --Lcell=10.0e-2 \
+	  --Temperature=1e-9 \
+	  --delta1=0 --delta2=0 --delta3=0 \
+	  --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0 
+
+1.5 cm seems to give largest advance
+2.5 cm gives zero delay
+10cm slow light again
+
+fast light with decay in 1.5 cm cell
+PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=0
+
+
+PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=3.0e7 --E2o=3.1e5 --E3o=3.0e7 --E4o=6.3e3
+10 cm first fast light  then around 8 cm transition to slow light
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/plot_fields_propagation_I4.gp
new file mode 120000
index 0000000..134a380
--- /dev/null
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/plot_fields_propagation_I4.gp
@@ -0,0 +1 @@
+../plot_fields_propagation_I4.gp
\ No newline at end of file
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/pp_I2.m
index 73487e0..4a1ed9b 100644..120000
--- a/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/pp_I2.m
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/fast_to_slow_switch/pp_I2.m
@@ -1,67 +1 @@
-Nlevels_no_dopler_with_z_4wm
-
-%% field 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, I2_out_1); colorbar
-xlabel('z (cm)')
-ylabel('t (uS)')
-zlabel('I_2')
-title('I_2')
-
-xskip=1;
-yskip=10;
-map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip);
-
-
-
-print('-color','fields_propagation_I2.eps')
-
-
-
-%% fields before and after the cell
-figure(2)
-%set(gca,'fontsize',30); 
-plot( ...
-	t_1,I2_out_1(:,1),'.-;before;', "linewidth", 4,  ...
-	t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ...
-	)
-xlabel('t (uS)')
-ylabel('I_2 (1/s)^2')
-title('I_2 before and after cell')
-legend('location', 'southwest');
-
-[b, a]=butter(3, 0.05);
-I2_out_after=I2_out_1(:,end);
-I2_out_after_filtered=filtfilt(b,a,I2_out_after);
-settling_time=0.8; %uS
-t_good_indx=t_1> min(t_1 + settling_time);
-[m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx));
-delay_time=t_1(max_pos_after)-t_1(max_pos_before);
-printf('Second field delay time  = %f uS\n',delay_time);
-
-%set(gca,'fontsize',40); 
-set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts
-print('-color','fields_before_after_cell_I2.eps')
-
-figure(4)
-I2_max_in=max(I2_out_1(t_good_indx,1));
-I2_max_out=max(I2_out_1(t_good_indx,end));
-I2_in_norm=(I2_out_1(:,1))/I2_max_in;
-I2_out_norm=(I2_out_1(:,end))/I2_max_out;
-tmin=-0.05;
-tmax=0.05;
-indx=(t_1>=tmin & t_1<=tmax);  % soom in in time to this region
-plot( ...
-	t_1(indx),I2_in_norm(indx),'.-;before;', "linewidth", 4,  ...
-	t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ...
-	)
-xlim([tmin,tmax],'manual');	
-xlabel('t (uS)')
-ylabel('I_2')
-title('I_2 before and after cell normalized')
-set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts
-print('-color','probe_before_after_cell_I2_normalized.eps')
-
+../pp_I2.m
\ No newline at end of file
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I2.gp
index b6721d4..77b4870 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I2.gp
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I2.gp
@@ -12,4 +12,4 @@ set ylabel "t ({/Symbol m}S)"
 set zlabel "I_2 (1/S)"
 set nokey
 #set view map
-splot [0:1.5][-0.4:0.4] 'I2.dat' 
+splot [0:][-0.2:0.2] 'I2.dat' 
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I4.gp
new file mode 100644
index 0000000..3849e2f
--- /dev/null
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/plot_fields_propagation_I4.gp
@@ -0,0 +1,15 @@
+set terminal postscript portrait enhanced color solid size 5,3.5
+set output 'fields_propagation_I4.eps'
+set dgrid3d 100,100 qnorm 4
+set pm3d map
+#set contour
+set hidden3d
+set palette rgb 10,13,31 negative
+
+
+set xlabel "z (cm)"
+set ylabel "t ({/Symbol m}S)"
+set zlabel "I_2 (1/S)"
+set nokey
+#set view map
+splot [0:][-0.2:0.2] 'I4.dat' 
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m
index 45d4913..3e5c174 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pp_I2.m
@@ -1,6 +1,6 @@
-Nlevels_no_dopler_with_z_4wm
+Nlevels_no_dopler_with_z_4wm;
 
-%% field propagation
+%% field I2  propagation
 z_1=z_1*100; % z in cm
 t_1=t_1*1e6; % time now measured in uS
 figure(1)
@@ -13,33 +13,50 @@ xlabel('z (cm)')
 ylabel('t (uS)')
 zlabel('I_2')
 title('I_2')
+print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps')
 
-xskip=1;
-yskip=10;
-%map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip);
+desired_x_size=200;
+desired_y_size=200;
+xskip=ceil(length(z_1)/desired_x_size);
+yskip=ceil(length(t_1)/desired_y_size);
+map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip);
 
 
+%% field I4  propagation
+figure(2)
+
+imagesc(z_1, t_1, I4_out_1); colorbar
+tmin=-0.4;
+tmax= 0.4;
+ylim([tmin,tmax],'manual');	
+xlabel('z (cm)')
+ylabel('t (uS)')
+zlabel('I_4')
+title('I_4')
+print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I4.eps')
+map2dat('I4.dat',z_1,t_1, I4_out_1, xskip, yskip);
+
 
-print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps')
 
 
 
 %% fields before and after the cell
-figure(2)
+figure(3)
 %set(gca,'fontsize',30); 
 plot( ...
-	t_1,I2_out_1(:,1),'.-;before;', "linewidth", 4,  ...
-	t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ...
+	t_1,I2_out_1(:,1),'.-;I_2 before;', "linewidth", 4  ...
+	,t_1,I2_out_1(:,end), '-;I_2 after;', "linewidth", 4 ...
+	,t_1,I4_out_1(:,end), '-;I_4 after;', "linewidth", 4 ...
 	)
 xlabel('t (uS)')
-ylabel('I_2 (1/s)^2')
-title('I_2 before and after cell')
-legend('location', 'southwest');
+ylabel('I (1/s)^2')
+title('Fields before and after cell')
+legend('location', 'northeast');
 
 [b, a]=butter(3, 0.05);
 I2_out_after=I2_out_1(:,end);
 I2_out_after_filtered=filtfilt(b,a,I2_out_after);
-settling_time=0.8; %uS
+settling_time=0.01; %uS
 t_good_indx=t_1> min(t_1 + settling_time);
 [m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx));
 delay_time=t_1(max_pos_after)-t_1(max_pos_before);
@@ -47,25 +64,29 @@ printf('Second field delay time  = %f uS\n',delay_time);
 
 %set(gca,'fontsize',40); 
 %set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts
-print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell_I2.eps')
+print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell.eps')
 
 figure(4)
 I2_max_in=max(I2_out_1(t_good_indx,1));
 I2_max_out=max(I2_out_1(t_good_indx,end));
+I4_max_out=max(I4_out_1(t_good_indx,end));
 I2_in_norm=(I2_out_1(:,1))/I2_max_in;
 I2_out_norm=(I2_out_1(:,end))/I2_max_out;
-tmin=-0.05;
-tmax=0.05;
+I4_out_norm=(I4_out_1(:,end))/I4_max_out;
+tmin=-.05;
+tmax=.05;
 indx=(t_1>=tmin & t_1<=tmax);  % soom in in time to this region
 plot( ...
-	t_1(indx),I2_in_norm(indx),'.-;before;', "linewidth", 4,  ...
-	t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ...
+	t_1(indx),I2_in_norm(indx),'.-;I_2 before;', "linewidth", 4,  ...
+	t_1(indx),I2_out_norm(indx), '-;I_2 after;', "linewidth", 4 ...
+	,t_1(indx),I4_out_norm(indx), '-;I_4 after;', "linewidth", 4 ...
 	)
 legend('location', 'southeast');
 xlim([tmin,tmax],'manual');	
+ylim([0.5,1],'manual');	
 xlabel('t (uS)')
-ylabel('I_2')
-title('I_2 before and after cell normalized')
+ylabel('I normalized')
+title('Fields before and after cell normalized')
 %set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts
-print('-color','-depsc2', '-tight','-S200,120', 'probe_before_after_cell_I2_normalized.eps')
+print('-color','-depsc2', '-tight','-S200,120', 'fields_before_after_cell_normalized.eps')
 
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/Makefile
new file mode 100644
index 0000000..95ebb0a
--- /dev/null
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/Makefile
@@ -0,0 +1,56 @@
+###  -*- make -*-
+### This makefile can be used to build and run the XMDS examples
+
+
+XSIL_FILES = Nlevels_no_dopler_with_z_4wm.xsil
+M_FILES = $(patsubst %.xsil,%.m,$(XSIL_FILES))
+GNUPLOT_FILES = $(wildcard *.gp)
+
+XSIL2GRAPHICS = xsil2graphics
+
+# fast light
+# PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4
+# slow light EIT
+#PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0
+#Fast light to Slow light switch
+PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=2e7 --E2o=3e3 --E3o=6e6 --E4o=3e2
+
+all: $(XSIL_FILES) Nlevels_no_dopler_with_z_4wm.xsil $(M_FILES) plot png
+
+Nlevels_no_dopler_with_z_4wm.xsil: ../Nlevels_no_dopler_with_z_4wm.run
+	$< $(PARAMS) | grep "Time elapsed for simulation is:" > exact_analysis_execution_time.txt
+
+%.m: %.xsil
+	$(XSIL2GRAPHICS) $<
+
+plot: $(M_FILES) $(GNUPLOT_FILES) 
+	octave pp_I2.m
+	gnuplot plot_fields_propagation_I2.gp
+
+clean:
+	rm -f $(CC_FILES) $(RUN_FILES) $(M_FILES) $(XSIL_FILES) *.wisdom.fftw3 *.dat octave-core *.wisdom *.pdf
+	rm -f $(pdf_targets)
+	rm -f $(eps_targets)
+
+real_clean: clean
+	rm -f $(png_targets)
+
+eps_targets  = $(wildcard *.eps)
+pdf_targets  = $(eps_targets:%.eps=%.pdf)	
+png_targets  = $(pdf_targets:%.pdf=%.png)
+
+pdf: $(pdf_targets) 
+
+$(pdf_targets): %.pdf : %.eps
+	cat $< | ps2eps -B > __tt.eps
+	epspdf __tt.eps $@
+	rm -f __tt.eps
+	#ps2eps -B $< | epspdf $< $@
+
+png: pdf $(png_targets) 
+
+$(png_targets): %.png : %.pdf
+	convert -density 300 $< $@
+
+.PRECIOUS: %.run %.xsil %.m
+.PHONY: all clean
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/map2dat.m b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/map2dat.m
new file mode 120000
index 0000000..14fae30
--- /dev/null
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/map2dat.m
@@ -0,0 +1 @@
+../map2dat.m
\ No newline at end of file
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp
new file mode 120000
index 0000000..e11c569
--- /dev/null
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp
@@ -0,0 +1 @@
+../plot_fields_propagation_I2.gp
\ No newline at end of file
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I4.gp
new file mode 120000
index 0000000..134a380
--- /dev/null
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I4.gp
@@ -0,0 +1 @@
+../plot_fields_propagation_I4.gp
\ No newline at end of file
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m
new file mode 120000
index 0000000..4a1ed9b
--- /dev/null
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m
@@ -0,0 +1 @@
+../pp_I2.m
\ No newline at end of file
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/slow_light/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/slow_light/plot_fields_propagation_I4.gp
new file mode 120000
index 0000000..134a380
--- /dev/null
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/slow_light/plot_fields_propagation_I4.gp
@@ -0,0 +1 @@
+../plot_fields_propagation_I4.gp
\ No newline at end of file