def test_transposed(self):
section = self.section_3()
section['scanningMode'] = 0b01100000
metadata = empty_metadata()
grid_definition_template_12(section, metadata)
expected = self.expected(1, 0)
self.assertEqual(metadata, expected)
def test_incompatible_grid_extent(self):
section = self.section_3()
section['X2'] += 100
metadata = empty_metadata()
with self.assertRaisesRegexp(iris.exceptions.TranslationError,
'grid'):
grid_definition_template_12(section, metadata)
def test_negative_y(self):
section = self.section_3()
section['scanningMode'] = 0b00000000
metadata = empty_metadata()
with self.assertRaisesRegexp(iris.exceptions.TranslationError,
'-y scanning'):
grid_definition_template_12(section, metadata)
def test_transposed(self):
section = self.section_3()
section['scanningMode'] = 0b01100000
metadata = empty_metadata()
grid_definition_template_12(section, metadata)
expected = self.expected(1, 0)
self.assertEqual(metadata, expected)
def test_negative_y(self):
section = self.section_3()
section['scanningMode'] = 0b00000000
metadata = empty_metadata()
with self.assertRaisesRegexp(iris.exceptions.TranslationError,
'-y scanning'):
grid_definition_template_12(section, metadata)
def test_spherical(self):
section = self.section_3()
section['shapeOfTheEarth'] = 0
metadata = empty_metadata()
grid_definition_template_12(section, metadata)
expected = self.expected(0, 1)
cs = expected['dim_coords_and_dims'][0][0].coord_system
cs.ellipsoid = iris.coord_systems.GeogCS(6367470)
self.assertEqual(metadata, expected)
def test_spherical(self):
section = self.section_3()
section['shapeOfTheEarth'] = 0
metadata = empty_metadata()
grid_definition_template_12(section, metadata)
expected = self.expected(0, 1)
cs = expected['dim_coords_and_dims'][0][0].coord_system
cs.ellipsoid = iris.coord_systems.GeogCS(6367470)
self.assertEqual(metadata, expected)
def test_scale_workaround(self):
section = self.section_3()
section['scaleFactorAtReferencePoint'] = 1065346526
metadata = empty_metadata()
grid_definition_template_12(section, metadata)
expected = self.expected(0, 1)
# A float32 can't hold exactly the same value.
cs = expected['dim_coords_and_dims'][0][0].coord_system
cs.scale_factor_at_central_meridian = 0.9996012449264526
self.assertEqual(metadata, expected)
def test_scale_workaround(self):
section = self.section_3()
section['scaleFactorAtReferencePoint'] = 1065346526
metadata = empty_metadata()
grid_definition_template_12(section, metadata)
expected = self.expected(0, 1)
# A float32 can't hold exactly the same value.
cs = expected['dim_coords_and_dims'][0][0].coord_system
cs.scale_factor_at_central_meridian = 0.9996012449264526
self.assertEqual(metadata, expected)
def test_di_tolerance(self):
# Even though Ni * Di doesn't exactly match X1 to X2 it should
# be close enough to allow the translation.
section = self.section_3()
section['X2'] += 1
metadata = empty_metadata()
grid_definition_template_12(section, metadata)
expected = self.expected(0, 1)
x = expected['dim_coords_and_dims'][1][0]
x.points = np.linspace(293000, 299000.01, 4)
self.assertEqual(metadata, expected)
def test_di_tolerance(self):
# Even though Ni * Di doesn't exactly match X1 to X2 it should
# be close enough to allow the translation.
section = self.section_3()
section['X2'] += 1
metadata = empty_metadata()
grid_definition_template_12(section, metadata)
expected = self.expected(0, 1)
x = expected['dim_coords_and_dims'][1][0]
x.points = np.linspace(293000, 299000.01, 4)
self.assertEqual(metadata, expected)
def test(self):
section = self.section_3()
metadata = empty_metadata()
grid_definition_template_12(section, metadata)
expected = self.expected(0, 1)
self.assertEqual(metadata, expected)
def test_incompatible_grid_extent(self):
section = self.section_3()
section['X2'] += 100
metadata = empty_metadata()
with self.assertRaisesRegexp(iris.exceptions.TranslationError, 'grid'):
grid_definition_template_12(section, metadata)
def test(self):
section = self.section_3()
metadata = empty_metadata()
grid_definition_template_12(section, metadata)
expected = self.expected(0, 1)
self.assertEqual(metadata, expected)
def test_negative_x(self):
section = self.section_3()
section["scanningMode"] = 0b11000000
metadata = empty_metadata()
with self.assertRaisesRegexp(iris.exceptions.TranslationError, "-x scanning"):
grid_definition_template_12(section, metadata)
