diff options
Diffstat (limited to 'fortran/navy_four_levels/FourLevelPulseProp_v2_Double.f95')
| -rw-r--r-- | fortran/navy_four_levels/FourLevelPulseProp_v2_Double.f95 | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/fortran/navy_four_levels/FourLevelPulseProp_v2_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_v2_Double.f95 index 3f6ce1c..045f4a3 100644 --- a/fortran/navy_four_levels/FourLevelPulseProp_v2_Double.f95 +++ b/fortran/navy_four_levels/FourLevelPulseProp_v2_Double.f95 @@ -92,7 +92,7 @@ ! beta=2*pi*3e6 !in Hz W41=0 - W42=0. + W42=1 W43=0 W32=1. W31=1. @@ -142,7 +142,7 @@ Om1peak=1 ! write (*,*)'Enter peak scaled Rabi frequency for the probe at entrance of cell' ! read (*,*)Om2peak - Om2peak=1e-2 + Om2peak=.01 ! write (*,*)'Enter maximum detuning in MHz for dispersion lineshape plot' ! read (*,*) delmax delmax=0 @@ -155,7 +155,7 @@ ! write (*,*)'Enter peak scaled Rabi frequency for the coupling field at entrance of cell' ! read (*,*)Omcpeak - Omcpeak=0. + Omcpeak=0.1 ! ! First plot the dispersion lineshape @@ -183,7 +183,7 @@ ! write (*,*)'Enter detuning of center frequency of the pump pulse in MHz' ! read (*,*)del2_prop !del_prop is the detuning used for the propagation - del2_prop=0; + del2_prop=0. del2_prop=2*pi*1e6*del2_prop/beta !Now dimensionless ! write (*,*)'Enter detuning of center frequency of the probe pulse in MHz' ! read (*,*)del1_prop !del_prop is the detuning used for the propagation @@ -197,7 +197,7 @@ tp=beta*tp !Now dimensionless ! write (*,*)'Enter length of cell in m' ! read (*,*)Lcell - Lcell=100; + Lcell=1 Lcell=beta*Lcell/c !Now dimensionless t_start=secnds(0.E0) write (*,*)'t_start = ',t_start @@ -210,7 +210,7 @@ dt=tmax/tpts dz=dt !(no c because we're dimensionless) ! write (*,*)'tp = ',tp - Nframe=zpts+int(6*tp/dt)+1 !Change the number 4 to anything you want to see longer pulse evolution + Nframe=zpts+int(4*tp/dt)+1 !Change the number 4 to anything you want to see longer pulse evolution if (Nframe.ge.Nframemax) write (*,*)'Error!!!!Nframe>Nframemax' ! write (*,*)'Nframe,tpts = ',Nframe,tpts ! @@ -328,7 +328,7 @@ j7=-0.25*ci*dt j8=W41*dt - k1=1-(W32+W31)*dt + k1=1-(W32+W31+W34)*dt k2=W43*dt k3=-0.25*ci*dt k4=-0.25*ci*dt |