diff options
Diffstat (limited to 'fortran/navy_four_levels/FourLevelPulseProp_Double.f95')
| -rw-r--r-- | fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 18 |
1 files changed, 12 insertions, 6 deletions
diff --git a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 index 870e035..4e0cfc9 100644 --- a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 +++ b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 @@ -99,15 +99,18 @@ program FourLevelPulseProp_v3_Double ! ! default values for peak amplitudes real*8 Om1peak_def,Om2peak_def,Om3peak_def ! default values for field peak amplitudes - parameter( Om1peak_def=0.0025*(2*pi*1e6)/Gamma_super ) - parameter( Om2peak_def=6*(2*pi*1e6)/Gamma_super ) - parameter( Om3peak_def=0.0*(2*pi*1e6)/Gamma_super ) + parameter( Om1peak_def=2*0.0025*(2*pi*1e6)/Gamma_super ) + parameter( Om2peak_def=2*6*(2*pi*1e6)/Gamma_super ) + parameter( Om3peak_def=2*0.0*(2*pi*1e6)/Gamma_super ) ! -------------------------------------------------------------------------- ! Atomic numbers (based on Rubidium 87) ! Atomic numbers (based on Rubidium 87) lambda=794.7e-9 ! Wavelength in [m] Ndensity=1e10*(1e6) ! Density in m^-1 ! atomic decay parameters in units of Gamma_super + ! by Eugeniy 2011/08/26 + ! Note extra factor of 2 compared to xmds version, + ! we need to double check this. 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 @@ -132,11 +135,14 @@ program FourLevelPulseProp_v3_Double ! Atomic parameters write (*,*)'New version from moved folder' - alpha1=3*Ndensity*lambda*lambda/(2*pi) + ! by Eugeniy 2011/08/26 + ! Note 4*pi in the denominator vs original 2*pi and 8*pi of xmds version. + ! We need to double check this. + alpha1=3*Ndensity*lambda*lambda/(4*pi) alpha1tilde=alpha1*c/Gamma_super - alpha2=3*Ndensity*lambda*lambda/(2*pi) + alpha2=3*Ndensity*lambda*lambda/(4*pi) alpha2tilde=alpha2*c/Gamma_super - alphac=3*Ndensity*lambda*lambda/(2*pi) + alphac=3*Ndensity*lambda*lambda/(4*pi) alphactilde=alphac*c/Gamma_super |