Open luaschwarz opened 1 year ago
This is currently done instantiating in world.neighbourhood via lpjml.grid.get_neighbourhood(), which looks like the following:
earth.neighbourhood
# <xarray.LPJmLData 'neighbourhood' (cell: 9, neighbour: 8)>
# array([[27411, -9999, -9999, -9999, -9999, -9999, -9999, -9999],
# [27410, 27412, -9999, -9999, -9999, -9999, -9999, -9999],
# [27411, 27413, -9999, -9999, -9999, -9999, -9999, -9999],
# [27412, 27414, -9999, -9999, -9999, -9999, -9999, -9999],
# [27413, 27415, -9999, -9999, -9999, -9999, -9999, -9999],
# [27414, -9999, -9999, -9999, -9999, -9999, -9999, -9999],
# [27417, -9999, -9999, -9999, -9999, -9999, -9999, -9999],
# [27416, 27418, -9999, -9999, -9999, -9999, -9999, -9999],
# [27417, -9999, -9999, -9999, -9999, -9999, -9999, -9999]])
# Coordinates:
# * cell (cell) int64 27410 27411 27412 27413 ... 27415 27416 27417 27418
# longitude (cell) float64 7.75 7.75 7.75 7.75 7.75 7.75 7.75 7.75 7.75
# latitude (cell) float64 51.25 51.75 52.25 52.75 ... 58.25 58.75 59.25
# * neighbour (neighbour) int64 0 1 2 3 4 5 6 7