diff options
Diffstat (limited to 'fortran/navy_four_levels/FourLevelPulseProp_Double.f95')
| -rw-r--r-- | fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 14 |
1 files changed, 3 insertions, 11 deletions
diff --git a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 index 096f293..600a6ab 100644 --- a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 +++ b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 @@ -26,14 +26,14 @@ program FourLevelPulseProp_v3_Double ! Version 3: Make a step function in the coupling field. character*150 fname - integer nmat,npts,Nfrac,Nframe,Nframemax,NSkip,NWrite,tpts,zpts + integer nmat,npts,Nframe,Nframemax,NSkip,NWrite,zpts parameter (nmat=3,npts=100) !matrix size, number of detuning points in dispersion curve !tpts is the number of temporal points in the cell !zpts is the number of z output slices - 1 parameter (zpts=100) ! Number of slices in z direction parameter (Nframemax=2000000) parameter (NWrite=1) !number of frames to actually write - integer i,j,k,m,n + integer m,n complex*16 a1,a2,a3,a4,a5,a6 complex*16 b1,b2,b3,b4,b5,b6,b7 complex*16 c1,c2,c3,c4,c5,c6 @@ -55,10 +55,7 @@ program FourLevelPulseProp_v3_Double 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) real*8 t,z(zpts) - complex*16 yplot(nmat,npts) - complex*16 Imat(nmat) complex*16 Om1(zpts),Om2(zpts),Omc(zpts),Om_vac(zpts) @@ -70,14 +67,10 @@ program FourLevelPulseProp_v3_Double complex*16 rho41_last(zpts),rho42_last(zpts),rho43_last(zpts),rho44_last(zpts) !No Om_last because we never need the previous spatial point - complex*16 L(nmat,nmat),Linv(nmat,nmat),Ltemp(nmat,nmat) common/para/ga12,W21 - real*8 d !used by NR Routines - integer indx(nmat) - ! Fundamental numbers ci=cmplx(0.,1.) pi=acos(-1.0) @@ -298,7 +291,6 @@ program FourLevelPulseProp_v3_Double NSkip=int(NFrame/NWrite) fname='MovieParameters4level.txt' ! File name to save parameters -3 format(a150) open(9,file=fname) write (9,133)Nframe,zpts,beta,NSkip,dt 133 format(1x,i10,',',i5,',',f12.2,',',i5,',',f12.2) @@ -334,7 +326,7 @@ program FourLevelPulseProp_v3_Double rho42_last(m)=rho42(m) rho43_last(m)=rho43(m) rho44_last(m)=rho44(m) -345 end do + end do do m=zpts,2,-1 z(m)=float(m-1)*dz |