def test17(self):
instance = QGmodel()
self.assertEqual(instance.parameters.position_x, 0. | units.m)
self.assertEqual(instance.parameters.position_y, 0. | units.m)
instance.parameters.position_y = 123456. | units.m
instance.parameters.position_x = 543216. | units.m
self.assertEqual(instance.parameters.position_y, 123456. | units.m)
self.assertEqual(instance.parameters.position_x, 543216. | units.m)
x = instance.grid[0, 0, 0].x
y = instance.grid[0, 0, 0].y
self.assertEqual(y, 123456. | units.m)
self.assertEqual(x, 543216. | units.m)
i, j = instance.get_index_of_position(x, y)
self.assertEqual(i, 1)
self.assertEqual(j, 1)
instance.stop()
def test12(self):
instance = QGmodel()
instance.parameters.Ly = instance.parameters.Ly / 16
instance.parameters.Lx = instance.parameters.Lx / 16
Nx, Ny, Nm = instance.grid.shape
ix = numpy.arange(Nx) + 1
iy = numpy.arange(Ny) + 1
onesx = numpy.ones(Nx)
onesy = numpy.ones(Ny)
i1 = numpy.outer(ix, onesx).flatten()
j1 = numpy.outer(onesy, iy).flatten()
x = instance.grid.x.flatten()
y = instance.grid.y.flatten()
i, j = instance.get_index_of_position(x, y)
self.assertEquals(i - i1, 0.)
self.assertEquals(j - j1, 0.)
instance.stop()
