def test_2d_point_fields(self):
# Create a field that looks like this,
#
# (0) --- (1) --- (2)
# | | |
# | 0 | 1 |
# | | |
# (3) --- (4) --- (5)
g = RasterField((2, 3), (1, 2), (0, 0), indexing="ij")
data = np.arange(6)
g.add_field("point_var_1", data, centering="point")
g.add_field("point_var_2", data * 2, centering="point")
data = np.arange(2)
g.add_field("cell_var", data, centering="zonal")
fields = g.get_point_fields()
self.assertEqual(len(fields), 2)
self.assertTrue("point_var_1" in fields)
self.assertTrue("point_var_2" in fields)
self.assert_field_values(g, "point_var_1", [0, 1, 2, 3, 4, 5])
self.assert_field_values(g, "point_var_2", [0, 2, 4, 6, 8, 10])
def test_2d_point_fields(self):
# Create a field that looks like this,
#
# (0) --- (1) --- (2)
# | | |
# | 0 | 1 |
# | | |
# (3) --- (4) --- (5)
g = RasterField((2, 3), (1, 2), (0, 0), indexing="ij")
data = np.arange(6)
g.add_field("point_var_1", data, centering="point")
g.add_field("point_var_2", data * 2, centering="point")
data = np.arange(2)
g.add_field("cell_var", data, centering="zonal")
fields = g.get_point_fields()
self.assertEqual(len(fields), 2)
self.assertTrue("point_var_1" in fields)
self.assertTrue("point_var_2" in fields)
self.assert_field_values(g, "point_var_1", [0, 1, 2, 3, 4, 5])
self.assert_field_values(g, "point_var_2", [0, 2, 4, 6, 8, 10])
def test_2d_point_fields (self):
# Create a field that looks like this,
#
# (0) --- (1) --- (2)
# | | |
# | 0 | 1 |
# | | |
# (3) --- (4) --- (5)
g = RasterField ((2,3), (1,2), (0, 0), indexing='ij')
data = np.arange (6)
g.add_field ('point_var_1', data, centering='point')
g.add_field ('point_var_2', data*2, centering='point')
data = np.arange (2)
g.add_field ('cell_var', data, centering='zonal')
fields = g.get_point_fields ()
self.assertEqual (len (fields), 2)
self.assertTrue (fields.has_key ('point_var_1'))
self.assertTrue (fields.has_key ('point_var_2'))
self.assert_field_values (g, 'point_var_1', [0, 1, 2, 3, 4, 5])
self.assert_field_values (g, 'point_var_2', [0, 2, 4, 6, 8, 10])