def constant_test():
"Model Constant"
# compute mass flux
melt_rate = config.get_number("frontal_melt.constant.melt_rate", "m second-1")
grid = create_grid()
geometry = create_geometry(grid)
inputs = PISM.FrontalMeltInputs()
water_flux = PISM.IceModelVec2S(grid, "water_flux", PISM.WITHOUT_GHOSTS)
water_flux.set(0.0)
inputs.geometry = geometry
inputs.subglacial_water_flux = water_flux
model = PISM.FrontalMeltConstant(grid)
model.init(geometry)
model.update(inputs, 0, 1)
check_model(model, melt_rate)
assert model.max_timestep(0).infinite()
