Fixed bugs in conduction part of dU_dt. Linear modes matches with
BalbusSaur. Both for stable and unstable modes. For some reason, it
converges only at 1st order. Needs more debugging. Ugh.