def run(self): self.model = PISM.Blatter(self.grid, self.Mz, self.mg_levels, self.coarsening_factor) inputs = PISM.StressBalanceInputs() inputs.geometry = self.geometry inputs.basal_yield_stress = self.yield_stress inputs.enthalpy = self.enthalpy self.model.update(inputs, True)