Improving upon our purely dynamical work, we present three-dimensional simulations of the atmospheric circulation on Earth-like (exo)planets and hot Jupiters using the GFDL-Princeton Flexible Modeling System (FMS). As the first steps away from the dynamical benchmarks of Heng, Menou & Phillipps (2011), we add dual-band radiative transfer and dry convective adjustment schemes to our computational setup. Our treatment of radiative transfer assumes stellar irradiation to peak at a wavelength shorter than and distinct from that at which the exoplanet re-emits radiation ("shortwave" versus "longwave"), and also uses a two-stream approximation. Convection is mimicked by adjusting unstable lapse rates to the dry adiabat. The bottom of the atmosphere is bounded by a uniform slab with a finite thermal inertia. For our models of hot Jupiters, we include an analytical formalism for calculating temperature-pressure profiles, in radiative equilibrium, which accounts for the effect of collision-induced absorption via a single parameter. We discuss our results within the context of: the predicted temperature-pressure profiles and the absence/presence of a temperature inversion; the possible maintenance, via atmospheric circulation, of the putative high-altitude, shortwave absorber expected to produce these inversions; the angular/temporal offset of the hot spot from the substellar point, its robustness to our ignorance of hyperviscosity and hence its utility in distinguishing between different hot Jovian atmospheres; and various zonal-mean flow quantities. Our work bridges the gap between three-dimensional simulations which are purely dynamical and those which incorporate multi-band radiative transfer, thus contributing to the construction of a required hierarchy of three-dimensional theoretical models.

The complexity of atmospheric modelling and its inherent non-linearity, together with the limited amount of data of exoplanets available, motivate model intercomparisons and benchmark tests. In the geophysical community, the Held-Suarez test is a standard benchmark for comparing dynamical core simulations of the Earth's atmosphere with different solvers, based on statistically-averaged flow quantities. In the present study, we perform analogues of the Held-Suarez test for tidally-locked exoplanets with the GFDL-Princeton Flexible Modeling System (FMS) by subjecting both the spectral and finite difference dynamical cores to a suite of tests, including the standard benchmark for Earth, a hypothetical tidally-locked Earth, a "shallow" hot Jupiter model and a "deep" model of HD 209458b. We find qualitative and quantitative agreement between the solvers for the Earth, tidally-locked Earth and shallow hot Jupiter benchmarks, but the agreement is less than satisfactory for the deep model of HD 209458b. Further investigation reveals that closer agreement may be attained by arbitrarily adjusting the values of the horizontal dissipation parameters in the two solvers, but it remains the case that the magnitude of the horizontal dissipation is not easily specified from first principles. Irrespective of radiative transfer or chemical composition considerations, our study points to limitations in our ability to accurately model hot Jupiter atmospheres with meteorological solvers at the level of ten percent for the temperature field and several tens of percent for the velocity field. Direct wind measurements should thus be particularly constraining for the models. Our suite of benchmark tests also provides a reference point for researchers wishing to adapt their codes to study the atmospheric circulation regimes of tidally-locked Earths/Neptunes/Jupiters.