def METHOD_NAME(self):
    reader_osc = reader_oscillating.Reader(
            'x_sea_water_velocity', amplitude=.5,
            period_seconds=3600*6,
            zero_time=datetime.now())
    reader_global = reader_global_landmask.Reader()
    # Three different stranding options, with
    # expected final status and position
    options = ['stranding', 'previous', 'none']
    status = ['stranded', 'active', 'active']
    lons_backward = [16.157, 16.087, 16.167]
    lons_forward = [16.198, 16.115, 16.167]
    for i, option in enumerate(options):
        for direction in ['forward', 'backward']:
            if direction == 'forward':
                lons = lons_forward
                time_step = 900
            else:
                lons = lons_backward
                time_step = -900
            o = OceanDrift(loglevel=50)
            o.set_config('general:coastline_action', option)
            o.add_reader([reader_osc, reader_global])
            # Adding northwards drift
            o.set_config('environment:constant:y_sea_water_velocity', .1)
            #o.seed_elements(lon=13.35, lat=68.0, radius=0,
            o.seed_elements(lon=16.12, lat=68.5, radius=0,
                            time=reader_osc.zero_time)
            o.run(duration=timedelta(hours=22), time_step=time_step)
            #o.animation()
            #o.plot ()
            print('Testing stranding: %s, %s' % (option, direction))
            if len(o.elements) == 1:
                el = o.elements
            else:
                el = o.elements_deactivated
            self.assertEqual(o.status_categories[int(el.status)], status[i])
            self.assertIsNone(np.testing.assert_array_almost_equal(
                el.lon, lons[i], 3))