diff options
| -rw-r--r-- | faraday/task_builder.m | 6 | ||||
| -rwxr-xr-x | faraday/task_solver.m | 19 |
2 files changed, 19 insertions, 6 deletions
diff --git a/faraday/task_builder.m b/faraday/task_builder.m index 46dcd64..22c346f 100644 --- a/faraday/task_builder.m +++ b/faraday/task_builder.m @@ -34,9 +34,3 @@ for i=1:Nsteps save(filename_task, 'detuning_freq', 'Ep', 'psi_el', 'B_field', 'theta', 'phi'); endfor return; - - -%[psr_rad_smEp_pos_el]=psr_vs_detuning(detuning_freq, Ep, psi_el, B_field, theta, phi) ; -[psr_rad_smEp_pos_el]=faraday_vs_B(detuning_freq, Ep, psi_el, B_fields, theta, phi) ; - -plot(B_fields, psr_rad_smEp_pos_el) diff --git a/faraday/task_solver.m b/faraday/task_solver.m new file mode 100755 index 0000000..2b68bd3 --- /dev/null +++ b/faraday/task_solver.m @@ -0,0 +1,19 @@ +#!/usr/bin/octave -qf +% reads parameters file (argv1) and output results of the calculations into output file (argv(2) + +arg_list = argv (); + +if (nargin < 2) + error("need at least two parameters"); +end +params_file=arg_list{1}; +output_fname=arg_list{2}; + +load(params_file); + +[xi_linear, xi_left, xi_right]=propagation_problem(detuning_freq, Ep, psi_el, B_field, theta, phi) + +save(output_fname, 'detuning_freq', 'Ep', 'psi_el', 'B_field', 'theta', 'phi', ... + 'xi_linear', 'xi_left', 'xi_right' ... + ); + |