Supplementary code for "Effect of Uncertainty in Water Vapor Continuum Absorption on CO2 Forcing, Longwave Feedback, and Climate Sensitivity"

Abstract

<h3>This directory contains the code needed to reproduce the findings and figures of "Effect of Uncertainty in Water Vapor Continuum Absorption on CO2 Forcing, Longwave Feedback, and Climate Sensitivity".</h3><p> </p><h4>To reproduce the data published with the manuscript:</h4><h4>spectral_olr.nc:</h4><ul><li>unzip modified_continuum_input_files.zip in the folder arts-cat-data/model/mt_ckd_4.0/</li><li>calculate spectral OLR using konrad and ARTS for different experiments (run_konrad_continuum.py). The variables are the continuum input file, 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.</li><li>reformat and save data (convert_data.py)</li></ul><h4>opacity_emission_level.nc:</h4><ul><li>calculate the optical depth of each vertical layer and each considered absorption species for different surface temperatures (calc_optical_depth.py)</li><li>vertically integrate optical depth (integrate_optical_depth.py)</li><li>calculate the emission level for each temperature (calc_emission_level.py)</li><li>merge the emission levels into one array (merge_emission_level.py)</li><li>reformat and save data (convert_data.py)</li></ul><h4>continuum_ref_abs_coef.nc</h4><ul><li>calculate absorption coefficient (proportional to optical depth) for reference conditions and save them (calc_abs_coef.py)</li><li>reformat and save data (convert_data.py)</li></ul><h4>modified_continuum_input_files.zip</h4><ul><li>read original data on self and foreign continuum from MT_CKD 4.0 and scale them (rescale_continuum.py)</li><li>the modified continuum files have to be located in the folder arts-cat-data/model/mt_ckd_4.0/ to run run_konrad_continuum.py</li></ul><h4> </h4><h4>To reproduce the figures:</h4><h4>Figure 1 & 2:</h4><ul><li>calculate the absorption cross-section of the self continuum (calc_cross_section.py)</li><li>plot opacity and cross sections (plot_Fig1_Fig2.py)</li></ul><h4>Figure 3 & 4:</h4><ul><li>calculate spectrally resolved and integrated feedbacks, forcing and plot them (plot_Fig3_4.py)</li></ul><h4>Figure 5:</h4><ul><li>calculate emission fraction (calc_emission_fraction.py)</li><li>read opacity, emission fraction, spectral OLR, and spectral feedbacks and plot them (plot_Fig5.py)</li></ul><h4>Figure 6:</h4><ul><li>calculate absorption coefficient (proportional to optical depth) for reference conditions and save them (calc_abs_coef.py)</li><li>read absorption coefficients, calculate scaling factors and plot them (plot_Fig6.py)</li></ul><p>This work was financially supported by the US National Science Foundation (award AGS-1916908) and by NOAA (award NA20OAR4310375).</p&gt