We present results of our pulsar population synthesis of normal pulsars from the Galactic disk using our previously developed Monte-Carlo code. From our studies of observed radio pulsars that have clearly identifiable core and cone components, in which we fit the polarization sweep as well as the pulse profiles in order to constrain the viewing geometry, we develop a model describing the luminosity and ratio of radio core-to-cone peak fluxes. In this model, short period pulsars are more cone-dominated. We explore models of neutron star evolution with and without magnetic field decay, and with different initial period distributions. We present preliminary results including simulated population statistics that are compared with the observed radio pulsar population. The evolved neutron star populations resulting from this simulation can be used to model distributions of gamma-ray pulsars for comparison to Fermi results.