summaryrefslogtreecommitdiff
path: root/fortran
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2011-08-23 13:47:27 -0400
committerEugeniy Mikhailov <evgmik@gmail.com>2011-08-23 13:47:27 -0400
commitfd1062474d4fc61c943b9640ab51ac62d0819abb (patch)
treee649a36323690e1fd1dd359a5f2383e90ef719f7 /fortran
parent9214d37cb22ac28fc64c152dd6b222481dbcb351 (diff)
downloadNresonances-fd1062474d4fc61c943b9640ab51ac62d0819abb.tar.gz
Nresonances-fd1062474d4fc61c943b9640ab51ac62d0819abb.zip
Moved forward the definition for the peak Rabi frequencies
Diffstat (limited to 'fortran')
-rw-r--r--fortran/navy_four_levels/FourLevelPulseProp_Double.f9526
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)