summaryrefslogtreecommitdiff
path: root/fortran
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2011-08-26 11:43:30 -0400
committerEugeniy Mikhailov <evgmik@gmail.com>2011-08-26 11:43:30 -0400
commit28582268a8cda10d477c92a7c15f5f144898e4e4 (patch)
tree4323e4635421e360dadbd298edf9c71a927ed003 /fortran
parent0c6c1aa2dab4cc09ba6636f6ff4fd360f8891fde (diff)
downloadNresonances-28582268a8cda10d477c92a7c15f5f144898e4e4.tar.gz
Nresonances-28582268a8cda10d477c92a7c15f5f144898e4e4.zip
delc --> del3. The coupling field is just a field 3 now
Diffstat (limited to 'fortran')
-rw-r--r--fortran/navy_four_levels/FourLevelPulseProp_Double.f9512
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