added logging, started harmonicThermal method
fixed absorption imaging

incorporated spatial sigma variation (due to intensity) in error calc.
need to put in the main absorption imaging, must reshape array