summaryrefslogtreecommitdiff
path: root/fortran/navy_four_levels/FourLevelPulseProp_Double.f95
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2011-08-23 01:01:57 -0400
committerEugeniy Mikhailov <evgmik@gmail.com>2011-08-23 01:01:57 -0400
commita09a012fef408991f15be5bdfd97115b693185bc (patch)
treebe394ab79f2a28b894231056d0ced82de97c11ee /fortran/navy_four_levels/FourLevelPulseProp_Double.f95
parent08e22ea22ba897a1dc67b90a94148b60b0a02b80 (diff)
downloadNresonances-a09a012fef408991f15be5bdfd97115b693185bc.tar.gz
Nresonances-a09a012fef408991f15be5bdfd97115b693185bc.zip
No more then given number of step along time is outputed
Diffstat (limited to 'fortran/navy_four_levels/FourLevelPulseProp_Double.f95')
-rw-r--r--fortran/navy_four_levels/FourLevelPulseProp_Double.f9512
1 files changed, 9 insertions, 3 deletions
diff --git a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95
index 6ae00e7..da5d8fd 100644
--- a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95
+++ b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95
@@ -307,7 +307,9 @@ program FourLevelPulseProp_v3_Double
fname='Movie4level_EndPoints.dat' ! File name to save endpoints
open(10,file=fname)
+ iSkipped=0;
do n=1,Nframe
+ iSkipped = iSkipped+1
t=tmin+float(n-1)*dt
Om1(1)=Om1peak*exp(-(t-tshift)**2/(tp*tp))
Om2(1)=Om2peak
@@ -411,13 +413,17 @@ program FourLevelPulseProp_v3_Double
rho44_last(m)=rho44(m)
-
- write (9,120) t,z(m),cdabs(Om1(m)),cdabs(Om2(m)),cdabs(Omc(m)),cdabs(Om_vac(m))
+ if ( iSkipped.ge.NSkip ) then
+ write (9,120) t/Gamma_super,z(m)/Gamma_super*c,cdabs(Om1(m)),cdabs(Om2(m)),cdabs(Omc(m)),cdabs(Om_vac(m))
+ end if
end do
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(Om1(zpts)),cdabs(Om_vac(zpts)) !EndPoint File
+ if ( iSkipped.ge.NSkip ) then
+ write (10,139) t/Gamma_super,cdabs(Om1(zpts)),cdabs(Om_vac(zpts)) !EndPoint File
+ iSkipped = 0
+ end if
Omold=Om2(zpts)
Omold_vac=Om_vac(zpts)