diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2010-06-24 03:20:43 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2010-06-24 03:20:43 +0000 |
| commit | be16b3f4016341d4fac53eb9ccf352514a393c8b (patch) | |
| tree | 928115a024709f21bc303ff9d1b754a946b504cf | |
| parent | 4503717f89973c9d2eba9512ed3d16e6a5027303 (diff) | |
| download | multi_mode_eit-be16b3f4016341d4fac53eb9ccf352514a393c8b.tar.gz multi_mode_eit-be16b3f4016341d4fac53eb9ccf352514a393c8b.zip | |
tweak in ellipticity parameters
| -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) |