def initPhysics(self): config = self.config self.basal = PISM.IceBasalResistancePlasticLaw( config.get("plastic_regularization") / PISM.secpera, config.get_flag("do_pseudo_plastic_till"), config.get("pseudo_plastic_q"), config.get("pseudo_plastic_uthreshold") / PISM.secpera) if PISM.optionsIsSet("-ssa_glen"): self.ice = PISM.CustomGlenIce(self.grid.com, "", config) B_schoof = 3.7e8 # Pa s^{1/3}; hardness self.ice.setHardness(B_schoof) else: self.ice = PISM.GPBLDIce(self.grid.com, "", config) self.ice.setFromOptions() self.enthalpyconverter = PISM.EnthalpyConverter(config) if PISM.getVerbosityLevel() > 3: self.enthalpyconverter.viewConstants(PETSc.Viewer.STDOUT())