These folders contain both the files needed to re-run the SAM simulations used in Garrett et al. (20XX), the post-processed data used to calculate the tables/figures, and the scripts to create the post-processed data. This archive does not contain the 2D and 3D output from SAM. As such, some of the scripts will not work since there is no data from which to calculate. Users can either follow the workflow in full, after running SAM from a restart, or recalculate the figures/data using some of the post processed data. The following steps will walk you through both methods.

1) Update the json files in SCRIPTS/namelists so that “dSAM” points to the appropriate SAM output directory. Update the “dWrite” to the folder you wish to save post-processed data (currently ../STATISTICS/).

*** If re-running SAM, follow steps below Otherwise jump to step 6)***

2) Run SAM using the binary files in RESTART, and the grid/namelist files in SAM_SETUP. Convert the com3D binary files to netcdf and ensure the “OUT_3D” folder is located under the base SAM output directory specified in step 1).

3) Run “calc3DVariables.py” in SCRIPTS. This will calculate the MSE, SSE, QVSAT which will be attached to the pre-existing 3D netcdf files and mean profiles which will be saved to

4) Run “calcPerimeterValues.py” in SCRIPTS. This will create netCDF files containing the cloud edge values of MSE for each 3D file. These will be saved in “{dSAM}/OUT_EDGE”. 

5) Run “calcDomainStatistics.py” and “calcPerimeterStatistics.py” in SCRIPTS. These will calculate the mean MSE/SSE values of the domain and cloud edge values, which are found in Fig. 3. 

*** Start here if not re-running SAM. Data is already found in STATISTICS***

6) Run “calc2Dhistogram.py” in SCRIPTS. This will calculate the convective stream function for fig. 2. Output will be saved to “STATISTICS”.

7) Run the “EASE_Scaling.ipynb” and “pauluisFigs.ipynb” notebooks in SCRIPTS. This will create the tables and figures used in the paper.

