def METHOD_NAME(sample_model, output_dir):
    METHOD_NAME = sample_model['model']
    rel_start_pos = sample_model['release_start_pos']
    rel_end_pos = sample_model['release_end_pos']
    METHOD_NAME.start_time = datetime(2015, 5, 14, 0)
    METHOD_NAME.cache_enabled = True
    METHOD_NAME.uncertain = True
    N = 10  # a line of ten points
    line_pos = np.zeros((N, 3), dtype=np.float64)
    line_pos[:, 0] = np.linspace(rel_start_pos[0], rel_end_pos[0], N)
    line_pos[:, 1] = np.linspace(rel_start_pos[1], rel_end_pos[1], N)
    start_pos = (-164.01696, 72.921024, 0)
    METHOD_NAME.spills += surface_point_line_spill(1,
                                             start_position=start_pos,
                                             release_time=METHOD_NAME.start_time,
                                             end_position=start_pos)
    release = Release(custom_positions=line_pos,
                             release_time=METHOD_NAME.start_time)
    METHOD_NAME.spills += Spill(release)
    METHOD_NAME.movers += c_ice_mover
    METHOD_NAME.outputters += IceGeoJsonOutput([c_ice_mover])
    METHOD_NAME.rewind()
    return METHOD_NAME