def METHOD_NAME(self):
    '''
    test weather_elements removes amount of mass specified by timeseries
    '''
    time_step = 900
    total = self.sc.spills[0].get_mass()
    model_time = self.sc.spills[0].release_time
    self.prepare_test_objs()
    self.b.prepare_for_model_run(self.sc)
    print(f"{self.b.water}")
    from gnome.environment import  Water
    self.b.water = Water()
    assert self.sc.mass_balance['observed_beached'] == 0.0
    while (model_time < self.b.active_range[1] + timedelta(seconds=time_step)):
        amt = self.sc.mass_balance['observed_beached']
        self.release_elements(time_step, model_time, self.environment)
        self.step(self.b, time_step, model_time)
        if not self.b.active:
            assert self.sc.mass_balance['observed_beached'] == amt
        else:
            # check total amount removed at each timestep
            assert self.sc.mass_balance['observed_beached'] > amt
        model_time += timedelta(seconds=time_step)
        # check - useful for debugging issues with recursion
        assert np.isclose(total,
                          self.sc.mass_balance['observed_beached'] + self.sc['mass'].sum())
    # following should finally hold true for entire run
    assert np.allclose(total,
                       self.sc.mass_balance['observed_beached'] +
                       self.sc['mass'].sum(), atol=1e-6)
    # volume units
    total_mass = self.b._get_mass(self.substance,
                                  self.b.timeseries['value'].sum(),
                                  self.b.units)
    assert np.isclose(self.sc.mass_balance['observed_beached'],
                      total_mass)