def METHOD_NAME(at, name, unit_str):
    grid = RasterModelGrid((4, 4))
    units(grid, name, at=at, units=unit_str)
    assert grid[at][name] == pytest.approx(0.0)
    assert grid.field_units(at, name) == unit_str
    assert grid[at].units[name] == unit_str