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)
self.ice = PISM.CustomGlenIce(self.grid.com, "", config)
self.enthalpyconverter = PISM.EnthalpyConverter(config)
def initPhysics(self):
secpera = PISM.secpera
self.ice = PISM.CustomGlenIce(self.grid.com,"",config)
self.ice.setHardness(DEFAULT_CONSTANT_HARDNESS_FOR_SSA);
self.basal = PISM.IceBasalResistancePlasticLaw( config.get("plastic_regularization") / secpera,
config.get_flag("do_pseudo_plastic_till"),
config.get("pseudo_plastic_q"),
config.get("pseudo_plastic_uthreshold") / secpera)
self.basal.printInfo(1,self.grid.com)
self.enthalpyconverter = PISM.EnthalpyConverter(config)
def initPhysics(self):
config = self.config
linear_q = 1.
self.basal = PISM.IceBasalResistancePlasticLaw(
config.get("plastic_regularization") / PISM.secpera,
True, # do not force pure plastic
linear_q,
config.get("pseudo_plastic_uthreshold") / PISM.secpera)
# The following are irrelevant because we force linear rheology
# and don't use the enthalpyconverter
self.ice = PISM.CustomGlenIce(self.grid.com, "", config)
self.enthalpyconverter = PISM.EnthalpyConverter(config)
def initPhysics(self):
config = self.config
do_pseudo_plastic = True
plastic_q = 0.
self.basal = PISM.IceBasalResistancePlasticLaw(
config.get("plastic_regularization") / PISM.secpera,
do_pseudo_plastic, plastic_q,
config.get("pseudo_plastic_uthreshold") / PISM.secpera)
self.ice = PISM.CustomGlenIce(self.grid.com, "", config)
self.ice.setHardness(B_schoof)
# irrelevant
self.enthalpyconverter = PISM.EnthalpyConverter(config)
def initPhysics(self):
config = self.config
#// The following is irrelevant because tauc=0
linear_q = 1.
self.basal = PISM.IceBasalResistancePlasticLaw(
config.get("plastic_regularization") / PISM.secpera,
True, #// do not force a pure-plastic law
linear_q,
config.get("pseudo_plastic_uthreshold") / PISM.secpera)
#// Use constant hardness
self.ice = PISM.CustomGlenIce(self.grid.com, "", config)
self.ice.setHardness(B0)
self.ice.setExponent(glen_n)
#// Enthalpy converter is irrelevant for this test.
self.enthalpyconverter = PISM.EnthalpyConverter(config)
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())
