From 25e8ec6e381a1c4896c83ce3d35af7918132b2da Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Fri, 13 May 2016 11:51:38 -0400
Subject: Simulation time span is a program argument now

---
 MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

diff --git a/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds b/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds
index eac3d26..ece5454 100644
--- a/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds
+++ b/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds
@@ -84,14 +84,16 @@
 				double probability_v; // will be used as p(v) in Maxwell distribution
 
 				//  ---------  Simulation parameters  -------------------------
-				const double t_span = 2e-5; // span of the simulation time lattice
 			]]>
 		</globals>
 		<validation kind="run-time"/> <!--allows to put ranges as variables-->
 		<benchmark />
     <arguments>
 			<!-- ============== Calculation parameters =================== -->
+			<!--Note timeSpan/numTimePoints should be less than 1/(highest frequency or detuning) in simulation-->
+			<!--May be some heuristic is needed-->
       <argument name="numTimePoints" type="integer" default_value="10000" />
+      <argument name="timeSpan" type="real" default_value="2e-5" />
 
 			<!-- Rabi frequency divided by 2 in [1/s] -->
       <argument name="EdLo" type="real" default_value="2*1.5*(2*M_PI*1e6)" />
@@ -163,7 +165,7 @@
 						* lattice="10000"   domain="(-1e-6, 1e-6)" 
 							works for Doppler averaging in up to 400K for Rb when lasers are zero detuned
 			 -->
-			<dimension name="t"   lattice="numTimePoints"   domain="(-t_span/2, t_span/2)" />
+			<dimension name="t"   lattice="numTimePoints"   domain="(-timeSpan/2, timeSpan/2)" />
 			<dimension name="v"   lattice="3"   domain="(V_maxwell_min, V_maxwell_max)" />
 		</transverse_dimensions>
 	</geometry>
@@ -324,7 +326,7 @@
 
 					double delta_one_ph = delta_one_ph0;
 					double delta_two_ph = delta_two_ph0;
-					delta_two_ph += abs(Gamma_super)*1*((t)/(t_span/2)); // two photon sweep to mimic EIT shapes
+					delta_two_ph += abs(Gamma_super)*1*((t)/(timeSpan/2)); // two photon sweep to mimic EIT shapes
 
 					Gab=Gamma_super+i*( Split0 + delta_one_ph + delta_two_ph + DopplerShiftPr); // probe right detuning
 					GAB=Gamma_super+i*(-Split0 + delta_one_ph + delta_two_ph + DopplerShiftPr); // probe left detuning
-- 
cgit v1.2.3