From e17a65ede463a43cb59f59ef92a6e6d2a4bca245 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Sat, 10 Jul 2010 13:02:20 +0000 Subject: added special file for plotting --- psr/ouput_psr_vs_detuning_combo.m | 36 ++++++++++++++++++++++++++++++++++++ 1 file changed, 36 insertions(+) create mode 100644 psr/ouput_psr_vs_detuning_combo.m diff --git a/psr/ouput_psr_vs_detuning_combo.m b/psr/ouput_psr_vs_detuning_combo.m new file mode 100644 index 0000000..154effc --- /dev/null +++ b/psr/ouput_psr_vs_detuning_combo.m @@ -0,0 +1,36 @@ +function 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) + + %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + figure(6); + tn_cl='blue'; + sm_cl='green'; + lg_cl='magenta'; + gt_cl='red'; + + plot ( ... + detuning_freq, psr_rad_tnEp_pos_el, '-;tn,pos el;', "linewidth",3, 'color',tn_cl + , detuning_freq, psr_rad_tnEp_neg_el, '-;tn,neg el;', "linewidth",1, 'color',tn_cl + , detuning_freq, psr_rad_smEp_pos_el, '-;sm,pos el;', "linewidth",3, 'color',sm_cl + , detuning_freq, psr_rad_smEp_neg_el, '-;sm,neg el;', "linewidth",1, 'color',sm_cl + , detuning_freq, psr_rad_lgEp_pos_el, '-;lg,pos el;', "linewidth",3, 'color',lg_cl + , detuning_freq, psr_rad_lgEp_neg_el, '-;lg,neg el;', "linewidth",1, 'color',lg_cl + , detuning_freq, psr_rad_grEp_pos_el, '-;gr,pos el;', "linewidth",3, 'color',gt_cl + , detuning_freq, psr_rad_grEp_neg_el, '-;gr,neg el;', "linewidth",1, 'color',gt_cl + ); + + str_title=sprintf("PSR. B field=%.5f Gauss, ellipticity=%.2f rad, theta=%.2f, phi=%.2f", B_field, psi_el, theta, phi); + title(str_title); + xlabel('detuning, MHz'); + ylabel('PSR, radians'); + fname=data_file_name(output_dir, 'PSR.','pdf', B_field, theta,phi,psi_el); + print(fname); + + ret=true; +endfunction + -- cgit v1.2.3