diff options
Diffstat (limited to 'fortran/navy_four_levels')
| -rw-r--r-- | fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 index 746b1e0..f3b695d 100644 --- a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 +++ b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 @@ -323,7 +323,7 @@ ! write (*,*)'Enter file name to save endpoints' open(10,file=fname) - do 60 n=1,Nframe + do n=1,Nframe ! write (*,*)'n = ',n t=float(n-1)*dt Om1(1)=Om1peak*exp(-(t-tshift)**2/(tp*tp)) @@ -331,7 +331,7 @@ Omc(1)=Omcpeak Om_vac(1)=Om1(1) - do 345 m=1,zpts + do m=1,zpts rho11_last(m)=rho11(m) rho12_last(m)=rho12(m) rho13_last(m)=rho13(m) @@ -348,9 +348,9 @@ rho42_last(m)=rho42(m) rho43_last(m)=rho43(m) rho44_last(m)=rho44(m) -345 continue +345 end do - do 50 m=zpts,2,-1 + do m=zpts,2,-1 z(m)=float(m)*dz Om1(m)=a1*Om1(m-1)+a2*rho31_last(m)+a3*rho31_last(m-1) @@ -430,14 +430,14 @@ if (mod(n,Nskip).eq.0) write (9,120) t,z(m),cdabs(Om1(m)),cdabs(Om2(m)),cdabs(Omc(m)) -50 continue + 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 Omold=Om2(zpts) Omold_vac=Om_vac(zpts) -60 continue + end do close(9) close(10) 139 format(1x,f12.6,',',F12.6,',',F12.6) |