diff options
| -rw-r--r-- | psr/psr.m | 7 |
1 files changed, 4 insertions, 3 deletions
@@ -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) |