With three years of sky-survey data from the Fermi LAT, significant gamma-ray pulsations have been detected from ~40 millisecond pulsars (MSPs) using radio timing solutions from across the globe. We have fit the radio and gamma-ray pulse profiles of these MSPs using geometric versions of slot gap and outer gap gamma-ray emission models and radio cone and core models. For MSPs with radio and gamma-ray peaks aligned in phase we also explore low-altitude slot gap gamma-ray models and caustic radio models. The best-fit parameters provide constraints on the viewing geometries and emission sites. While the exact pulsar magnetospheric geometry is unknown, we can use the increased number of known gamma-ray MSPs to look for significant trends in the population which average over these uncertainties. Portions of this research performed at NRL are sponsored by NASA DPR S-15633-Y.