for->parfor
update limit on polarization parameters
add comments
output a struct instead of separate parameters
adjust to new output of irf_ebsp
correct calculation of PF
Merge branch 'master' of github.com:irfu/irfu-matlab