diff options
| -rw-r--r-- | xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds b/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds index aff6c02..1a43fd9 100644 --- a/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds +++ b/xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds @@ -55,7 +55,7 @@ // --------- Atom and cell properties ------------------------- // range of Maxwell distribution atomic velocities - const double Temperature=1.0; // cell temperature in [K] ! make sure it is not equal to zero! + const double Temperature=1.0e-3; // cell temperature in [K] ! make sure it is not equal to zero! // TODO: looks like Temperature > 10 K knocks solver, I am guessing detunings are too large and thus it became a stiff equation const double mass = (86.909180527 * 1.660538921e-27); // atom mass in [kg] // above mass expression is written as (expression is isotopic_mass * atomic_mass_unit) @@ -280,7 +280,7 @@ <sequence> <!--For this set of conditions ARK45 is faster than ARK89--> <!--ARK45 is good for small detuning when all frequency like term are close to zero--> - <integrate algorithm="ARK45" tolerance="1e-5" interval="10.0e-2"> + <integrate algorithm="ARK45" tolerance="1e-5" interval="20.0e-2"> <!--RK4 is good for large detunings when frequency like term are big, it does not try to be too smart about adaptive step which ARK seems to make too small--> <!--When ARK45 works it about 3 times faster then RK4 with 1000 steps--> <!--<integrate algorithm="RK4" steps="100" interval="1.5e-2">--> |