Glass production has a great industrial importance and is associated with many technological challenges. Control related problems concern especially the last part of the process, so called glass conditioning. Molten glass is gradually cooled down in a long ceramic channels called forehearths during glass conditioning. The glass temperature in each zone of the forehearth should be precisely adjusted according to the assumed profile. Due to cross-couplings and unmeasured disturbances, traditional control systems based on PID controllers, often do not ensure sufficient control quality. This problem is the main motivation for the research presented in the paper. A Model Predictive Control algorithm is proposed for the analysed process. It is assumed the dynamic model for each zone of the forehearth is identified on-line with the Modulating Functions Method. These continuous-time linear models are subsequently used for two purposes: for the predictive controller tuning, measurable disturbances compensation and for a static set point optimisation. Proposed approach was tested using Partial Differential Equation model to simulate two adjacent zones of the forehearth. The experimental results proved that it can be successfully applied for the aforementioned model.