diff options
| -rw-r--r-- | fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 48 |
1 files changed, 24 insertions, 24 deletions
diff --git a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 index ab340d3..b4229a0 100644 --- a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 +++ b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 @@ -184,30 +184,30 @@ ! pause ! do 110 n=1,Nframe - do 100 m=1,zpts - Om1(m)=cmplx(0.,0.) - Om2(m)=cmplx(0.,0.) - Omc(m)=cmplx(0.,0.) - - Om_vac(m)=cmplx(0.,0.) - rho11(m)=cmplx(1.,0.) !Change this to change the initial condition - rho12(m)=cmplx(0.,0.) - rho13(m)=cmplx(0.,0.) - rho14(m)=cmplx(0.,0.) - rho21(m)=cmplx(0.,0.) - rho22(m)=cmplx(0.,0.) - rho23(m)=cmplx(0.,0.) - rho24(m)=cmplx(0.,0.) - rho31(m)=cmplx(0.,0.) - rho32(m)=cmplx(0.,0.) - rho33(m)=cmplx(0.,0.) - rho34(m)=cmplx(0.,0.) - rho41(m)=cmplx(0.,0.) - rho42(m)=cmplx(0.,0.) - rho43(m)=cmplx(0.,0.) - rho44(m)=cmplx(0.,0.) - -100 continue + do m=1,zpts + Om1(m)=cmplx(0.,0.) + Om2(m)=cmplx(0.,0.) + Omc(m)=cmplx(0.,0.) + + Om_vac(m)=cmplx(0.,0.) + rho11(m)=cmplx(1.,0.) !Change this to change the initial condition + rho12(m)=cmplx(0.,0.) + rho13(m)=cmplx(0.,0.) + rho14(m)=cmplx(0.,0.) + rho21(m)=cmplx(0.,0.) + rho22(m)=cmplx(0.,0.) + rho23(m)=cmplx(0.,0.) + rho24(m)=cmplx(0.,0.) + rho31(m)=cmplx(0.,0.) + rho32(m)=cmplx(0.,0.) + rho33(m)=cmplx(0.,0.) + rho34(m)=cmplx(0.,0.) + rho41(m)=cmplx(0.,0.) + rho42(m)=cmplx(0.,0.) + rho43(m)=cmplx(0.,0.) + rho44(m)=cmplx(0.,0.) + + end do !110 continue |