Arctic sea ice box model with 2-layer atmosphere

toy_model.py contains the full model equations, physical constants used, and integration scheme. 2D_sweep_revisions.py is a driver file that calls functions from toy_model to sweep over two different model parameters simultaneously (i.e., testing all pairwise parameter combinations) using multi-threading. The default parameter values listed in 2D_sweep_revisions can also be used to run a single model simulation. All constants have been updated to match those used in the final version of the accepted manuscript.