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")