Time-lapse video for delay-time kernel

Author: Li-Yu Kan (1) & Li Zhao (2)
(1) Department of Earth Sciences, National Central University, Taoyuan, Taiwan
(2) School of Earth and Space Sciences, Peking University, Beijing, China

A movie showing the Vp and Vs delay-time kernel. Made with mc kernel.


Source depthStation DistanceLink
500 km105 deghttps://sgylab.earth.ncu.edu.tw/kernel_video/PS_dep500_105deg.mp4
100 km150 deghttps://sgylab.earth.ncu.edu.tw/kernel_video/PS_dep100_150deg.mp4
20 km60 deghttps://sgylab.earth.ncu.edu.tw/kernel_video/PS_dep20_60deg.mp4

High-pass 100s filtered (zero phase)

Source depthStation DistanceLink
500 km105 deghttps://sgylab.earth.ncu.edu.tw/kernel_video/PS_dep500_105deg_hp.mp4
500 km
130 deghttps://sgylab.earth.ncu.edu.tw/kernel_video/PS_dep500_130deg_hpR.mp4
100 km150 deghttps://sgylab.earth.ncu.edu.tw/kernel_video/PS_dep100_150deg_hp.mp4
100 km180 deghttps://sgylab.earth.ncu.edu.tw/kernel_video/PS_dep100_180deg_hp.mp4
20 km105 deg(not finished)

Draw the outline of the text in GMT

While we choose a bright color for filling the text, a black outline is usually used to emphasize it (or it is unnoticeable). For plotting text on the canvas, we would immediately come up with pstext. But when we want to depict the outline of the text, there is no intuitive way or independent arguments for this. The -G and -W of pstext are actually for the background textbox rather than the text itself. Here, we provide two methods to achieve this: append =<pen> to the font specification in pstext or use psxy -Sl to plot the outline of the text.

閱讀更多Draw the outline of the text in GMT

Plot SAC with pssac in the cartesian time axis

photo_2019-03-20 18.19.24

After the release of GMT5.3.0, pssac has been included into the GMT package. PSSAC operates a more convenient way than pswiggle or psxy to plot with SAC file (a common file format used in seismology community). We no longer have to convert SAC files into XY files and plot them by sorting them manually. PSSAC provides several types of regular profile like sorting by azimuth, epicentral distance, or user-defined profile through -E arguments. However, there is no direct argument in -E to plot such a figure with a cartesian time of X-axis (or Y if given -Q). In this article, we provide the trick for plotting SAC file in an axis of time.

閱讀更多Plot SAC with pssac in the cartesian time axis