Ignoring PNG.
Added ability to do pol-intensity-weighted average.
Added flux % estimates for Attard 2009, sent to me by Nicholas Chapman.
Now differentiating between close and far SFRs.  Plotting close as Xs, not plotting far-away ones.
SHARP bolo vector scaling is now the same as SCUBA and Hertz, since I now I have the data for IRAS4A and IRAS4B.  I commented out the part about I = 999 for sources without intensity data.
Calculated KS statistics for all four distributions, ruling out the tightly aligned scenario.
For the referee: taking the difference between the lg-sm consistency angles after weighting by total intensity and then by polarized intensity.