diff options
| -rw-r--r-- | fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 12 |
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) |