Fixed issue with boundary conditions -- I needed the RHS to be fourier coefficients for each mode at the temperature boundaries.