In light of irreversibility of atmospheric motion, the memory function obtained by utilizing the previous information fully obtained from the observational data is introduced. The selfmemorization equation of the heat conduction equation is induced, and the stability and characteristics of which are studied, taking respectively Richardson scheme (RS) and DuFort-Frankel scheme (DS) as the retrospective order p=1. The calculation results indicate that the numerical integral is diffused sometimes in the multi-time model, due to the fact that the memory function is determined by the observational data via the special mathematics arithmetic, which makes the difference scheme designed previously change, and is a smooth time factor in itself.