generalized and fixed ground levels decay

Now each of 4 ground levels decay into other 3.

1 files changed, 10 insertions, 4 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 8c11b5a..001b424 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
@@ -328,14 +328,20 @@
 					complex EpLac = conj(EpLa);
 					complex EpRac = conj(EpRa);
 
+					// ground level repopulation rate, i.e. arrival from other levels decay
+					// note that a level repopulate itself, thus there is 4 in 
+					// a term like drbb_dt = ... + Repop - 4*gbc*rbb;
+					// which otherwise decays only in 3 other ground levels
+					complex Repop = Gamma_super*raa + gbc*(rbb+rcc+rBB+rCC);
+
 					draa_dt = -i*EpRac*rab+i*EpRa*rba-i*EdLac*rac+i*EdLa*rca-4*Gamma_super*raa
 					          -i*EpLac*rAB+i*EpLa*rBA-i*EdRac*rAC+i*EdRa*rCA;
 
-					drbb_dt =  i*EpRac*rab-i*EpRa*rba+Gamma_super*raa-gbc*rbb+gbc*rcc;
-					drBB_dt =  i*EpLac*rAB-i*EpLa*rBA+Gamma_super*raa-gbc*rBB+gbc*rCC;
+					drbb_dt =  i*EpRac*rab-i*EpRa*rba + Repop - 4*gbc*rbb;
+					drBB_dt =  i*EpLac*rAB-i*EpLa*rBA + Repop - 4*gbc*rBB;
 
-					drcc_dt =  i*EdLac*rac-i*EdLa*rca+Gamma_super*raa-gbc*rcc+gbc*rbb;
-					drCC_dt =  i*EdRac*rAC-i*EdRa*rCA+Gamma_super*raa-gbc*rCC+gbc*rBB;
+					drcc_dt =  i*EdLac*rac-i*EdLa*rca + Repop - 4*gbc*rcc;
+					drCC_dt =  i*EdRac*rAC-i*EdRa*rCA + Repop - 4*gbc*rCC;
 
 					drab_dt = -Gab*rab+i*EpRa*(rbb-raa)+i*EdLa*rcb;
 					drca_dt = -Gca*rca+i*EdLac*(raa-rcc)-i*EpRac*rcb;