diff options
Diffstat (limited to 'psr/psr_vs_detuning_combo.m')
| -rw-r--r-- | psr/psr_vs_detuning_combo.m | 65 |
1 files changed, 57 insertions, 8 deletions
diff --git a/psr/psr_vs_detuning_combo.m b/psr/psr_vs_detuning_combo.m index 18ad5cc..b2ea9d3 100644 --- a/psr/psr_vs_detuning_combo.m +++ b/psr/psr_vs_detuning_combo.m @@ -1,6 +1,7 @@ 1; data_dir='results/'; +output_dir='results/'; N_detun_steps=1000; detuning_p_min=-200.0; %detuning_p_max=-detuning_p_min; @@ -29,8 +30,15 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, psi_el, theta, phi); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el + , output_dir ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % 0.1 G magnetic field,, 30 degree ellipticity @@ -40,8 +48,15 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el + , output_dir ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % 0.0001 G magnetic field,, 30 degree ellipticity @@ -51,8 +66,15 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % 1.0 G magnetic field,, 30 degree ellipticity @@ -62,8 +84,15 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % -0.1 G magnetic field,, 30 degree ellipticity @@ -73,8 +102,15 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % -0.0001 G magnetic field,, 30 degree ellipticity @@ -84,8 +120,15 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % -1.0 G magnetic field,, 30 degree ellipticity @@ -95,6 +138,12 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) - +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); |