Non-equilibrium fractionation during ice cloud formation in iCAM5: evaluating the common parameterization of supersaturation as a linear function of temperature

This archive includes data and python scripts to create the figures in Duetsch, M., Blossey, P. N., Steig, E. J., and Nusbaumer, J. M. (2019). Non-equilibrium fractionation during ice cloud formation in iCAM5: evaluating the common parameterization of supersaturation as a linear function of temperature. Submitted to J. Adv. Model. Earth Sy. <strong>Model</strong> Model output data are saved in model.tar.gz For all simulations, monthly averages of the following variables are saved in h0 files (used for Figures 3, 6, 8, 9):<br> T: temperature (Si_real simulations only)<br> PS: surface pressure (Si_real simulations only)<br> U: zonal wind (Si_real simulations only)<br> V: meridional wind (Si_real simulations only)<br> LANDFRAC: fraction of surface area covered by land<br> PRECT_H2O: total precipitation rate for H2O<br> PRECT_HDO: total precipitation rate for HDO<br> PRECT_H218O: total precipitation rate for H218O<br> H2OV: H2O mixing ratio for vapor<br> HDOV: HDO mixing ratio for vapor<br> H218OV: H218O mixing ratio for vapor<br> H2OI: H2O mixing ratio for cloud ice<br> HDOI: HDO mixing ratio for cloud ice<br> H218OI: H218O mixing ratio for cloud ice<br> H2OL: H2O mixing ratio for cloud liquid<br> HDOL: HDO mixing ratio for cloud liquid<br> H218OL: H218O mixing ratio for cloud liquid For the control and Si_real microphysical sensitivity simulations, 6-hourly averages of the following variables are saved in h1 files (used for Figures 4, 5, 7):<br> PRECL_H2O: Large-scale precipitation rate for H2O (Si_real simulations only)<br> PRECL_HDO: Large-scale precipitation rate for HDO<br> PRECL_H218O: Large-scale precipitation rate for H218O<br> PRECL_SAT: Large-scale precipitation rate for Si tracer (Si_real simulations only)<br> PRECL_TMP: Large-scale precipitation rate for T tracer (Si_real simulations only)<br> PRECL_RVD: Large-scale precipitation rate for R^D tracer<br> PRECL_RVO: Large-scale precipitation rate for R^18O tracer To limit the size of the data set, only the first 5 years of the simulations are saved. <strong>Measurements</strong> Ice core measurements for present-day climate and last glacial maximum are saved in icecores_PD.txt and icecores_LGM.txt, respectively. LGM values are averaged from 19000 BCE to 16000 BCE, PD values are averaged from 1000 BCE to 2000 CE. Antarctic surface snow measurements from Masson-Delmotte et al. (2008) are available at https://doi.org/10.1594/PANGAEA.681697 <strong>Python scripts</strong> Scripts run with python 3.6 Required packages:<br> - numpy<br> - matplotlib<br> - netCDF4<br> - datetime<br> - Basemap<br> - copy<br> - intergrid<br> - cmocean