summaryrefslogtreecommitdiff
path: root/xmds2
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2012-01-12 15:30:56 -0500
committerEugeniy Mikhailov <evgmik@gmail.com>2012-01-12 15:34:09 -0500
commit36d1b3aa9e3b2335ea7a33ab776f3515083e4aa3 (patch)
tree178c660c450bced6e1bd9004cdafcd892dd71a2a /xmds2
parent416702a44bdc76cf49e9cd3b980fdf73f50d8caf (diff)
downloadNresonances-36d1b3aa9e3b2335ea7a33ab776f3515083e4aa3.tar.gz
Nresonances-36d1b3aa9e3b2335ea7a33ab776f3515083e4aa3.zip
Doppler feature added: --Ndens now is taken in account
Diffstat (limited to 'xmds2')
-rw-r--r--xmds2/Nlevels_with_doppler_with_z_4wm/Nlevels_with_doppler_with_z_4wm.xmds8
1 files changed, 7 insertions, 1 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 c63870e..d8fd6a4 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
@@ -50,7 +50,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 sequence -> filter element to use --Ndens parameter
+ 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
@@ -287,6 +288,11 @@
</vector>
<sequence>
+ <filter>
+ <![CDATA[
+ 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)]
+ ]]>
+ </filter>
<!--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="Lcell">