This repository contains codes to reproduce the results in: Rudolph, M.L., Manga, M., Walker, M., and Rhoden, A.R. (in review). Cooling crusts create concomittant cryovolcanic conduits.
Thermal and stress evolution of icy moons, with application to Europa and Enceladus. Requirements:
- MATLAB 2020a or later due to the use of tiled layouts.
- Fabio Crameri's scientific colormaps
- Parallel computing toolbox
Description of files:
- The equations solved by the code are derived in the Supplemental Material:
2021GL094421_Supplemental_Material.pdf
- The main driver program that sweeps parameter space is
sweep_parameter_space.m
. This will produce about 16GB of output as .mat files. - The file
thickening_ice_shell.m
solves for the stresses in an initially stress-free thickening ice shell with no heat input. - Helper subroutines are located in the subdirectory
core/
- Files in the
benchmarks/
subdirectory can be used to reproduce the results in Nimmo (2004).