Increasing the fault dip in test.py

[rzinke]

Hello,

I am testing out these routines using the test.py script. I wanted to model an east-striking strike-slip fault, so I modified part of the test.py file to:

#-----------------------------------
# 2. Make a dislocation
#-----------------------------------

xc = np.array([0.0])
yc = np.array([0.0])
depth = np.array([2.12133])
length = np.array([10.0])
width = np.array([6.0])
dip = np.array([90.])  # was originally 45.
strike = np.array([90.0])  # was originally 0.0

ss = np.array([1.0])
ds = np.array([0.0])
ts = np.array([0.0])

mu = 30.0e9
nu = 0.25

However, when I run the script, I get the following:

(csi) MT-109587:test rzinke$ python3 test.py 
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 0, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 1, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 2, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 3, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 4, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 5, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 6, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 7, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 8, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 9, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 10, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 11, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 12, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 13, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 14, subfault 0)
10.000000 6.000000 -0.878655 
Warning: Unphysical model (station 15, subfault 0)
10.000000 6.000000 -0.878655 
...

And the displacement and stress fields all turn out to zeros.

If it helps, I am running commit 279f14a on MacOS Catalina.

[rzinke]

I see that it has to do with the fault width, but why does going below ~7 km (depth + width) raise a warning?