diff options
| -rw-r--r-- | fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 index 39fdfe7..7a9b8dd 100644 --- a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 +++ b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 @@ -62,7 +62,7 @@ program FourLevelPulseProp_v3_Double complex*16 l1,l2,l3,l4,l5 complex*16 Omold, Omold_vac - real*8 alpha1,alpha2,alpha1tilde,alpha2tilde,alphac,alphactilde,Gamma_super,c,delmax,del1_prop,del2_prop,delc_prop + real*8 alpha1,alpha2,alpha1tilde,alpha2tilde,alphac,alphactilde,Gamma_super,c,delmax,del1_prop,del2_prop,del3_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 @@ -157,8 +157,8 @@ program FourLevelPulseProp_v3_Double ! Now that the user has an idea of the dispersion, do the full propagation problem ! Detuning of center frequency of the coupling pulse in MHz - delc_prop=0 - delc_prop=2*pi*1e6*delc_prop/Gamma_super !Now dimensionless + del3_prop=0 + del3_prop=2*pi*1e6*del3_prop/Gamma_super !Now dimensionless ! Detuning of center frequency of the pump pulse in MHz del2_prop=0. @@ -270,7 +270,7 @@ program FourLevelPulseProp_v3_Double e4=0.25*ci*dt e5=-0.25*ci*dt - f1=1-(ga14-ci*(del2_prop-del1_prop-delc_prop))*dt + f1=1-(ga14-ci*(del2_prop-del1_prop-del3_prop))*dt f2=0.25*ci*dt f3=-0.25*ci*dt f4=0.25*ci*dt @@ -284,13 +284,13 @@ program FourLevelPulseProp_v3_Double gc6=0.25*ci*dt gc7=0.25*ci*dt - h1=1-(ga24+ci*delc_prop)*dt + h1=1-(ga24+ci*del3_prop)*dt h2=0.25*ci*dt h3=0.25*ci*dt h4=0.25*ci*dt h5=0.25*ci*dt - i1=1-(ga34-ci*(del2_prop-delc_prop))*dt + i1=1-(ga34-ci*(del2_prop-del3_prop))*dt i2=-0.25*ci*dt i3=0.25*ci*dt i4=0.25*ci*dt |