diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2011-07-12 22:22:19 -0400 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2011-07-12 22:31:09 -0400 |
| commit | 6314afd82280074bbd71c74d17d872e1a802484d (patch) | |
| tree | 881d0f13960fe362db250eec79a69dd40b190767 /fortran | |
| parent | 8aabf9626b71c93160af9091dfc4a895ee1ae5c6 (diff) | |
| download | Nresonances-6314afd82280074bbd71c74d17d872e1a802484d.tar.gz Nresonances-6314afd82280074bbd71c74d17d872e1a802484d.zip | |
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
Diffstat (limited to 'fortran')
| -rw-r--r-- | fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 14 |
1 files changed, 7 insertions, 7 deletions
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) |