computedT;

T(iny,inx) = T(iny,inx) + dT; 
setTBCs;