Numerical experiments are performed on the chaotic movement of a simplified cli mate model. Considering various factors related to the predictability, we numeri cally compare the harmonized multi-time difference scheme (HMTDS) with the varia tional assimilation. It is shown that the HMTDS corrects remarkably the model bi as, especially the differential equation bias, which is just the drawback of the variational assimilation.