Scripts and data for "A Strong Constraint on Radiative Forcing of Well-mixed Greenhouse Gases"

All-sky Radiative Fluxes Computed by GRTcode using ERA5 Reanalysis and CMIP7 Input4MIPs as Inputs for 2001 to 2024. These results are spectrally integrated as longwave broadband fluxes for assessible data size. Spectrally resolved results are avaiable upon request to [e-mail]. GRTcode is available at https://github.com/fengzydy/GRTCODE.git Using annual-mean, monthly-mean surface temperature, humidity, ozone concentration, and cloud conditions from the ERA5 reanalysis dataset (Hersbach et al., 2020), a set of global-scale line-by-line experiments perturbing the concentrations of well-mixed greenhouse gases was conducted for the period 2001–2024. These experiments include: Control: A control experiment using a time series of global-mean, annual-mean concentrations of WMGHGs (Nicolls et al., 2025, as included in control.nc). x_PI: Identical to Control, except that an individual WMGHG species, x, is held fixed at its 1850 concentration (as included in PI.nc). x_2xPI, x_3xPI, x_4xPI: Similar to x_PI, except that concentration x is two, three and four times of 1850 concentration. x_4x: Similar to x_PI, except that concentration x is four times of its annual-mean concentration. Data under 'hres' folder is at 2.5x2.5 degree resolution, and is used to construct the regression coefficients. Other data is at 7.5x7.5 degree resolution and is used to evaluate the regression model. The scripts for using the regression model are provided as 'OLR_2_forcing' for matlab and python environment. Scripts for generating regression coefficents, key figures, and IRF data are included as plot_ERA5LBL.m Data used for generating Figures 1(a), 2, and 4 directly is included in a .xlsx file.