Small additions and fix for 0 flux and power density points

fit contains simple trajectory corrections.

fix for points with 0 flux or power density

simpler polarization calculations (to be highlighted more in v2.1 release)