from parameters.material import wetdiabase
from core.constants import consts
from rupturotops import wave_forms
assert consts
from core.debug import _DEBUG
assert _DEBUG
import numpy as np
# Define the standard parameter data structure.
params = Data()
# What material should we use? Look in the parameters/material file
# to see the options and to see how to define a new material.
params.material = wetdiabase
# Setup the solution domain, first we define the cell spacings
delta_x = []
count = 2000
domain_length = 10.0
for i in range(count):
if i % 10 == 0:
delta_x.append(domain_length / count)
else:
delta_x.append(domain_length / count)
width = np.sum(delta_x)
params.delta_x = np.array(delta_x)
# plotting parameters
