this works.  needs a lot of cleaning up.  Also needs to make sure the variance for broad band is calculated correctly.
made changes that are compatible with photom_vs_spectral.py