diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2010-06-16 19:35:56 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2010-06-16 19:35:56 +0000 |
| commit | 45f970c010cd2981f86eb4b7ce3cd25f9947e96f (patch) | |
| tree | 43ad327822a75be66e7587b4f024961600f09aa2 /psr | |
| parent | 9709621e8c71e9a86cd6e3e47edc3d8b8e6becd9 (diff) | |
| download | multi_mode_eit-45f970c010cd2981f86eb4b7ce3cd25f9947e96f.tar.gz multi_mode_eit-45f970c010cd2981f86eb4b7ce3cd25f9947e96f.zip | |
small change in initial conditions
Diffstat (limited to 'psr')
| -rw-r--r-- | psr/field_description.m | 2 | ||||
| -rw-r--r-- | psr/psr.m | 9 |
2 files changed, 4 insertions, 7 deletions
diff --git a/psr/field_description.m b/psr/field_description.m index 738ebe3..e9fa227 100644 --- a/psr/field_description.m +++ b/psr/field_description.m @@ -1,7 +1,7 @@ 1; %EM field definition -Ed=10.12; %drive +Ed=0.002; %drive Edc=conj(Ed); Ep=0.8*Ed; %probe Epc=conj(Ep); @@ -29,7 +29,7 @@ field_description; detuning_p=0; N_detun_steps=100; %detuning_p_min=-B_field*gmg*4; % span +/-4 Zeeman splitting -detuning_p_min=-10040.0; +detuning_p_min=-40.0; detuning_p_max=-detuning_p_min; detuning_freq=zeros(1,N_detun_steps+1); kappa_p =zeros(1,N_detun_steps+1); @@ -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=-3/180*pi; +psi_el=-2/180*pi; % we define light as linearly polarized % where phi is angle between light polarization and axis x @@ -133,11 +133,8 @@ save '/tmp/problem_definition.mat' problems_cell_array atom_properties detuning fprintf (stderr, "now really hard calculations begin\n"); fflush (stderr); % once we define all problems the main job is done here -%kappa_p=cellfun( @susceptibility_steady_state_at_freq, problems_cell_array); -%kappa_p=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array); [xi_linear, xi_left, xi_right]=parcellfun(2, @susceptibility_steady_state_at_freq, problems_cell_array); -%relative_transmission_vs_detuning=parcellfun(2, @total_relative_transmission, problems_cell_array); -%relative_transmission_vs_detuning=cellfun(@total_relative_transmission, 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; |