def test__rotated_pole(self):
cs = RotatedGeogCS(grid_north_pole_latitude=75.3,
grid_north_pole_longitude=54.321,
ellipsoid=self.default_ellipsoid)
test_cube = self._make_test_cube(cs=cs)
grid_definition_template_1(test_cube, self.mock_grib)
self._check_key("latitudeOfSouthernPole", -75300000)
self._check_key("longitudeOfSouthernPole", 234321000)
self._check_key("angleOfRotation", 0)
def test__grid_points(self):
test_cube = self._make_test_cube(x_points=[1, 3, 5, 7],
y_points=[4, 9])
grid_definition_template_1(test_cube, self.mock_grib)
self._check_key("longitudeOfFirstGridPoint", 1000000)
self._check_key("longitudeOfLastGridPoint", 7000000)
self._check_key("latitudeOfFirstGridPoint", 4000000)
self._check_key("latitudeOfLastGridPoint", 9000000)
self._check_key("DxInDegrees", 2.0)
self._check_key("DyInDegrees", 5.0)
def test__shape_of_earth_spherical(self):
ellipsoid = GeogCS(1.23)
cs = RotatedGeogCS(grid_north_pole_latitude=90.0,
grid_north_pole_longitude=0.0,
ellipsoid=ellipsoid)
test_cube = self._make_test_cube(cs=cs)
grid_definition_template_1(test_cube, self.mock_grib)
self._check_key('shapeOfTheEarth', 1)
self._check_key('scaleFactorOfRadiusOfSphericalEarth', 0)
self._check_key('scaledValueOfRadiusOfSphericalEarth', 1.23)
def test__fail_rotated_pole_nonstandard_meridian(self):
cs = RotatedGeogCS(grid_north_pole_latitude=90.0,
grid_north_pole_longitude=0.0,
north_pole_grid_longitude=22.5,
ellipsoid=self.default_ellipsoid)
test_cube = self._make_test_cube(cs=cs)
with self.assertRaisesRegexp(
TranslationError,
'not yet support .* rotated prime meridian.'):
grid_definition_template_1(test_cube, self.mock_grib)
def test__shape_of_earth_flattened(self):
ellipsoid = GeogCS(semi_major_axis=1.456, semi_minor_axis=1.123)
cs = RotatedGeogCS(grid_north_pole_latitude=90.0,
grid_north_pole_longitude=0.0,
ellipsoid=ellipsoid)
test_cube = self._make_test_cube(cs=cs)
grid_definition_template_1(test_cube, self.mock_grib)
self._check_key('shapeOfTheEarth', 7)
self._check_key('scaleFactorOfEarthMajorAxis', 0)
self._check_key('scaledValueOfEarthMajorAxis', 1.456)
self._check_key('scaleFactorOfEarthMinorAxis', 0)
self._check_key('scaledValueOfEarthMinorAxis', 1.123)
def test__scanmode(self):
grid_definition_template_1(self.test_cube, self.mock_grib)
self._check_key('iScansPositively', 1)
self._check_key('jScansPositively', 1)
def test__grid_shape(self):
test_cube = self._make_test_cube(x_points=np.arange(13),
y_points=np.arange(6))
grid_definition_template_1(test_cube, self.mock_grib)
self._check_key('Ni', 13)
self._check_key('Nj', 6)
def test__template_number(self):
grid_definition_template_1(self.test_cube, self.mock_grib)
self._check_key('gridDefinitionTemplateNumber', 1)
def test__scanmode_reverse(self):
test_cube = self._make_test_cube(x_points=np.arange(7, 0, -1))
grid_definition_template_1(test_cube, self.mock_grib)
self._check_key('iScansPositively', 0)
self._check_key('jScansPositively', 1)
