# to increase model resolution, increase refine (refine = 2 doubles resolution)
refine = 1
nt = 7000*refine # Number of timesteps
nx1 = 600*refine+1 # number of grid point in left block
nx2 = 600*refine+1 # number of grid point in right block
nx = nx1+nx2 # total for x coordinates
nby0 = 200*refine+1 # for 000 and 100 blocks
fault_surface = 51.
fault_bottom = 42.
surf_height = 10.
[x1, y1, x2, y2, block000_length, block100_length] = read_surface_data(nx1, nx2, fault_bottom,
fault_surface, surf_height)
p = fdfault.problem('normalfault')
# set rk and fd order
p.set_rkorder(4)
p.set_sbporder(4)
# set time step info
p.set_nt(nt)
p.set_cfl(0.3)
p.set_ninfo(100*refine)
# set number of blocks and coordinate information
p.set_nblocks((2,1,1))
p.set_nx_block(([nx1, nx2], [nby0], [1]))
plt.rc('font', **font)
plt.rc('text', usetex=True)
#-------- Grid set up and import surface data ----------
refine = mesh.refine
nt = mesh.nt * refine # Number of timesteps
nx1 = mesh.nx1 * refine + 1 # number of grid point in left block
nx2 = mesh.nx2 * refine + 1 # number of grid point in right block
nx = nx1 + nx2 # total for x coordinates
nby0 = mesh.nby0 * refine + 1 # for 000 and 100 blocks
nby1 = mesh.nby1 * refine + 1 # for 000 and 100 blocks
ny = nby0 + nby1
[x1, y1, x2, y2, block010_length,
block110_length] = read_surface_data(nx1, nx2)
import pickle
mydict = {'x1': x1, 'x2': x2, 'y1': y1 - 70, 'y2': y2 - 70}
output = open('surface.pkl', 'wb')
surf_height = max([max(y1), max(y2)])
arg = str(50)
problem_name = 'longterm_' + arg
p = fdfault.problem(problem_name)
# set rk and fd order
p.set_rkorder(4)
p.set_sbporder(4)
from get_stress import lithostatic_stress import mesh as mesh import matplotlib.pyplot as plt #-------- Grid set up and import surface data ---------- refine = mesh.refine nt = mesh.nt*refine # Number of timesteps nx1 = mesh.nx1*refine+1 # number of grid point in left block nx2 = mesh.nx2*refine+1 # number of grid point in right block nx = nx1+nx2 # total for x coordinates nby0 = mesh.nby0*refine+1 # for 000 and 100 blocks nby1 = mesh.nby1*refine+1 # for 000 and 100 blocks ny = nby0 + nby1 [x1, y1, x2, y2, block010_length, block110_length] = read_surface_data(nx1, nx2) surf_height = max([max(y1), max(y2)]) problem_name = 'litho_deep_initial' p = fdfault.problem(problem_name) # set rk and fd order p.set_rkorder(4) p.set_sbporder(4) # set time step info p.set_nt(nt) p.set_cfl(0.3) p.set_ninfo(100*refine) # set number of blocks and coordinate information