Computation of longwave radiative flux and vertical heating rate with 4A-Flux v1.0 as an integral part of the radiative transfer code 4A/OP v1.5

Based on advanced spectroscopic databases, line-by-line and layer-by-layer radiative transfer codes numerically solve the radiative transfer equation with very high accuracy. Taking advantage of its pre-calculated optical depth lookup table, the fast and accurate radiative transfer model Automatized Atmospheric Absorption Atlas OP-erational (4A/OP) calculates the transmission and radiance spectra for a user-defined layered atmospheric model. Here, we present a module called 4A-Flux, which is developed and implemented into 4A/OP in order to include the calculation of the clear-sky longwave radiative flux profiles and heating rate profiles at a very high spectral resolution. Calculations are performed under the assumption of local thermodynamic equilibrium, a plane-parallel atmosphere, and specular reflection on the surface. The computation takes advantage of pre-tabulated exponential integral functions that are used instead of a classic angular quadrature. Furthermore, the sub-layer variation of the Planck function is implemented to better represent the emission of layers with a high optical depth. Thanks to the implementation of 4A-Flux, 4A/OP models have participated in the Radiative Forcing Model Intercomparison Project (RFMIP-IRF) along with other state-of-theart radiative transfer models. 4A/OP hemispheric flux profiles are compared to other models over the 1800 representative atmospheric situations of RFMIP, yielding an outgoing longwave radiation (OLR) mean difference between 4A/OP and other models of -0.148 Wm(-2) and a standard deviation of 0.218 Wm(-2), showing a good agreement between 4A/OP and other models. 4A/OP is applied to the Thermodynamic Initial Guess Retrieval (TIGR) atmospheric database to analyze the response of the OLR and vertical heating rate to several perturbations of temperature or gas concentration. This work shows that 4A/OP with 4A-Flux module can successfully be used to simulate accurate flux and heating rate profiles and provide useful sensitivity studies including sensitivities to minor trace gases such as HFC134a, HCFC22, and CFC113. We also highlight the interest for the modeling community to extend intercomparison between models to comparisons between spectroscopic databases and modeling to improve the confidence in model simulations.