Starting Things Out

I've decided to start this page to document the various problems and solutions I encounter during standard data analysis and mapping for geophysics problems. Most of what I will share is either Python, Matlab, or GMT (the esoteric nature of GMT is particularly vexing and could be an entire  page on its own).

A bit about myself, I am a Research Scientist at the Pacific Northwest Seismic Network, working on seismological and geodetic problems. I received my PhD at Scripps Institution of Oceanography in 2013 and have been at the PNSN ever since. 

Comments

Post a Comment

Popular posts from this blog

Plotting USGS 3d Elevation Data in GMT

Image
This post is about plotting high resolution digital elevation model (DEM) data in GMT for high-quality maps. Every time I need to do this, I end up googling around, wasting a bunch of time, so this post mostly serves as a reference for me. First things first, the USGS provides elevation data at numerous resolutions, the most common being 1 arcsecond (roughly 30 meters) and 1/3 arcsecond (roughly 10 meters). Other resolutions exist, but these are the two you will probably use most. Getting the elevation data can be done through the USGS's National Map web service ( https://viewer.nationalmap.gov/basic/#productGroupSearch ). The process is rather simple, just zoom into the area you care about (you can also enter coordinates) and then under " Elevation Products (3DEP)", click which one you want. After clicking 'Find Products', you can go to the Products tab and toggle on and off the footprints of the different maps. They are done in 1 degree squares, so you will need
Read more

Culling data based on coastlines in GMT

Image
This is a problem that took me way too long to figure out, but it is a simple one line command. The problem arose as I was taking a finite fault solution for the Kaikoura earthquake from a different author and forward modeling the surface displacements onto a set grid of points (in my case, either 10, 20 or 40 km). From there, I wanted to invert for a possible slip model based on this new grid of 'synthetic' GPS stations, therefore, I had to exclude any stations/grid points over water. The solution for this is the gmtselect command, and its used as follows: Since my data was in a separate file (with columns lon, lat, east, north, up), I first start with the awk command within my bash script awk '{print $1, $2, $3, $4, $5}' file.txt I then add on the gmtselect command awk '{print $1, $2, $3, $4, $5}' file.txt | gmtselect -R -Di -Ns/k/s/k/s > subset.txt See, easy. The culled solutions are now in the file 'subset.txt' with the same column order
Read more