Code runs and produces sensible results.  Still no E(b-V)