diff options
Diffstat (limited to 'psr')
| -rw-r--r-- | psr/psr.m | 14 |
1 files changed, 8 insertions, 6 deletions
@@ -32,7 +32,8 @@ N_detun_steps=100; detuning_p_min=-200.0; detuning_p_max=-detuning_p_min; detuning_p_max=1000; -detuning_freq=zeros(1,N_detun_steps+1); +detuning_freq=linspace(detuning_p_min,detuning_p_max,N_detun_steps); +%detuning_freq=zeros(1,N_detun_steps+1); kappa_p =zeros(1,N_detun_steps+1); kappa_m =zeros(1,N_detun_steps+1); detun_step=(detuning_p_max-detuning_p_min)/N_detun_steps; @@ -109,10 +110,11 @@ psi_el=-5/180*pi; fprintf (stderr, "tuning laser in forloop to set conditions vs detuning\n"); fflush (stderr); -for detuning_p_cntr=1:N_detun_steps+1; +for detuning_p_cntr=1:length(detuning_freq); wp0=w_pf1-w_sf2; %Fg=2 -> Fe=1 %wd=w_pf1-w_hpf_ground; - detuning_p=detuning_p_min+detun_step*(detuning_p_cntr-1); + %detuning_p=detuning_p_min+detun_step*(detuning_p_cntr-1); + detuning_p=detuning_freq(detuning_p_cntr); wp=wp0+detuning_p; light_positive_freq=[ wp]; % we calculate dc and negative frequiencies as well as amplitudes @@ -129,7 +131,7 @@ for detuning_p_cntr=1:N_detun_steps+1; %kappa_p(detuning_p_cntr)=susceptibility_steady_state_at_freq( atom_field_problem); - detuning_freq(detuning_p_cntr)=detuning_p; + %detuning_freq(detuning_p_cntr)=detuning_p; endfor save '/tmp/problem_definition.mat' problems_cell_array atom_properties detuning_freq ; @@ -140,9 +142,9 @@ fflush (stderr); %[xi_linear, xi_left, xi_right]=cellfun( @susceptibility_steady_state_at_freq, problems_cell_array); %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 ; +save '/tmp/xi_vs_detuning.mat' detuning_freq xi_linear xi_left xi_right E_field_pos_freq E_field_probe ; -output_psr_results_vs_detuning.m; +output_psr_results_vs_detuning; elapsed_time = etime (clock (), t0) |