P-wave velocity and density models of Montserrat

P-wave velocity model and density model obtained by joint inversion of seismic traveltimes and gravity anomalies as described in Paulatto et al. (2019).The traveltime data come from the SEA-CALIPSO experiment and are described in Paulatto et al. (2012). The gravity data are a compilation of land and marine gravity data. The code used for the inversion (jif3D, https://sourceforge.net/projects/jif3d/) is described in Moorkamp et al. (2010). The authors welcome reuse of the models and data. For any questions contact Michele Paulatto ([e-mail]).<br><b>Included files:</b>monji_ttcx_1d0dd.nc = traveltime datasetmonji_ga5_grad_f500.nc = gravity anomaly datasetresult.tomo.inv.nc = output slowness fileresult.grav.inv.nc = ourput density anomaly fileresult.dens.inv.nc = output density file (density anomaly + reference 1d density model)paperzsections_2_4_6.pdf = plots of horizontal cross-sectionspvsection_23_23.5.pdf = plots of vertical cross-sections<br><br><b>Coordinates and units</b>The model coordinates are defined over a local Cartesian coordinate system with origin at longitude -62.4 degrees, latitude 16.5 degrees. Easting, northing and depth are in units of metres from the origin point. Depth is referenced to an elevation of 2000 m above mean sea level and is positive downward (zero Depth corresponds to 2000 m elevation). The Slowness field is in units of s/m and the density is in units of g/cm³.<br>Northing bounds: 0 / 50000Easting bounds: 0 / 45000Depth bounds: 0 / 14500 (-2000 / 12500 )<br><b>Reading the models</b>The contents can be loaded in Python using the netcdf4 library. See example code snippet below.<br> from netCDF4 import Dataseta = Dataset('result.tomo.inv.nc', mode='r')northing = a.variables['Northing'][:]easting = a.variables['Easting'][:]depth = a.variables['Depth'][:]slowness = a.variables[Slowness][:,:,:]a.close()<br>b = Dataset('result.dens.inv.nc', mode='r')northing = b.variables['Northing'][:]easting = b.variables['Easting'][:]depth = b.variables['Depth'][:]density = b.variables[Density][:,:,:]b.close()<br>c = Dataset('result.grav.inv.nc', mode='r')northing = c.variables['Northing'][:]easting = c.variables['Easting'][:]depth = c.variables['Depth'][:]density_a = c.variables[Density][:,:,:]c.close()<br>