diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2011-07-16 00:16:31 -0400 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2011-07-16 00:16:31 -0400 |
| commit | c9cd26c2f61ebf64daf0b8ec067b2fe689260c54 (patch) | |
| tree | b9c56b8d71d0cc1f131ea8e5875e16bc2a7f3dcf /fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | |
| parent | fba5cedc1dfe1d1c3989fd39413ae171e9fce6b0 (diff) | |
| download | Nresonances-c9cd26c2f61ebf64daf0b8ec067b2fe689260c54.tar.gz Nresonances-c9cd26c2f61ebf64daf0b8ec067b2fe689260c54.zip | |
Add tmin and move pulse center to t=0
Diffstat (limited to 'fortran/navy_four_levels/FourLevelPulseProp_Double.f95')
| -rw-r--r-- | fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 11 |
1 files changed, 6 insertions, 5 deletions
diff --git a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 index 2873180..096f293 100644 --- a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 +++ b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 @@ -52,7 +52,7 @@ program FourLevelPulseProp_v3_Double real*8 dt,dz,eta real*8 W12,W21,W31,W32,W41,W42,W43,W34,ga12,ga13,ga14,ga23,ga24,ga34 real*8 Ga2,Ga4,Om_crit - real*8 Lcell,Om1peak,Om2peak,Omcpeak,pi,tmax,tp,tshift,t_end,t_start,t_elapsed + real*8 Lcell,Om1peak,Om2peak,Omcpeak,pi,tmax,tmin,tp,tshift,t_end,t_start,t_elapsed real*8 tpeak,tpeak_vac real*8 epsil,hbar,lambda real*8 del(npts) @@ -153,17 +153,18 @@ program FourLevelPulseProp_v3_Double ! ! Set up initial pulse. ! - tshift=2*tp + tshift=0 write (*,*) 'peak center at the cell begining i.e. tshift = ', tshift tmax=Lcell*20 !Length of time to pass cell (no c because we're dimensionless) - tmax=130.0 + tmax=60.0 + tmin = -tmax write(*,*) 'tmax = ', tmax write(*,*) 'Lcell = ', Lcell dz=Lcell/(zpts-1) !(no c because we're dimensionless) ! It is crucial that dz = dt in unitless coordinates, there is built in ! cancellation of some term on grid because of it. See Simon's note. dt = dz - Nframe=tmax/dz+1 + Nframe=(tmax-tmin)/dz+1 if (Nframe.ge.Nframemax) write (*,*)'Error!!!!Nframe>Nframemax' ! ! Initialize matrices @@ -310,7 +311,7 @@ program FourLevelPulseProp_v3_Double open(10,file=fname) do n=1,Nframe - t=float(n-1)*dt + t=tmin+float(n-1)*dt Om1(1)=Om1peak*exp(-(t-tshift)**2/(tp*tp)) Om2(1)=Om2peak Omc(1)=Omcpeak |