def __init__(self,
experiment,
model=1,
mode=1,
Mx=None,
Mz=15,
semianalytic=True,
initials="ABC",
executable=None):
self.model = model
self.mode = mode
self.experiment = experiment
self.initials = initials
self.semianalytic = semianalytic
if executable:
self.executable = executable
if mode == 3:
self.Mx = Mx
else:
self.Mx = 2 * MISMIP.N(self.mode) + 1
self.My = 3
if self.experiment == "2b":
self.Lz = 7000
else:
self.Lz = 6000
def x(mismip_mode, N=None):
if mismip_mode in (1, 2):
return np.linspace(-MISMIP.L(), MISMIP.L(),
2 * MISMIP.N(mismip_mode) + 1)
return np.linspace(-MISMIP.L(), MISMIP.L(), N)