def test_age_model_runs(self):
"Check if AgeModel runs"
ice_thickness = PISM.model.createIceThicknessVec(self.grid)
u = PISM.IceModelVec3(self.grid, "u", PISM.WITHOUT_GHOSTS,
self.grid.z())
v = PISM.IceModelVec3(self.grid, "v", PISM.WITHOUT_GHOSTS,
self.grid.z())
w = PISM.IceModelVec3(self.grid, "w", PISM.WITHOUT_GHOSTS,
self.grid.z())
ice_thickness.set(4000.0)
u.set(0.0)
v.set(0.0)
w.set(0.0)
model = PISM.AgeModel(self.grid, None)
input_options = PISM.process_input_options(ctx.com, ctx.config)
model.init(input_options)
inputs = PISM.AgeModelInputs(ice_thickness, u, v, w)
dt = PISM.util.convert(1, "years", "seconds")
model.update(0, dt, inputs)
model.age().dump(self.output_file)
grid = create_dummy_grid()
ice_thickness = PISM.model.createIceThicknessVec(grid)
u = PISM.IceModelVec3()
u.create(grid, "u", PISM.WITHOUT_GHOSTS)
v = PISM.IceModelVec3()
v.create(grid, "v", PISM.WITHOUT_GHOSTS)
w = PISM.IceModelVec3()
w.create(grid, "w", PISM.WITHOUT_GHOSTS)
ice_thickness.set(4000.0)
u.set(0.0)
v.set(0.0)
w.set(0.0)
model = PISM.AgeModel(grid, None)
input_options = PISM.process_input_options(ctx.com, ctx.config)
model.init(input_options)
inputs = PISM.AgeModelInputs(ice_thickness, u, v, w)
dt = PISM.util.convert(1, "years", "seconds")
model.update(0, dt, inputs)
model.age().dump("age.nc")
