From 6314afd82280074bbd71c74d17d872e1a802484d Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Tue, 12 Jul 2011 22:22:19 -0400 Subject: Calculate canonic EIT case * field 1 is weak and has Gaussian temporal profile * field 2 is strong * field 3 is off * vacuum propagated field start similar to field 1 * cell length enlarged to see pulse delay --- fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) (limited to 'fortran') diff --git a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 index 406b1b1..eaa22ef 100644 --- a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 +++ b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 @@ -139,10 +139,10 @@ ! ! write (*,*)'Enter peak scaled Rabi frequency for the pump at entrance of cell' ! read (*,*)Om1peak - Om1peak=1 + Om1peak=0.01 ! write (*,*)'Enter peak scaled Rabi frequency for the probe at entrance of cell' ! read (*,*)Om2peak - Om2peak=.01 + Om2peak=1 ! write (*,*)'Enter maximum detuning in MHz for dispersion lineshape plot' ! read (*,*) delmax delmax=0 @@ -197,7 +197,7 @@ tp=beta*tp !Now dimensionless ! write (*,*)'Enter length of cell in m' ! read (*,*)Lcell - Lcell=1 + Lcell=100 Lcell=beta*Lcell/c !Now dimensionless t_start=secnds(0.E0) write (*,*)'t_start = ',t_start @@ -373,10 +373,10 @@ do 60 n=1,Nframe ! write (*,*)'n = ',n t=float(n-1)*dt - Om1(1)=Om1peak - Om2(1)=Om2peak*exp(-(t-tshift)**2/(tp*tp)) + Om1(1)=Om1peak*exp(-(t-tshift)**2/(tp*tp)) + Om2(1)=Om2peak Omc(1)=Omcpeak - Om_vac(1)=Om2(1) + Om_vac(1)=Om1(1) if (int(n/10).eq.0) write(fname,130)'Movie',n if (int(n/10).ge.1.and.int(n/100).eq.0) write (fname,131)'Movie',n if (int(n/10).ge.1.and.int(n/100).gt.0) write (fname,132)'Movie',n @@ -489,7 +489,7 @@ 50 continue if (cdabs(Om2(zpts)).gt.cdabs(Omold)) tpeak=t if (cdabs(Om_vac(zpts)).gt.cdabs(Omold_vac)) tpeak_vac=t - write (10,139) t,cdabs(Om2(zpts)),cdabs(Om_vac(zpts)) !EndPoint File + write (10,139) t,cdabs(Om1(zpts)),cdabs(Om_vac(zpts)) !EndPoint File Omold=Om2(zpts) Omold_vac=Om_vac(zpts) -- cgit v1.2.3