From be16b3f4016341d4fac53eb9ccf352514a393c8b Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 24 Jun 2010 03:20:43 +0000 Subject: tweak in ellipticity parameters --- psr/psr.m | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/psr/psr.m b/psr/psr.m index 66191bb..3aefad8 100644 --- a/psr/psr.m +++ b/psr/psr.m @@ -91,7 +91,7 @@ phi=pi*2/8; % theta is angle between lab z axis (light propagation direction) and magnetic field axis (z') theta=0; % psi_el is the ellipticity parameter (phase difference between left and right polarization) -psi_el=-2/180*pi; +psi_el=-5/180*pi; % we define light as linearly polarized % where phi is angle between light polarization and axis x @@ -100,6 +100,7 @@ psi_el=-2/180*pi; [E_field_lab_pos_freq.x, E_field_lab_pos_freq.y] = rotXpolarization(phi, E_field_lab_pos_freq.linear); % we add required ellipticity E_field_lab_pos_freq.x*=exp(I*psi_el); + E_field_lab_pos_freq.y*=exp(-I*psi_el); E_field_lab_pos_freq.z=E_field_zero; E_field_pos_freq=xyz_lin2atomic_axis_polarization(theta, E_field_lab_pos_freq); @@ -137,10 +138,10 @@ fflush (stderr); [xi_linear, xi_left, xi_right]=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array); %[xi_linear, xi_left, xi_right]=cellfun( @susceptibility_steady_state_at_freq, problems_cell_array); -save '/tmp/xi_vs_detuning.mat' detuning_freq xi_linear xi_left xi_right ; %save '/tmp/relative_transmission_vs_detuning.mat' detuning_freq relative_transmission_vs_detuning; +save '/tmp/xi_vs_detuning.mat' detuning_freq xi_linear xi_left xi_right E_field_pos_freq ; -output_xi_results; +output_psr_results; elapsed_time = etime (clock (), t0) -- cgit v1.2.3