# 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.))