eet the technical requirements of IMT-2020, new features are captured in 5G channel models such as support of frequencies up to 100 GHz and large bandwidth, 3D modeling, support of large antenna arrays, blockage modeling, and spatial consistency [4,5]. The 3D modeling describes the channel propagation in azimuth and elevation directions between the transmitting and the receiving antennas. It is more complete and accurate rela- tive to the 2D modeling which only considers the propagation characteristics in the azimuth direction. Multi-antenna techniques capable of exploiting the elevation dimension have been developed for LTE since Rel-13 and are considered very important in 5G, which include modeling the elevation angles of departure and arrival, and their correlation with other para- meters [7]. For a base station (access node) equipped with columns of active antenna arrays, vertical analog beamforming can be applied prior to the power amplifier (PA); thus the electrical downtilt can be dynamically adjusted over time. This feature allows the base sta- tion to dynamically adjust its cell coverage depending on the user distribution in the cell, for example, analog beamforming can be directed toward UEs congregated at the same location. Nonetheless, such beamforming is still fundamentally cell specific (i.e., spatial separation of individual UE is not possible) and operate on a long timescale, although time- varying base station antennas are modeled having one or multiple antenna panels, where an antenna panel has one or multiple antenna elements placed vertically, horizontally or in a 2D array within each panel. As a result, 3D channel modeling is required for performance evaluation of full-dimension MIMO. The full-dimension MIMO involves precoding/beam- forming exploiting both horizontal and vertical degrees of freedom on a small timescale in a frequency-selective manner. We start our study with the 3D antenna models in mmWave bands where implementation of large antenna arrays is feasible in the gNB or the UE. Le