# sigma2 stress
if "-S2" in sArg:
fSigma2 = float(sArg[3:])
# theta, in degrees
if "-T" in sArg:
fTheta = float(sArg[2:])
# filename is mandatory
if not bFilename:
sys.exit('Usage: python3 fpqSlipTendency.py -I<inputfilename>')
# get the x,y nodes & calculate angles from the file data
nodelist = fpq.getNodes(fn)
xnodelist = nodelist[0]
ynodelist = nodelist[1]
segangle = fpq.getSegAngles(xnodelist, ynodelist)
# get the cmocean 'thermal' colour map, scaled from 0-100%
start = 0.0
stop = 1.0
number_of_lines = 101
cm_subsection = np.linspace(start, stop, number_of_lines)
colours = [cmo.cm.thermal(x) for x in cm_subsection]
# slip tendency
# calculate max. possible slip tendency for these stresses
nAlpha = np.arange(1, 180, 1)
Tsmax = max(
((fSigma1 - fSigma2) * np.sin(nAlpha * DEG2RAD) * np.cos(nAlpha * DEG2RAD))
/ (fSigma1 * np.cos(nAlpha * DEG2RAD)**2.0 +
fSigma2 * np.sin(nAlpha * DEG2RAD)**2.0))
def plot_folds(linear_filename, bbox):
import fracpaq as fpq
import matplotlib.pylab as plt
import sys
import numpy as np
linear = gpd.read_file(linear_filename)
folds = linear[linear['feature'].str.contains("Fold")]
folds = folds.dropna(subset=['geometry'])
f = open('fold.txt', 'w')
for inf, fold in folds.iterrows():
for coords in fold['geometry'].coords:
ostr = "{}\t{}\t".format(coords[0], coords[1])
f.write(ostr)
f.write('\n')
f.close()
fn = 'fold.txt'
# defaults
CM2INCHES = 0.3937
bGrid = False
bEqualArea = False
nBinWidth = 5
sColour = 'C0'
xSize = 15.0 * CM2INCHES
ySize = 15.0 * CM2INCHES
# get nodes and calculate angles from nodes
nodelist = fpq.getNodes(fn)
xnodelist = nodelist[0]
ynodelist = nodelist[1]
segangle = fpq.getSegAngles(xnodelist, ynodelist)
nSegs = len(segangle)
if (nSegs > 0):
# bin the data and find maximum per bin
nBins = int(round(360 / nBinWidth))
segangleDoubled = np.zeros(len(segangle))
segangleDoubled = np.copy(segangle)
segangleDoubled = np.concatenate(
[segangleDoubled, segangleDoubled + 180.0])
n, b = plt.histogram(segangleDoubled, nBins)
nMax = max(n)
# plot the segment angle distribution
plt.figure(figsize=(xSize, ySize))
plt.subplot(111, projection='polar')
coll = fpq.rose(segangle,
bins=nBins,
bidirectional=True,
eqarea=True,
color=sColour)
plt.xticks(np.radians(range(0, 360, 45)),
['0', '45', '90', '135', '180', '215', '270', '315'])
if (nMax > 6):
plt.rgrids(range(0, int(round(nMax * 1.1)),
int(round((nMax * 1.1) / 5))),
angle=330)
plt.ylim(0, int(round(nMax * 1.1)))
plt.title('Segment strikes, n=%i' % nSegs)
plt.savefig("folds_rose.pdf")
#folds.plot(figsize=(7,7),edgecolor='#ff0000',linewidth=0.3)
#print('Plotted %5d segments & angles' % nSegs)
#plt.savefig('folds.pdf')
plt.title('Fold axial traces')
plt.show()