Previous studies have shown that numerical simulations of the Madden-Julian oscillation (MJO) are very sensitive to the vertical distribution of diabatic heating. Since atmospheric diabatic heating is generally difficult to estimate, the vertical diabatic heating profile associated with the MJO is not well known. Judged by its propagation properties and spatial structure, the MJO is reasonably well simulated by the nine-layer Goddard Laboratory for Atmospheres (GLA) general circulation model. Although only a simulation the model MJO may provide an indication of the vertical diabatic heating profile associated with the real oscillation. The diabatic heating structure of the model MJO is illustrated with composite charts made for those times when this low-frequency mode reaches its maximum and minimum amplitudes. These composite charts compare the vertically integrated diabatic heating with potential functions, the vertical distribution of diabatic heating with the east-west mass flux function in the tropics, and the vertical profiles of diabatic heating at the centers of maximum and minimum MJO amplitude. Three interesting features of the model MJO's diabatic heating are revealed: (1) the maximum heating rate of this low-frequency mode is located over the Asian monsoon region and its amplitude is about a half of the maximum value of the seasonal mean heating rate in this region, (2) the vertical diabatic heating rate profile has a maximum at 500 mbar and resembles the seasonal mean total heating profile, and (3) the total diabatic heating is for the most part composed of the latent heat released by cumulus convection. ¿ American Geophysical Union 1993