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