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)
