diff options
| -rw-r--r-- | faraday/tasks_pp.m | 19 |
1 files changed, 19 insertions, 0 deletions
diff --git a/faraday/tasks_pp.m b/faraday/tasks_pp.m new file mode 100644 index 0000000..bc81b07 --- /dev/null +++ b/faraday/tasks_pp.m @@ -0,0 +1,19 @@ + +fnames=glob('results/*.mat'); + +Nsteps=length(fnames) +B_field=zeros(1,Nsteps); +xi_linear=zeros(1,Nsteps); +xi_left=zeros(1,Nsteps); +xi_right=zeros(1,Nsteps); + +% read the information from resulting files +for i=1:Nsteps + d=load(fnames{i}); + B_field(i)=d.B_field; + xi_linear(i)=d.xi_linear; + xi_left(i)=d.xi_left; + xi_right(i)=d.xi_right; +endfor + +plot(B_field, real(xi_right+xi_left)); |