In the case of the irradiance calculation, we had to first find the horizon angle from each of the >1 million pixels in each of 12 directions sampled evenly from the full 360 degrees, as well as calculate the slope angle and aspect at each pixel. Following our parallelization strategy, these calculations were submitted as independent jobs to the batch scheduler, each of which writing an intermediate file based on the output of the calculation. Then a final “master” job was submitted that waited on those previous jobs to finish, took the intermediate files as inputs, and produced the final irradiance raster we needed as a feature for our regression models. The process for calculating MUHA was similar, though it involved inverting the DEM and calculating “negative horizon angles” as well, since these are indicative of a propensity for wind scouring.