diff options
| -rw-r--r-- | MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds | 16 |
1 files changed, 16 insertions, 0 deletions
diff --git a/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds b/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds index 001b424..4c4262a 100644 --- a/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds +++ b/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds @@ -172,10 +172,26 @@ <![CDATA[ // Initial (at starting 'z' position) electromagnetic field does not depend on detuning // as well as time + + // constant input + /* EdL=EdLo; EdR=EdRo; EpL=EpLo; EpR=EpRo; + */ + + // amplitude noise modulated + double r1 = 2.0*(1.0*rand()/RAND_MAX) - 1.0; // random number -1 .. +1 + double r2 = 2.0*(1.0*rand()/RAND_MAX) - 1.0; // random number -1 .. +1 + double r3 = 2.0*(1.0*rand()/RAND_MAX) - 1.0; // random number -1 .. +1 + double common_ampl_noise_frac = 1+.0*r1; + double one_phot_noise_frac = .0*r2; + double two_phot_noise_frac = .3*r3; + EdL=EdLo * common_ampl_noise_frac * exp(i*one_phot_noise_frac); + EdR=EdRo * common_ampl_noise_frac * exp(i*one_phot_noise_frac); + EpL=EpLo * common_ampl_noise_frac * exp(i*one_phot_noise_frac) * exp(i*two_phot_noise_frac); + EpR=EpRo * common_ampl_noise_frac * exp(i*one_phot_noise_frac) * exp(i*two_phot_noise_frac); // Temporary pulse in all fields /* |