This directory contains the code needed to reproduce the findings and figures of "Effect of Uncertainty in Water Vapor Continuum Absorption on Radiative Forcing, Longwave Feedback and Climate Sensitivity". To reproduce the data published with the manuscript: spectral_olr.nc: - calculate spectral OLR using konrad and ARTS for different experiments (run_konrad_continuum.py). The variables are the varied continuum species, scaling factor, and surface temperature to calculate the spectral longwave feedback, as well as perturbations of the surface temperature (+/- 1K) to calculate the surface/atmospheric feedback, and of the CO2 concentration (doubling) to calculate the radiative forcing. - reformat and save data (convert_data.py) opacity_emission_level.nc: - calculate the optical depth of each vertical layer and each considered absorption species for different surface temperatures (calc_optical_depth.py) - vertically integrate optical depth (integrate_optical_depth.py) - calculate the emission level for each temperature (calc_emission_level.py) - merge the emission levels into one array (merge_emission_level.py) - reformat and save data (convert_data.py) To reproduce the figures: Figure 1: - calculate the absorption cross-section of the self continuum (calc_cross_section.py) - plot opacity and cross sections (plot_Fig1.py) Figure 2: - calculate spectrally resolved feedbacks, forcing and plot them (plot_Fig2.py) Figure 3: - calculate spectrally integrated feedbacks, forcing, as well as climate sensitivity and plot them (plot_Fig3.py) Figure 4: - calculate emission fraction (calc_emission_fraction.py) - read opacity, emission fraction, spectral OLR, and spectral feedbacks and plot them (plot_Fig4.py)

This work was financially supported by the US National Science Foundation (award AGS-1916908) and by NOAA (award NA20OAR4310375).