diff options
Diffstat (limited to 'fortran/navy_four_levels')
| -rw-r--r-- | fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 22 |
1 files changed, 15 insertions, 7 deletions
diff --git a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 index e00e49a..83612a1 100644 --- a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 +++ b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 @@ -68,6 +68,8 @@ program FourLevelPulseProp_v3_Double real*8 alpha1,alpha2,alpha1tilde,alpha2tilde,alphac,alphactilde,Gamma_super,c,delmax,del1_prop,del2_prop,delc_prop real*8 dt,dz,Ndensity real*8 W12,W21,W31,W32,W41,W42,W43,W34,ga12,ga13,ga14,ga23,ga24,ga34 + real*8 R41, R42, R31, R32 ! splitting ratios of decays + real*8 G3, G4 ! total decay rate for given level i (Gi) real*8 Ga2,Ga4,Om_crit real*8 tpeak,tpeak_vac real*8 epsil,hbar,lambda @@ -95,14 +97,20 @@ program FourLevelPulseProp_v3_Double lambda=794.7e-9 ! Wavelength in [m] Ndensity=1e10*(1e6) ! Density in m^-1 ! atomic decay parameters in units of Gamma_super - W41=0 - W42=1 - W43=0 - W32=1. - W31=1. - W21=.001 + G3 = 2*2.7 *(2*pi*1e6) / Gamma_super ! decay rate for level 3 + G4 = 2*3.0 *(2*pi*1e6) / Gamma_super ! decay rate for level 4 + R41 = 0.5 ! splitting ratios for level i->j for Rij + R42 = 0.5 + R31 = 0.5 + R32 = 0.5 + W41 = R41*G4 + W42 = R42*G4 + W43 = 0 ! transition 4 <-> 3 is forbiden + W34 = 0 + W32 = R32*G3 + W31 = R31*G3 + W21=.01 * (2*M_PI*1e6) / Gamma_super W12=W21 - W34=0 ga12=0.5*(W21+W12) ga13=0.5*(W31+W12+W32) ga14=0.5*(W41+W42+W12) |