1#!/usr/bin/env bash 2# 3# Check full moment tensor, especially with big isotropic components. 4# 5# The data files and original script are provided by Carl Tape, 6# available from https://github.com/carltape/compearth, licensed under MIT license. 7# 8# Also see https://github.com/GenericMappingTools/gmt/issues/661. 9# 10ps=full_moment_tensor.ps 11gmt set PS_MEDIA 21.5ix9.5i MAP_TICK_LENGTH 0.0c 12 13gmt psxy -JH0/2.8i -R-30/30/-90/90 -T -K -P -X0.5i -Y0.5i > $ps 14for dfile in fullmt_ipts1_iref1 \ 15 fullmt_ipts1_iref2 \ 16 fullmt_ipts1_iref3 \ 17 fullmt_ipts1_iref4 \ 18 fullmt_ipts1_iref5 \ 19 fullmt_ipts2_iref3; do 20 gmt psbasemap -J -R -Bxa10f5g10 -Bya10f5g10 -Bwesn+g200 -K -O >> $ps 21 gmt psmeca ${dfile} -J -R -Sm0.45i/8p -L0.5p,0/0/0 -G255/0/0 -N -K -O >> $ps 22 gmt psxy -J -R -T -K -O -X3.5i >> $ps 23done 24 25gmt psxy -R -J -O -T >> $ps 26gmt psconvert -A -P -Tf $ps 27rm gmt.conf gmt.history 28