During the boreal summer, the intraseasonal oscillation has the tendency to propagate northward from the Indian Ocean into the land area in South Asia, and to propagate northwestward in the Philippine Sea toward Taiwan and Southern China. This study investigates the ability of seven GCMs, namely, DNM, JMA, SUNYA, NCAR, PNNL, ECHAM, NTU, in simulating this phenomenon. Most of the GCM underestimate the intraseasonal variability in the regions, e.g., the Bay of Bengal and the western North Pacific, where the active intraseasonal fluctuation is observed. Spectral analysis of the outgoing longwave radiation in the above two regions reveals that none of the above models exhibits a spectral peak around 30-60 days as in the observation. None of the models is able to reproduce realistic norhtward propagation. Some of them do not even exhibit any preferred propagation route. The PNNL and ECHAM are the two models that have limited capability in simulating the northward propagation of the intraseasonal oscillation. The SVD analysis was applied to the simulated data to retrieve the most recurrent coupled feature of the OLR and 850-hPa streamfunction. The ECHAM model is the only model that produces the first and second singular vectors resembling those documented in the observation. The pattern however exhibits limited northward propagation in the southern Bay of Bengal, which is much less significant than the observed. The deficiency of the model in simulating realistic surface heat flux in the Indian Ocean and the South Asia may result in the limited capability in simulating the northward propagation.