diff options
| -rw-r--r-- | fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 26 |
1 files changed, 20 insertions, 6 deletions
diff --git a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 index eec2884..2230cb0 100644 --- a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 +++ b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 @@ -24,17 +24,32 @@ program FourLevelPulseProp_v3_Double ! Four Level Notes: Valid to first order in dt ! ! Version 3: Make a step function in the coupling field. + + ! --------------------------------------------------------------------------- + ! SYSTEM DEFAULTS + ! --------------------------------------------------------------------------- + integer zpts + parameter (zpts=10) ! Number of slices in z direction + + ! default values for peak amplitudes + real*8 Om1peak_def,Om2peak_def,Omcpeak_def ! default values for field peak amplitudes + parameter (Om1peak_def=0.01,Om2peak_def=5.00,Omcpeak_def=0.0) + ! --------------------------------------------------------------------------- character*150 fname - integer nmat, npts, Nframe, tpts_max, NSkip, NWrite, zpts, zptsmax + integer nmat, npts, Nframe, tpts_max, NSkip, NWrite, zptsmax integer iSkipped; 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 (zptsmax=200) ! max Number of slices in z direction parameter (tpts_max=200) ! max Number of slices in t direction parameter (NWrite=200) !number of frames to actually write + + + real*8 Lcell,Om1peak,Om2peak,Omcpeak,pi,tmax,tmin,tp,tshift,t_end,t_start,t_elapsed + + integer m,n complex*16 a1,a2,a3,a4,a5,a6 complex*16 b1,b2,b3,b4,b5,b6,b7 @@ -54,7 +69,6 @@ program FourLevelPulseProp_v3_Double real*8 dt,dz,Ndensity real*8 W12,W21,W31,W32,W41,W42,W43,W34,ga12,ga13,ga14,ga23,ga24,ga34 real*8 Ga2,Ga4,Om_crit - 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 t,z(zptsmax) @@ -111,9 +125,9 @@ program FourLevelPulseProp_v3_Double ! User defined numbers ! Rabi frequency are unit less and scaled by Gamma_super i.e. Om1/Gamma_super -> Om1 - Om1peak=0.01 ! Field 1 peak scaled Rabi frequency for the pump at entrance of cell - Om2peak=5.00 ! Field 2 peak scaled Rabi frequency for the pump at entrance of cell - Omcpeak=0.0 ! Field 3 peak scaled Rabi frequency for the pump at entrance of cell + Om1peak=Om1peak_def ! Field 1 peak scaled Rabi frequency for the pump at entrance of cell + Om2peak=Om2peak_def ! Field 2 peak scaled Rabi frequency for the pump at entrance of cell + Omcpeak=Omcpeak_def ! Field 3 peak scaled Rabi frequency for the pump at entrance of cell ! Maximum detuning in MHz for dispersion lineshape plot delmax=0 Ga4=(W41+W42+W43) |