This archive contains codes and data files supporting the study by

Maia, J.S., Wieczorek, M.A. and Plesa, A.-C. entitled

"The Mantle Viscosity Structure of Venus"

submitted to Geophysical Research Letters.

It includes 2 directories:

-----
data
-----
admitcor.txt
Global and localized admittance and correlation curves for Venus as presented in Figure 1(b) of the manuscript.

The columns in the file correspond to:
SH_degree - spherical harmonic degree
global_correlation - spectral global correlation between gravity and topography data
global_admittance - spectral global admittance between gravity and topography data (mgal/km)
localized_correlation - spectral localized correlation between gravity and topography data
localized_admittance - spectral localized admittance between gravity and topography data (mgal/km)

=================
multitaper_localization.grd
netcdf file of multitaper localization grid.

=================
samples_nominal.txt
samples_freeslip.txt
samples_uniformload.txt
Samples from the posterior probability distribution of the mantle viscosity profile for the three inversion scenarios as presented in Figure 3.

=================
density_mass-sheet.grd
density_mass-sheet.sh
Distribution of mantle density anomalies:
- grid (netcdf) that can be directly plotted using, e.g., pygmt
- spherical harmonic coefficients that can be directly loaded into pyshtools


------------
codes
------------
create_figures.py

Create the most relevant figures from our study using the data made available in the 'data' directory

Package dependencies:
- numpy
- matplotlib
- pygmt: plot maps https://www.pygmt.org/latest/
- cmcrameri: scientific colormaps https://pypi.org/project/cmcrameri/

==================
Python scripts computed theoretical gravity and topography using the dynamic loading model
by Richards and Hager (1984) and Hager and Clayton(1989)

dynamic_loading.py: computes and saves to tex files the spherical harmonic coeffients of the predicted gravity (m/s2), topography (m), and mass-sheet mantle density anomalies (kg/m2)

dynamic_loading_functions.py: a collection of functions used to compute the model

multitaper_localization_functions.py: a collection of functions to load the data and localize it
using the multitaper approch by Wieczorek and Simons (2005,2007)


Package dependencies:
- numpy
- scipy
- sympy: https://www.sympy.org/
- pyshtools: https://shtools.github.io/SHTOOLS/
