def test_crop_grid():
"""Crop a grid."""
logger.info("Read grid...")
grd = Grid(EMEGFILE2)
zprop = GridProperty(EMEZFILE2, name="Zone", grid=grd)
logger.info("Read grid... done, NLAY is {}".format(grd.nlay))
logger.info("Read grid...NCOL, NROW, NLAY is {} {} {}".format(
grd.ncol, grd.nrow, grd.nlay))
grd.crop((30, 60), (20, 40), (1, 46), props=[zprop])
grd.to_file(join("TMP", "grid_cropped.roff"))
grd2 = Grid(join("TMP", "grid_cropped.roff"))
assert grd2.ncol == 31
