1) Added CT but problems with the zones in touch with the boundary.
2) Added preliminary conduction code.
3) Issues with robustness: Need to set GAMMA_MAX to 5 and X2_START to 0.001.