Abstract: A high precision orbit model based on numerical integration is introduced. Satellite's instantaneous position and velocity are derived from the position data measured by on-board GPS receiver, and then they are used as initial value to predict the orbit by high precision numerical integration. The six orbital elements which are used to describe the motion of satellite can be transformed to the satellite's position and velocity, and vice versa. The satellite's instantaneous orbital elements at time t₁ can be calculated from satellite's two position vectors and at time t₂ and t₂ respectively, then the satellite's instantaneous position and velocity at time t₂ can be obtained. There are many perturbing terms that have been considered in the satellite's high precision orbit model, including the aspherical gravitations of the earth, the perturbations that arise from the gravitational attraction of the sun and the moon, the solar radiation pressure and the atmospheric drag. The earth gravity model that is used is the EGM-96 model of degree and order of 360. Numerical methods are used to solve the equation of the satellite's motion. The variable order and stepsize method DE/DEABM are used for the numerical integration.This new method are applied to Terra/Aqua-MODIS satellite's orbit calculation. Satellite Terra's position vectors r₁ at time t₁ and r₂ at time t₂ are used to compute the instantaneous velocity vector. The result shows that the departure between the calculated observation and the GPS measurement is very small with the error of instantaneous speed less than±0.1 m/s. Then they are used as initial input to predict one of Terra/MODIS' further 10-minute orbit by applying the high precision orbit model. It shows that the satellite's position and velocity vectors and range which are predicted by this method are all consistent with the measurements by GPS. The following results are gained that the position error is between-20 m and +5 m, the velocity error is less than±0.2 m/s. The same method are used to compute Aqua/MODIS' orbit. The Aqua's position vectors error between the predictions from high precision orbit model and the GPS measurements is less than 40 m, the error of predicted velocity vectors is less than±0.4 m/s.Because of the high precision and the continuous stability in certain time, the high precision orbit model based on numerical integration can be used to verify the reliability of GPS measurement results if the initial input is correct, and substitute the wrong GPS measurement value. Further more, the new method is suitable for FY-3 satellite's orbit computation.