def test_prepare_subset_sdim(self):
for subset_sdim in self.get_subset_sdim():
for ss, k in self:
try:
prepared = ss._prepare_subset_sdim_(subset_sdim)
# check that a deepcopy has occurred
self.assertFalse(np.may_share_memory(prepared.uid, subset_sdim.uid))
except KeyError:
# the target has a rotated pole coordinate system. transformations to rotated pole for the subset
# geometry is not supported.
if isinstance(ss.sdim.crs, CFRotatedPole):
continue
else:
raise
self.assertEqual(prepared.crs, ss.sdim.crs)
# test nebraska against an unwrapped dataset specifically
nebraska = SpatialDimension.from_records([self.nebraska])
field = self.test_data.get_rd('cancm4_tas').get()
ss = SpatialSubsetOperation(field)
prepared = ss._prepare_subset_sdim_(nebraska)
self.assertEqual(prepared.wrapped_state, WrappableCoordinateReferenceSystem._flag_unwrapped)
