From the results of the ray-tracing simulation of 28 GHz frequency band over downtown of Ottawa and New York University campus, we propose 3D large-scale channel models for urban city which are applicable to a mmWave version of 3D spatial channel model. Due to the change of geographical topology dimension from 2D to 3D, the height information of the transmitter should be reflected to the parameters of the channel model such as line-of-sight probability, Ricean K factor, path loss and shadow fading accordingly. In the simulation results, the line-of-sight probability depends on the height of the transmitter compared to the 2D model. The dual-slope path loss model is still appropriate in 3D channel model. We also propose the linear standard deviation function of the shadow fading, which shows a large standard deviation and an increasing trend with the distance.