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)