From b1ba94456e3c6a7a042f77925d9b22a5f5dfaaaa Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Tue, 23 Aug 2011 11:14:29 -0400 Subject: Added output at z=0 and at the very last non skipped time point --- fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 18 +++++++++++++----- 1 file changed, 13 insertions(+), 5 deletions(-) (limited to 'fortran/navy_four_levels/FourLevelPulseProp_Double.f95') diff --git a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 index 3bb5eea..eec2884 100644 --- a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 +++ b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 @@ -308,7 +308,7 @@ program FourLevelPulseProp_v3_Double fname='Movie4level_EndPoints.dat' ! File name to save endpoints open(10,file=fname) - iSkipped=0; + iSkipped=NSkip; ! since we want initial time points to be at output as well do n=1,Nframe iSkipped = iSkipped+1 t=tmin+float(n-1)*dt @@ -414,15 +414,23 @@ program FourLevelPulseProp_v3_Double rho44_last(m)=rho44(m) - if ( iSkipped.ge.NSkip ) then + if ( iSkipped.ge.NSkip ) then ! Outputting every Nskip frame 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 + ! the very first z point is skipped in above cycle we compencete for it here + if ( iSkipped.ge.NSkip ) then ! Outputting every Nskip frame + write (9,120) t/Gamma_super,z(1)/Gamma_super*c,cdabs(Om1(1)),cdabs(Om2(1)),cdabs(Omc(1)),cdabs(Om_vac(1)) + end if + if (cdabs(Om2(zpts)).gt.cdabs(Omold)) tpeak=t if (cdabs(Om_vac(zpts)).gt.cdabs(Omold_vac)) tpeak_vac=t - if ( iSkipped.ge.NSkip ) then - write (10,139) t/Gamma_super,cdabs(Om1(zpts)),cdabs(Om_vac(zpts)) !EndPoint File + + if ( iSkipped.ge.NSkip ) then ! Outputting every Nskip frame + ! Outputting Rabi frequency at the beginning and the end of + ! the cell for given time to the EndPoint File + write (10,139) t/Gamma_super, cdabs(Om1(1)), cdabs(Om_vac(1)), cdabs(Om1(zpts)), cdabs(Om_vac(zpts)) iSkipped = 0 end if Omold=Om2(zpts) @@ -432,7 +440,7 @@ program FourLevelPulseProp_v3_Double close(9) close(10) -139 format(1x,E15.9,',',E15.9,',',E15.9) +139 format(1x,E15.9,',',E15.9,',',E15.9,',',E15.9,',',E15.9) 120 format(E15.9,',',E15.9,',',E15.9,',',E15.9,',',E15.9,',',E15.9,',',E15.9,',',E15.9,',',E15.9) t_end=secnds(0.E0) t_elapsed=t_end-t_start -- cgit v1.2.3