summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds13
1 files changed, 5 insertions, 8 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 99d88a0..6455d22 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
@@ -80,9 +80,6 @@
// total decay of i-th level branching ratios. Rij branching of i-th level to j-th
//const double Rab=0.5, Rac=0.5;
-
- complex EdLac, EdRac, EpLac, EpRac;
-
// inner use variables
double probability_v; // will be used as p(v) in Maxwell distribution
@@ -256,11 +253,6 @@
<!--This sets boundary condition at all times and left border of z (i.e. z=0)-->
<dependencies>E_field_avgd</dependencies>
<![CDATA[
- EdLac = conj(EdLa);
- EdRac = conj(EdRa);
- EpLac = conj(EpLa);
- EpRac = conj(EpRa);
-
rbb = 0.25; rcc = 0.25; rBB = 0.25; rCC = 0.25;
raa = 0;
rcb = 0; rab = 0; rca = 0;
@@ -311,6 +303,11 @@
complex rAC=conj(rCA);
complex rBC=conj(rCB);
+ complex EdLac = conj(EdLa);
+ complex EdRac = conj(EdRa);
+ complex EpLac = conj(EpLa);
+ complex EpRac = conj(EpRa);
+
draa_dt = -i*EpRac*rab+i*EpRa*rba-i*EdLac*rac+i*EdLa*rca-2*g*raa
-i*EdRac*rAB+i*EdRa*rBA-i*EpLac*rAC+i*EpLa*rCA;