diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2016-05-13 11:51:38 -0400 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2016-05-13 11:51:38 -0400 |
| commit | 25e8ec6e381a1c4896c83ce3d35af7918132b2da (patch) | |
| tree | ef1f05287f0cb09cbeedaca8c5fa70f3a20d876f /MOR_5_levels_with_doppler_and_propagation | |
| parent | a0aeceac697bd1fde96a0f3c3db78d19bbe9c2c4 (diff) | |
| download | noisy_eit_xmds-25e8ec6e381a1c4896c83ce3d35af7918132b2da.tar.gz noisy_eit_xmds-25e8ec6e381a1c4896c83ce3d35af7918132b2da.zip | |
Simulation time span is a program argument now
Diffstat (limited to 'MOR_5_levels_with_doppler_and_propagation')
| -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 |