Additional software and data for "The glass ceiling convective regime and the origin and diversity of coronae on Venus"
Creators: name: Kerr, Madeleine
name: Stegman, Dave
Publisher: name: Zenodo
Year: 2025
Subjects: Venus
Mantle convection
planetary science

description: This repository stores additional software and data for the preprint "The glass ceiling convective regime and the origin and diversity of coronae on Venus" Kerr, MC; Stegman DR; Smrekar SE; Adams AC. 2025

The 2D cylindrical mantle convection code APSECT was run on the NERSC computer Perlmutter with the following dependencies:

"------------------------------------------------------------------------------- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.--     . version 2.4.0 (aspect-2.4, 9232772c3)--     . using deal.II 9.3.3--     .       with 64 bit indices and vectorization level 2 (256 bits)--     . using Trilinos 12.18.1--     . using p4est 2.3.2--     . running in OPTIMIZED mode--     . running with 128 MPI processes-----------------------------------------------------------------------------

Loading shared library </pscratch/sd/m/mkerr/code_ASPECT/aspect/benchmarks/entropy_adiabat/plugins/libentropy_plugins.so>

"

In the directory titled 'ASPECTparameter-files' there are ASPECT parameter files and input files used with APSECT version 2.4 and the entropy adiabat benchmark plugin which is configured and compiled in addition to ASPECT. These parameter files are adapted to Venus geometry, gravity, and viscosity structure from the parameter files used in "An entropy method for geodynamic modeling of phase transitions: Capturing sharp and broad transitions in a multi-phase assemblage" by Dannberg et al. 2022.

Additionally there is a table generated by HeFESTo in the 'ASPECTparameter-files/additional_inputs/' directory called "pyrtable1_s.aspect.txt" which is the entropy-pressure table which defined the thermodynamic properties of the composition. This table is the same one used in "An entropy method for geodynamic modeling of phase transitions: Capturing sharp and broad transitions in a multi-phase assemblage" by Dannberg et al. 2022

In the folder 'analysis', there is a series of modules and scripts which are used to analyze the output to the ASPECT model data and to compare the models to Venus observations. The power spectra data of the surface deformation along the Baltis Vallis canali, as published by Conrad and Nimmo (2023), are available on their Zenodo and are here as well. 

To run scripts and use these modules, the following python packages are needed:

Python v3.12.2 with packages Pyvista v0.43.7, Matplotlib v3.8.4, Numpy v1.26.4, Crameri scientific computing color maps v1.19, and Scipy v1.2.0. or newer.

To use with model output stored in a directory with the general structure 'global_dir/model_name/output/solution', start by modifying the global variables python file and run the 'write_heatflux_files.py' file first which generates new files that simplify further computation from the initially heavy dataset.

In the directory "plotted_data" is the tabular version of the data in the plots, excluding colormaps which can be made simply with the HeFESTo table provided and the scripts in 'analysis', or the visualizations, which can be made with pyvista or Paraview with the sample outputs in the ASPECTparameter-files directory or by running ASPECT with the given parameter file for a more complete output.

For questions, contact Madeleine Kerr: mkerr@ucsd.edu or Dave Stegman at dstegman@ucsd.edu