# print(np.shape(new_xx), np.shape(new_yy), np.shape(stress[0, :, :]))
# plt.pcolor(new_xx, new_yy, stress[2, :, :], vmax=0.)
# plt.colorbar()
# plt.show()
# sys.exit()
p.set_het_stress(stress)
# set material properties
# Same material properties used in Tectonic model
density = 2.7
first_lames_param = 30
shear_modulus = 30
p.set_material(
fdfault.material('elastic', density, first_lames_param, shear_modulus))
# # set interface types
p.set_iftype(0, 'locked')
p.set_iftype(1, 'slipweak')
p.set_iftype(2, 'locked')
p.set_iftype(3, 'locked')
# set slip weakening parameters for fault interface-1
p.add_pert(fdfault.swparam('constant', dc=0.4, mus=0.681, mud=0.25), 1)
# p.add_pert(fdfault.swparam('boxcar', 0., x0 = 65.0, dx = 0.25, mus = -0.20, c0 = 0.),1)
# p.add_pert(fdfault.swparam('boxcar', 0., x0 = 61.0, dx = 1.0, mus = 10000., c0 = 500.),1)
# # add cohesionion to top 5 km of the fault
# cohesion = np.zeros((nby1,1))
# zer = np.zeros((nby1,1))
p.set_bounds(
(0, 1, 0),
['absorbing', 'absorbing', 'none', 'absorbing', 'absorbing', 'free'])
# turn on artificial dissipation
#p.set_cdiss(0.1)
# set material
cs = 3.464
cp = 6.
rho = 2.67
p.set_material(
fdfault.material('elastic', rho, rho * (cp**2 - 2. * cs**2), rho * cs**2))
# set interface type
p.set_iftype(0, 'slipweak')
# set slip weakening parameters
p.add_pert(fdfault.swparam('constant', 0., 0., 0., 0., 0., 0.4, 0.677, 0.525))
p.add_pert(
fdfault.swparam('boxcar', 0., 0., 20., -17.55, 2.5, 0., 10000., 0., 10.))
p.add_pert(
fdfault.swparam('boxcar', 0., -17.55, 2.5, -7.5, 7.5, 0., 10000., 0., 10.))
p.add_pert(
fdfault.swparam('boxcar', 0., 17.55, 2.5, -7.5, 7.5, 0., 10000., 0., 10.))