Based on the previous work, a new three-dimension mid- and low-latitude theoretical ionospheric model in realistic geomagnetic fields is developed, named Three-Dimension Theoretical Ionospheric Model of the Earth in the Institute of Geology and Geophysics, Chinese Academy of Sciences (TIME3D-IGGCAS). This new model covers the mid- and low-latitude ionosphere and whole plasmasphere. It self-consistently solves the equations of mass continuity, motion and energy of electron and ions to give out the time-dependent three-dimensional structures of the main ionospheric and plasmaspheric parameters in realistic geomagnetic fields, including ion number densities of O + , H + , He + , NO + , O 2 + , N 2 + and electron; electron and ion temperature; and ion velocity vectors. We carry out simulations in March Equinox and in June Solstice, and compare the simulated results with that from IRI empirical model. TIME3D-IGGCAS can well reproduce the main ionospheric features in all simulations. We also simulate the ionospheric differences between different kinds of geomagnetic fields. The results suggest that the geomagnetic field configuration obviously affect the ionospheric plasma density, and the differences between NmF2 in realistic geomagnetic fields and that in tilted dipole fields can be larger than 40%.