def METHOD_NAME(self):
    o = PelagicEggDrift(loglevel=20)
    reader_nordic = reader_ROMS_native.Reader(o.test_data_folder() + '2Feb2016_Nordic_sigma_3d/Nordic-4km_SLEVELS_avg_00_subset2Feb2016.nc')
    reader_arctic = reader_netCDF_CF_generic.Reader(o.test_data_folder() + '2Feb2016_Nordic_sigma_3d/Arctic20_1to5Feb_2016.nc')
    ######################################################
    # Vertical interpolation is another issue to be fixed:
    reader_nordic.zlevels = reader_arctic.z
    ######################################################
    o.add_reader([reader_nordic, reader_arctic])
    o.set_config('environment:fallback:land_binary_mask', 0)
    o.set_config('environment:fallback:x_wind', 10)  # Some wind for mixing
    # Seed close to Nordic boundary
    o.seed_elements(lon=14.9, lat=71.1, radius=2000, number=100,
                    time=reader_nordic.start_time, z=0)
    o.set_config('vertical_mixing:timestep', 20)
    o.set_config('drift:advection_scheme', 'runge-kutta')
    o.run(steps=5, time_step=3600, time_step_output=3600)
    self.assertEqual(o.num_elements_active(), 100)
    self.assertEqual(o.num_elements_deactivated(), 0)
    self.assertAlmostEqual(o.elements.lat[0], 71.16, 1)
    #self.assertAlmostEqual(o.elements.z.min(), -41.92, 1)  # With past ROMS-masking
    self.assertAlmostEqual(o.elements.z.min(), -46.75, 1)
    self.assertAlmostEqual(o.elements.z.max(), -0.28, 1)
    #self.assertAlmostEqual(o.elements.lon.max(), 14.949, 2)
    self.assertAlmostEqual(o.elements.lon.max(), 14.915, 2)