From 3e8234c27875209afc0a9a170f129b5bd23dd1f1 Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Thu, 12 May 2016 16:50:28 -0400
Subject: added naive noise components

---
 .../Nlevels_with_MOR.xmds                                | 16 ++++++++++++++++
 1 file changed, 16 insertions(+)

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
 			/*
-- 
cgit v1.2.3