#!/bin/bash
#	$Id: meca_5.sh 12114 2013-09-03 19:19:00Z fwobbe $
#
# Check gmt psmeca for plotting beach balls

ps=meca_5.ps

size=1.0c

#THOUGHTS:
#echo 2.0 1.0 0.0  90 90  90 270  0  90 4 23 0 0 Dip slip    | gmt psmeca -Sc${size} -L1 -Gblack -R -J -K -O >> $ps
#echo 2.0 1.0 0.0  90 90  90  90  0 -90 4 23 0 0 Dip slip    | gmt psmeca -Sc${size} -L1 -Gblack -R -J -K -O >> $ps
#The previous two lines don't give the same result.. Maybe there could be a warning to not use the
#same strike twice?

###THE FOLLOWING DOES NOT WORK
###echo 2.0 1.0 0.0  90  4 270 5 0 0 Dip slip    | gmt psmeca -Sa${size}  -Gblack -R -J -K -O >> $ps
#the rake should be within -180 to 180, leaving it as 270 breaks the code.. should there be a warning?


#TESTING AKI OPTION -Sa  - PURE STRIKE SLIP
# Explocion
# - not possible
# Implocion
# - not possible
# Right lateral strike slip
echo 1.0 1.0 0.0  0 90  0 5 0 0 Strike slip | gmt psmeca -Sa${size} -L1 -Gblack -R-1/6/-1/2 -JM14c -P -Y20c -B2 -K > $ps
# Rotated strike slip
echo 1.0 0.0 0.0  45 90 180 5 0 0 Strike slip 45 deg | gmt psmeca -Sa${size} -L1 -Gblack -R -J -K -O >> $ps
# Dip slip
#echo 2.0 1.0 0.0  90 90 90 5 0 0 Dip slip    | gmt psmeca -Sa${size} -L1 -Gblack -R -J -K -O >> $ps
###THE FOLLOWING DOES NOT WORK
###echo 2.0 1.0 0.0  90  4 270 5 0 0 Dip slip    | gmt psmeca -Sa${size}  -Gblack -R -J -K -O >> $ps
echo 2.0 1.0 0.0  90  0 -90 5 0 0 Dip slip    | gmt psmeca -Sa${size} -L1 -Gblack -R -J -K -O >> $ps
# Dip slip
#echo 2.0 0.0 0.0   0 90  90 5 0 0 Dip slip    | gmt psmeca -Sa${size} -L1 -Gblack -R -J -K -O >> $ps
echo 2.0 0.0 0.0   0  0 -90 5 0 0 Dip slip    | gmt psmeca -Sa${size} -L1 -Gblack -R -J -K -O >> $ps
# Thrust
echo 3.0 1.0 0.0  90 45 90 5 0 0 Thrust      | gmt psmeca -Sa${size} -L1 -Gblack -R -J -K -O >> $ps
# Thrust
echo 3.0 0.0 0.0   0 45 90 5 0 0 Thrust      | gmt psmeca -Sa${size} -L1 -Gblack -R -J -K -O >> $ps
# Horizontal CLVD
# - not possible
# Horizontal CLVD
# - not possible
# Vertical CLVD
# - not possible
# Vertical CLVD
# - not possible

#TESTING HARVARD OPTION -Sc - two fault planes

# Right lateral Strike Slip

# Rotated strike Slip
# Explocion
# - not possible
# Implocion
# - not possible
# Right lateral strike slip
echo 1.0 1.0 0.0  90 90 180   0 90   0 4 23 0 0 Right Strike Slip | gmt psmeca -Sc${size} -L1 -Gblack -R -JM14c -Y-8c -P -B2 -K -O >> $ps
# Rotated strike slip
echo 1.0 0.0 0.0  45 90 180 135 90   0 4 23 0 0 Strike slip 45 deg | gmt psmeca -Sc${size} -L1 -Gblack -R -J -K -O >> $ps
# Dip slip
##The following two lines should look equal to the ones in the box above and below, but does not
echo 2.0 1.0 0.0  90 90  90 270  0  90 4 23 0 0 Dip slip    | gmt psmeca -Sc${size} -L1 -Gblack -R -J -K -O >> $ps
# Dip slip
echo 2.0 0.0 0.0   0 90  90 180  0  90 4 23 0 0 Dip slip    | gmt psmeca -Sc${size} -L1 -Gblack -R -J -K -O >> $ps
# Thrust
echo 3.0 1.0 0.0  90 45  90 270 45  90 4 23 0 0 Thrust      | gmt psmeca -Sc${size} -L1 -Gblack -R -J -K -O >> $ps
# Thrust
echo 3.0 0.0 0.0   0 45  90 180 45  90 4 23 0 0 Thrust      | gmt psmeca -Sc${size} -L1 -Gblack -R -J -K -O >> $ps
# Horizontal CLVD
# - not possible
# Horizontal CLVD
# - not possible
# Vertical CLVD
# - not possible
# Vertical CLVD
# - not possible

#TESTING HARVARD OPTION -Sm - moment tensor elements
#Following Dahlen and Tromp (1998) page 174
# Explocion
echo 0.0 1.0 0.0  1.0  1.0  1.0  0.0  0.0  0.0 23 0 0 Explosion   | gmt psmeca -Sm${size} -L1 -Gblack -R -J -Y-8c -P -B2 -K -O >> $ps
# Implocion
echo 0.0 0.0 0.0 -1.0 -1.0 -1.0  0.0  0.0  0.0 23 0 0 Implocion   | gmt psmeca -Sm${size} -L1 -Gblack -R -J -K -O >> $ps
# Right lateral strike slip
echo 1.0 1.0 0.0  0.0  0.0  0.0  0.0  0.0 -1.0 23 0 0 Strike slip | gmt psmeca -Sm${size} -L1 -Gblack -R -J -K -O >> $ps
# Rotated strike slip
echo 1.0 0.0 0.0  0.0  1.0 -1.0  0.0  0.0  0.0 23 0 0 Strike slip 45 deg | gmt psmeca -Sm${size} -L1 -Gblack -R -J -K -O >> $ps
# Dip slip
echo 2.0 1.0 0.0  0.0  0.0  0.0  1.0  0.0  0.0 23 0 0 Dip slip    | gmt psmeca -Sm${size} -L1 -Gblack -R -J -K -O >> $ps
# Dip slip
echo 2.0 0.0 0.0  0.0  0.0  0.0  0.0  1.0  0.0 23 0 0 Dip slip    | gmt psmeca -Sm${size} -L1 -Gblack -R -J -K -O >> $ps
# Thrust
echo 3.0 1.0 0.0  1.0 -1.0  0.0  0.0  0.0  0.0 23 0 0 Thrust      | gmt psmeca -Sm${size} -L1 -Gblack -R -J -K -O >> $ps
# Thrust
echo 3.0 0.0 0.0  1.0  0.0 -1.0  0.0  0.0  0.0 23 0 0 Thrust      | gmt psmeca -Sm${size} -L1 -Gblack -R -J -K -O >> $ps
# Horizontal CLVD
echo 4.0 1.0 0.0  1.0  1.0 -2.0  0.0  0.0  0.0 23 0 0 Horiz. CLVD | gmt psmeca -Sm${size} -L1 -Gblack -R -J -K -O >> $ps
# Horizontal CLVD
echo 4.0 0.0 0.0  1.0 -2.0  1.0  0.0  0.0  0.0 23 0 0 Horiz. CLVD | gmt psmeca -Sm${size} -L1 -Gblack -R -J -K -O >> $ps
# Vertical CLVD-L1 
echo 5.0 1.0 0.0 -2.0  1.0  1.0  0.0  0.0  0.0 23 0 0 Vertical CLVD | gmt psmeca -Sm${size} -L1 -Gblack -R -K -J -O >> $ps
# Vertical CLVD
echo 5.0 0.0 0.0  2.0 -1.0 -1.0  0.0  0.0  0.0 23 0 0 Vertical CLVD | gmt psmeca -Sm${size} -L1 -Gblack -R    -J -O >> $ps