Reproducibility Repository for "Microphysical parameter choices modulate ice content and relative humidity in the outflow of a warm conveyor belt"

This repository contains the code, data, and instructions to reproduce the results from our 2025 publication: "Microphysical parameter choices modulate ice content and relative humidity in the outflow of a warm conveyor belt"DOI (preprint): 10.5194/egusphere-2025-1816 The code (without data) will also available at my GitHub page -------------------------------------------------------------------- Instructions 1. Unzip the archive2. Navigate to the CODE/ directory3. Make sure you have Julia installed4. Install required packages by running: julia --project=. install_packages.jl5. Open the Julia REPL: julia --project=.6. Recreate the paper plots by running: include("create_plots.jl") plot_Fig2(true) plot_Fig3(true) ... continue with Fig4 ... Fig16, FigA1, SI1 ... SI7 Note: Fig. 1 is not included here because it requires access to ERA5 data. --------------------------------------------------------------------Repository Contents Julia Scripts:- general_plots.jl Plotting building block functions- plotting_functions.jl Definitions of plotting functions- create_plots.jl Generates paper figures (calls functions above)- do_analysis_*.jl Additional trajectory data analysis, produces filtered netCDF files- model.jl Functions for the RF models used in the paper- calcs.jl Calculation utilities, constants, module imports- statistical_functions.jl Statistical measure functions- thermodynamic_functions.jl Thermodynamic calculations- use_config.jl Loads useful packages for Julia REPL work Other Files:- cartopy_background.py Python script for creating background maps with borders- slurmfiles/ Output directory for cluster jobs run with SLURM --------------------------------------------------------------------Data The DATA/ directory contains trajectory data in the sens_exp_ID0XX/ directories, as well as the processed data: - PPE_vars_all.nc- PPE_vars_end_hours.nc- PPE_vars_accum.nc- PPE_vars_accum_hours.nc- PPE_vars_t600.nc These files contain multiple variables (e.g., pressure of glaciation, summary statistics) and are produced by the following scripts: Script Output file-------------------------------------------------do_analysis_all.jl PPE_vars_all.ncdo_analysis_after_hours.jl PPE_vars_end_hours.ncdo_analysis_accum.jl PPE_vars_accum.ncdo_analysis_accum_hours.jl PPE_vars_accum_hours.ncdo_analysis_t600.jl PPE_vars_t600.nc Running analysis: To recreate a netCDF file like PPE_vars_all.nc (may require large RAM): julia --project=. do_analysis_all.jl For multithreading: julia --project=. --threads X do_analysis_all.jl (replace X with number of threads; note that some lines may need to be uncommented, and stability is not guaranteed). --------------------------------------------------------------------Additional Notes - params.txt Contains information on the perturbed parameters.- Additional data from Oertel (2025) relevant to this study is available at: https://doi.org/10.35097/ecgs4f56mp3ymjmt Acknowledgements We thank all four reviewers for their constructive feedback which greatly improved the interpretation of our results and the overall quality of the paper. The PPE was built within subproject B8 of the Transregional Collaborative Research Center SFB/TRR 165 “Waves to Weather” (https://www.wavestoweather. de) funded by the German Research Foundation (DFG). AO is supported by the Italia–Deutschland science services network in weather and climate (IDEA-S4S; INVACODA, grant no. 4823IDEAP6). This Italian–German research network of universities, research institutes, and the Deutscher Wetterdienst is funded by the BMDV (Federal Ministry of Transport and Digital Infrastructure). The authors acknowledge support by the state of Baden-Württemberg through bwHPC. The PPE simulations were carried out on the supercomputer HoreKa at Karlsruhe Institute of Technology, Karlsruhe, which is funded by the Ministry of Science, Research and the Arts Baden-Württemberg, Germany, and the German Federal Ministry of Education and Research. This work was also funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)– TRR 301– Project-ID 428312742: “The tropopause region in a changing atmosphere”, sub-project B08 coordinated by Annette Miltenberger.