Reproducibility Repository for "Effects of Model Grid Spacing for Warm Conveyor Belt (WCB) Moisture Transport into the Upper Troposphere and Lower Stratosphere (UTLS) — Parts I and II"

This repository contains the code used to reproduce the results presented in the following publications submitted to ACP:“Effects of Model Grid Spacing for Warm Conveyor Belt (WCB) Moisture Transport into the Upper Troposphere and Lower Stratosphere (UTLS)” — Part I: Lagrangian Perspective, and Part II: Eulerian Perspective. Data Access The numerical weather prediction (NWP) data required for these analyses are not publicly available due to their large size. To reproduce the results, you will need:- ICON model output for both the global and nested simulations, and- The full trajectory dataset used in the studies. Alternatively: the finished analysis data, zipped size ca. 100GB. Please contact the authors to obtain access to these data. --- Prerequisites You need working Julia and Python environments. Navigate to the CODE/ directory. 1. Install all required Julia packages: julia --project=. install_packages.jl2. Install all required Python packages: python install_packages.py 3. Make sure you have CDO (Climate Data Operators) available as a module and loaded. On the MOGON-NHR Cluster: module load data/CDO/1.9.8-intel-2019b Setup 1. Open `PATH_DEFINITIONS.jl` and modify all directory paths to point to your local data and working directories.2. Navigate to CODE/EXECUTE_WCB_SELECTION/ and to merge model output and create WCB files (on a HPC cluster): sbatch DO_WCB_SELECTION.slurm Make sure to modify the SLURM script according to your cluster. This script executes do_WCB_selection.jl which calls WCB_selection.jl to generate all required intermediate files for subsequent trajectory and WCB analyses. You can also run each line from DO_WCB_SELECTION.slurm in the command line without SLURM if you have enough memory (~40GB) --- Recreating Figures for Part I To reproduce all figures and analysis presented in Part I, run:julia --project=. PART_I/make_plots.jl--- Recreating Figures for Part II To reproduce all figures and analysis presented in Part II, execute the following steps in order: 1. Create the pressure-level WCB masks: python PART_II/make_plev_WCB_mask.py2. Build Eulerian masks from trajectory data: julia --project=. PART_II/create_masks.jl3. Perform the pressure-level analysis: julia --project=. PART_II/do_plev_analysis.jl4. Generate all figures: julia --project=. PART_II/make_plots.jl Contact For data requests, questions, or reproducibility inquiries, please contact:Cornelis Schwenk ([e-mail], private: [e-mail])Annette Miltenberger ([e-mail])