MinVel

This program is used to calculate anharmonic p- and s-wave velocity and density for zero-porosity mineral aggregates. It is based on the work of Hacker and Abers (2004, updated in 2016) with additional minerals and optimized for related work. Running the code in either MATLAB or Python requires a mineral physics database, MineralPhysicsDatabase.nc (https://doi.org/10.5066/P9HN170G). Combined with a petrologic database (Sowers and Boyd, 2019) and thermal and geologic models, results from this code are used in the program RockProp.m, which accounts for rock porosity and seismic attenuation. These codes are part of an effort to produce a three dimensional national crustal model for use with seismic hazard studies (https://doi.org/10.5066/P9T96Q67).