This repository contains the data, scripts, SAM codes and files required to reproduce the results of the manuscript "The Unreasonable Efficiency of Total Rain Evaporation Removal in Triggering Convective Self-Aggregation" submitted to the Geophysical Research Letters (GRL). Brief description of project This project aims to examine the impact of rain evaporation removal or reduction in the planetary boundary layer (PBL) on convective self aggregation (CSA). Non-rotating radiative-convective equilibrium (RCE) simulations were conducted with the System for Atmospheric Modeling (SAM) cloud resolving model. Rain evaporation in the lowest 1 km was progressively reduced and the effect on CSA was investigated. The physical processes underlying this type of aggregation (referred to in the manuscript as no-evaporation CSA, or NE-CSA) were analyzed and described. The default SAM code base (version 6.10.8) can be downloaded from here: http://rossby.msrc.sunysb.edu/~marat/SAM.html Repository structure /SAM6.10.8 directory: contains the SAM fortran code and run time files required to conduct the simulations described in the manuscript /SAM6.10.8/SRC: contains the modified fortran code to run the experiments where rain evaporation was removed or reduced /SAM6.10.8/RCE_MM: contains the run time files required to run the RCE experiments /data directory: contains the model output data, in csv format /data/pw: contains data for the spatial snapshot of precipitable water /data/pw_range_timeseries: contains timeseries data for the PW range (moist minus dry patches) /data/cold_pools: contains data for the virtual potential temperature anomaly on 1st model level as proxy for cold pools /data/MSE_streamfunction: contains data for plotting the MSE and streamfunction plot /data/thetav_timeseries: contains the timeseries data for virtual potential temperature, potential temperature and water vapor anomalies (deviation from domain mean) for buoyancy analyses /data/qrad_profile_end: contains data for the vertical profiles of the steady state radiative cooling rates /data/qrad_profile_start: contains data for the vertical profiles of the early stage (first-10 days average) radiative cooling rates /data/w_timeseries: contains timeseries data for the vertical velocity in the dry regions, averaged over lowest 1, 2, 3 and 5 km /data/agg_index_timeseries: contains timeseries data for the aggregation index, measured using two indices: variance of precipitable water divided by mean ("org_pw"), and variance of column relative humidity ("crh_var") /data/thetav_profile: contains data for the vertical profiles of virtual potential temperature, potential temperature and water vapor anomalies averaged over five separate days /data/dry_subs_timeseries: contains data for plotting figures for the dry-subsidence feedback analyses (time-height arrays of w, \(q_{v}\), and \(q_{v}\) vertical flux) /data/Tq_moist_dry_timeseries: contains timeseries data of T and \(q_{v}\) of the dry and moist patches, averaged over the PBL (lowest 1 km) /scripts directory: contains the scripts required to plot all figures in the manuscript plot_cold_pools_pw.py: script to plot Figure 1 (RCE mean state PW and cold pools for the HomoRad and HomoFlx cases) plot_MSE_streamfunction.py: script to plot Figure 2 (MSE and streamfunction averaged over final 5 days) plot_thetav_timeseries.py: script to plot Figure 3 (\(\theta^{'}_{v}\), \(\theta^{'}\)and \(q_{v}^{'}\) timeseries for aggregated and scattered cases) plot_qrad_profiles_pw_w_timeseries.py: script to plot Figure 4 (early stage radiative cooling profiles, PW range and \(\overline{w}_{dry,PBL}\) timeseries) SI_plot_agg_index_timeseries.py: script to plot Figure S1 (aggregation index timeseries) SI_plot_cold_pools_pw_IntRadIntFlx.py: script to plot Figure S2 (RCE mean state PW and cold pools for the various IntRad and IntFlx cases) SI_plot_cold_pools_pw_various_PBL_heights.py: script to plot Figure S3 (RCE mean state PW and cold pools for various PBL heights) SI_plot_thetav_profiles.py: script to plot Figure S4 (\(\theta^{'}_{v}\), \(\theta^{'}\)and \(q_{v}^{'}\) profiles on five separate days) SI_plot_dry_subsidence_time_height.py: script to plot Figure S5 (time-height arrays of w, \(q_{v}\), and \(q_{v}\) vertical flux) SI_plot_Tq_moist_dry_timeseries.py: script to plot Figure S6 (timeseries of T and \(q_{v}\) of moist and dry patches, averaged over the PBL) SI_plot_thetav_timeseries.py: script to plot Figures S7 and S8 (\(\theta^{'}_{v}\), \(\theta^{'}\)and \(q_{v}^{'}\) timeseries for additional \(\alpha\)values, for aggregated and scattered cases) SI_plot_CDA_RDA_thetav_timeseries.py: script to plot Figure S9 (\(\theta^{'}_{v}\), \(\theta^{'}\)and \(q_{v}^{'}\) timeseries comparison for CDA vs. RDA) /scripts/applications directory: contains the utilities script that include helper functions for the plotting scripts