diff options
| -rw-r--r-- | MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds | 8 |
1 files 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 |