From 573a31eabfd5a3ba387bd61e74a83e4b1cb761d5 Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Wed, 11 Jan 2012 23:46:57 -0500
Subject: Doppler: new run time parameters

density of atoms, temperature, pulse width
---
 .../Nlevels_with_doppler_with_z_4wm.xmds            | 21 +++++++++++++--------
 1 file changed, 13 insertions(+), 8 deletions(-)

(limited to 'xmds2')

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 1e5b25c..d63704f 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
@@ -49,14 +49,11 @@
 				const double k_boltzmann= 1.3806505e-23; // Boltzmann knostant in [J/K]
 				const double lambda=794.7e-9; //wavelength in m
 				const double Kvec = 2*M_PI/lambda;  // k-vector
-				const double N=1e09*(1e6);    //number of particles per cubic m i.e. density
 				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*N*Gamma_super/8.0/M_PI;  // eta constant in the wave equation for Rabi frequency. Units are [1/(m 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)]
 				
 				//  ---------  Atom and cell properties -------------------------
 				// range of Maxwell distribution atomic velocities
-				const double Temperature=5; // 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)
 
@@ -85,7 +82,7 @@
 
 			]]>
 		</globals>
-		<validation kind="run-time"/> <!--allows to put ranges as wariables-->
+		<validation kind="run-time"/> <!--allows to put ranges as variables-->
 		<benchmark />
     <arguments>
 			<!-- Rabi frequency divided by 2 in [1/s] -->
@@ -97,10 +94,18 @@
       <argument name="delta1"  type="real" default_value="0.0" />
       <argument name="delta2"  type="real" default_value="0.0" />
       <argument name="delta3"  type="real" default_value="0.0" />
-			<!--Cell length-->
-			<argument name="Lcell"  type="real" default_value="1.5e-2" />
-			<!--Pulse duration/width-->
+			<!--Pulse duration/width [s] -->
 			<argument name="Pwidth"  type="real" default_value="0.1e-6" />
+			<!--  Atom and cell properties -->
+			<!--Cell length [m] -->
+			<argument name="Lcell"  type="real" default_value="1.5e-2" />
+			<!--Density of atoms [1/m^3] -->
+			<argument name="Ndens"  type="real" default_value="1e15" />
+			<!--Atoms temperature [K] -->
+			<!--TODO: looks like Temperature > 10 K knocks solver, 
+					 I am guessing detunings are too large and thus it became a stiff equation-->
+			<!--! make sure it is not equal to zero!-->
+			<argument name="Temperature"  type="real" default_value="5" />
     </arguments>
 		<bing />
 		<fftw plan="patient" />
-- 
cgit v1.2.3